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Abstract 

We present the new Coherent Exclusive Exponentiation (CEEX), the older Ex- 
clusive Exponentiation (EEX) and the semi-analytical Inclusive Exponentiation 
(IEX) for the process e + e~ — > ff + 717, / = /j, r, d, u, s, c, b with validity for centre 
of mass energies from r lepton threshold to ITeV, that is for LEP1, LEP2, SLC, 
future Linear Colliders, 6, c, r-factories etc. They are based on Yennie-Frautschi- 
Suura exponentiation. In CEEX effects due to photon emission from initial beams 
and outgoing fermions are calculated in QED up to second-order, including all inter- 
ference effects. Electroweak corrections are included in first-order, at the amplitude 
level. Beams can be polarized longitudinally and transversely, and all spin correla- 
tions are incorporated in an exact manner. EEX is more primitive, lacks initial-final 
interferences, but it is valuable for testing the newer CEEX. IEX provides us set 
of a sophisticated semi-analytical formulas for the total cross section and selected 
inclusive distributions which are mainly used for cross-checks of the MC results. 
We analyse numerical results at the Z-peak 189 GeV and 500 GeV for simple kine- 
matical cuts (comparisons with IEX) and for realistic experimental cuts. Physical 
precision and technical precision are determined for the total cross section and for 
the charge asymmetry. 
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1 Introduction 



At the end of LEP2 operation the total cross section for the process e~e + — > ff + wy 
will have to be calculated with the precision 0.2% — 1%, depending on event selection. 
The arbitrary differential distributions have to be also calculated with the corresponding 
precision. In future linear colliders (LC's) the precision requirement can be even more 
demanding. This is especially true for high luminosity linear colliders, like in the case of 
TESLA. The above new requirements necessitate development of the new calculational 
framework for the QED corrections and the construction of new dedicated MC programs. 
The present work is a part of an effort in this direction. 

The main limiting factor preventing us from getting more precise theoretical predic- 
tions for the e~e + —* ff+n'j process is higher-order QED radiative corrections (QED part 
of electroweak Standard Model). In order to achieve the 0.2% precision tag, the virtual 
corrections have to be calculated up to two-three loops and the multiple bremsstrahlung 
up to two-three hard photons, integrating exactly the multiphoton phase-space for the 
arbitrary event selection (phase-space limits). 

For any realistic kinematical cuts, one cannot get the precise theoretical predictions 
for e~e + — > / '/ + w-f at the above ambitious precision level without Monte Carlo (MC) 
event generators. It is therefore mandatory to formulate perturbative Standard Model 
(SM) calculations in a formulation friendly to their use within the Monte Carlo event 
generator. 

Let us stress that the Monte Carlo method is for us nothing more and nothing less 
than the numerical integration over the Lorentz invariant phase-space. It is therefore an 
exercise in the applied mathematics. In the present work we shall not, however, elaborate 
on the methods of the Monte Carlo phase-space integration and construction of the Monte 
Carlo event generator. This is delegated to ref. jl| which describes the new Monte Carlo 
event generator JCJC in which the matrix element of the present paper is implemented and 
all numerical results presented here are calculated using the latest version 4.13 of K.JC. 

In the present work we concentrate on the definition and construction of the matrix 
element for the process e~e + — > // within Standard Model. We shall especially address 
the problem of the higher-order QED corrections. This work is a continuation of two 
recent papers |2| and ||. 

1.1 Two types of QED matrix element and exponentiations 

In the JCJC Monte Carlo and in this paper we use two types of matrix element with two 
types of exponentiation: exclusive exponentiation nicknamed EEX and coherent exclu- 
sive exponentiation referred to as CEEX. Both are termed as "exclusive" as opposed to 
"inclusive", see also the discussion in (3J|. The exclusivity means that the procedure of 
exponentiation, that is summing up the infrared (IR) real and virtual contribution, within 
the standard perturbative scheme of quantum field theory, is done at the level of the fully 
differential (multiphoton) cross section, or even better, at the level of the scattering ma- 
trix element (spin amplitudes), before any phase-space integration over photon momenta 
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is done. 

The other "inclusive" exponentiation is an ad hoc procedure of summing up IR correc- 
tions after phase-space integration over photon momenta, that is for inclusive distributions. 
In spite of its weak theoretical basis the inclusive exponentiation is very commonly done 
routinely in all semi-analytical approaches like that in ref. ||. In Section |5.1| we shall 
come back to inclusive exponentiation and show how to justify it theoretically. 

The two exclusive exponentiations EEX and CEEX are well suited for the fully exclu- 
sive Monte Carlo event generators in which four-momenta of all final-state particles are 
available. Historically EEX was formulated for the first time in ref. for the initial-state 
radiation (ISR) and an improved version was presented in ref. It follows very closely 
the Yennie-Frautschi-Suura (YFS) exponentiation of the classical ref. ||. The extension 
of EEX to the final-state radiation (FSR) was done shortly thereafter ||, but it was 
actually never fully published. The computer program YFS3, in which EEX for FSR 
was implemented, was incorporated in KORALZ [|10J and some numerical results were 
published in ||, without actually giving details of the QED matrix element. The present 
work gives in fact the first full account of the EEX matrix element for ISR and FSR for 
the process e~e + — > ff + nj, f ^ e. This is to be contrasted with the situation for small 
angle Bhabha scattering (the well-known LEP-SLC luminosity process) where the EEX 
type matrix element was fully documented in refs. JllHIII]- 

CEEX is a new version of the exclusive exponentiation, generally more efficient for 
calculations beyond first-order, facilitating inclusion of full spin polarization, narrow reso- 
nances and any kind of interferences. Its first version, limited to first-order, was presented 
in ref. 0. In the present work we extend it to (still incomplete) second-order. 

Let us characterize briefly the main features of EEX and CEEX. EEX is formulated 
in terms of spin summed/averaged differential distributions, this is the source of some 
advantages and disadvantages which may be summarized as follows: 

• Differential distributions in practice are given analytically in terms of Mandelstam 
variables and scattering angles, they are therefore easy to inspect by human eye 
and to check correctness of certain important limits like leading-logarithmic and 
soft limits. 

• Analytical representation of the differential distributions allows for analytical phase- 
space integration and development of the semi-analytical formulas, which are useful 
for cross-check with the MC results. 

• Spin effects are difficult to add already at 0{a 1 ), because one is forced to calculate 
radiative corrections to spin density matrices, not an easy task. 

• Squaring sums of spin amplitudes from groups of Feynman diagrams leads to many 
interference terms which in the exponentiation procedure are handled analytically 
and individually. Because of that interference terms can be dealt with efficiently in 
EEX only for simple processes dominated by a small number of Feynman diagrams 
and only up to first-order. 
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CEEX is formulated in terms of spin amplitudes and this is also the source of some 
advantages and disadvantages: 

• Differential distributions are calculated out of spin amplitudes numerically - spin 
amplitudes are generally simpler /smaller objects, especially beyond Ola 1 ). 

• Since an analytical representation for differential distributions is not available semi- 
analytical integration over the phase-space is practically impossible. 

• Spin effects are added relatively easily, during numerical evaluation of the differential 
distributions out of spin amplitudes. Adding higher-order corrections does no make 
the treatment of spin polarization more difficult. 

• Inclusion of all kind of interference effects (among real photon emissions, many 
Feynman diagrams etc.) comes almost for free - it is done numerically in the 
process of summing and squaring various contributions to spin amplitudes. 

As we see CEEX has many advantages over EEX, so why do we keep EEX? There are 
important reasons: 

• Generally, CEEX is a relatively new invention, the older and more primitive but 
well established EEX is a useful reference for numerical tests of CEEX. 

• EEX is better suited for semi-analytical integration over the phase-space, and can 
be tested with these semi-analytical results. 

• In the present /C/C MC the 0(a 3 ) leading logarithmic corrections are available for 
EEX and are not yet available for CEEX. 

Summarizing, we see that it make sense to keep EEX as a backup solution even if we 
already rely on CEEX as a default and leading solution. 

1.2 Notation, terminology 

It is useful to introduce certain notation and terminology already at this stage. In partic- 
ular, the most common perturbative calculation (no exponentiation) is "order-by-order" . 
That is all terms beyond a certain order are set to zero. In Fig. |I| that means we end at 
certain row - at 0(a 2 ) we include the first three rows. Exponentiation is blurring this 
picture because a certain class of terms is summed up to infinite order and the meaning 
of the r — th order exponentiation is that we truncate to 0(a r ) the infrared (IR) finite 
components, the so-called /3's. On the other hand, in the leading-logarithmic approxi- 
mation the focus is on summing up first the contributions like a n L n and later a n L n_1 , 
that is in Fig. p] we sum up in column-wise order, neglecting terms far away from the 
first column which represents the so-called LL-approximation. Taking the actual value 
of a I < pi ~ 1/400 and of the big logarithm L = ln(s/m^) ~ 10, we discover quickly 
that in Fig. [I] the limiting line following the numerical importance of the terms is neither 
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Figure 1: QED perturbative leading and subleading corrections. Rows represent corrections in 
consecutive perturbative orders - first row is Born contribution. First column represents leading 
logarithmic (LL) approximation and second column depicts next-to-leading (NLL) approxima- 
tion. In the figure terms selected for (a) second and (b) third-order pragmatic expansion are 
limited with help of additional line. 



row-wise nor column-wise but diagonal- wise. This is why we shall often use 0(a r ) prag 
r — 1, 2, 3 approximation, depicted also in Fig. ^ in which we use (exponentiated or not) 
0(a r ) calculation in which we use incomplete sub-leading terms, in the sense of the LL 
approximation. Note that for the LL approximation we shall never use the strict collinear 
(zero pt) approximation. The LL approximation will be done at the level of the differen- 
tial distributions (or spin amplitudes) without forcing pt = on photons. Just to give 
a rough idea, the precision level of order 0.5 — 1% corresponds to 0(a 1 ) pra g, 0.1 — 0.5% 
to 0(a 2 ) prag and going below 0.05% will require 0(a 3 ) prag . The above is true for the 
exponentiated calculation. Lack of exponentiation makes the calculation less precise by a 
factor 2 — 5. The pure non-logarithmic terms of order 0(a 2 ) are negligible (< 10 -5 ) for 
any foreseeable practical application. 



1.3 Outline 

The outline of the paper is the following. In Section 2 we describe in detail the SM/QED 
matrix element for the exclusive exponentiation (EEX) based on the Yennie-Frautschi- 
Suura (YFS) work of ref. ||, that is of the type of matrix element defined for the first time 
in ref. ||. In Section 2 we describe the new second-order matrix element with coherent 
exclusive exponentiation (CEEX), which is the default matrix element in /C/C MC. Its 
first-order variant was given in ||, and is also defined here for the sake of completeness. 
In Section 3 we elaborate on how do we combine the electroweak corrections of refs. |5|,[14| 
with the QED corrections within EEX and CEEX. In Section 4 we discuss the differences 
between EEX and CEEX. In Section 5 we integrate analytically over the phase-space 
for the EEX matrix element in the case of very simple kinematical cuts. The resulting 
analytical results are used in Section 6 where numerical results from /C/C MC are presented. 
The most important task in Section 6 is, however, the determination of the physical 
and technical precision for the total cross section and charge asymmetry at the Z-peak, 
LEP2 and 500 GeV. In particular we discuss the contribution from the initial-final state 
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Figure 2: Kinematics of the process with multiple photon emission from the initial- and 
final-fermions in the annihilation process. 



interference (IFI) which is included in our new CEEX matrix element (IFI is neglected 
in EEX). In the last Section 7 we summarize our work. In Appendix A we define the 
Weyl-spinor techniques used in construction of CEEX multi-photon spin amplitudes. 



2 Amplitudes for Exclusive Exponentiation 

As it was already indicated, the role of the EEX matrix element described in this section 
is to provide a testing environment for the new more sophisticated matrix element of the 
CEEX class, which will be defined in the next section. 

The kinematics of the process e~e + — > f f + is depicted in fig. |2|. In the case of the 
EEX matrix element presented here we neglect the initial-final state interference (IFI). 
Consequently, we are allowed in the following to distinguish among photons emitted from 
the initial- and final-state fermions. The four-momentum 

n n' 

X = pi+p 2 -^2kj = q 1 + q 2 + ^k' l (1) 

j=i i=i 

of the s-channel virtual boson Z + 7 is then well defined. Let us denote the rest frame of 
X as XMS. 
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2.1 Master formula 

Denoting Lorentz invariant phase-space by 



d n Ups(P- Pl , P2 , ...,p n ) = fl ^ 6^ (p - f>) 

3=1 Pj V 3=1 J 



(2) 



we define for the process e (pi) + e + (p 2 ) — ¥ fill) + fioz) + n l(kj) + n'j(k[) the 0(a r ) 
total cross-section 



^eex = ^T^Tj / rfn+n ' Li P s Oi + P2; gi, 92, fei..., fe», M-, O PeLo r = 0, 1, 2, 3 

n=0 n'=0 ' ' 

(3) 

in terms of the fully differential multiphoton distribution 

n n' ( 

5=1 Z=l I 

+ ^ (3ii(X J pi J p 2j quq2,k j ) + ^ ^ Pi»P2» ft' ft' fc 



^ S F (k) 



~ (2) ~(2) 

) Pi> P2, 9i; <?2, kj, kk) -^-^ ^ff^X, Pi,P2,qi,q2,k,k 



y- l~> 2 nV^,Pi,P2, vi) ^, fi-j, fcfc; 
n>^>i Siik^Sjih) n/ J^ m >i SF(ki)S F (k m ) 



^tr" " ny ^ ei S I (kj)S I (k k )S I (ki 



(4) 



Let us explain the notation and physics content in the above expression. The YFS soft 
factors for real photons emitted from the initial- and final-state fermions read 

where electric charges of the electron and fermion / are Q e and Qf. The F-function in 
the exponential YFS form factor is defined as in ref. |J: 

Y f (n,p,p) =2Q 2 f aB(n,p,p) + 2Q 2 f am(p,p) 



2Q 2 a$t 



cftk i ( 2p — k 2p — k 



k 2 (2vr) 3 \2kp - k 2 2kp - k 2 
8 



The above form factor is infrared-finite and depends explicitly on the soft photon do- 
mains Q = Qi,Qp which includes (surrounds) the IR divergence point k = 0. We define 
Q(Q; k) = 1 for k £ and G(f2; k) = for k Q. Contributions from the real photons 
inside Q are summed to infinite-order and combined with the analogous virtual contribu- 
tions forming the exponential YFS form factor. In the Monte Carlo we generate photons 
k ^ Q characterized by the function Q(Q,k) = 1 — Q(Q,k). We require, as usual, that 
Qi and Qf are small enough (they can be chosen arbitrarily small) such that the total 
cross section as defined in eq. (|j) and any other physically meaningful observable do not 
depend on the actual choice of them, i.e. fij^ are dummy parameters in the calculation! 
If we neglect the initial-final state interference then we may choose Qj ^ Qp- Let us 
define fli with the k° < E m i n condition in the centre of the mass system of incoming e 
beams and flp with k° < E' min in the centre of the mass of the outgoing fermions //. The 
two domains differ because the Lorentz frames in which they are defined are different. 
The above choice is the easiest for the Monte Carlo generation but in the later discussion 
we shall describe in detail how do we implement the Qi = Qp option in our Monte Carlo. 
The actual YFS form factors for the above choices are well known 



(7) 



where 



Wn v , 2E min 1 2 a ( 1 7T 2 

Y e (n i; Pl ,p 2 ) = 7e In -^== + - le + Q e ~ [ - 2 + y 

v fn \ l 2Emin i 1 i n2 a ( 1 , n *\ 



_ 2 a/ 2 Pl p 2 \ o^2«/i 2 9l?2 



7 = 7e = 2Q - ln^-1 , 7 / = 2g 2 f - In ^ - 1 . (8) 
it \ mi / 7r \ mj 



2.2 Pure virtual corrections 

The perturbative QED matrix element is located in the /9-functions. The f3o function is 
"proportional" to the Born e~e + — > ff differential cross section da Bovn (s, $)/dQ and it 
contains (infrared-finite) corrections calculable order by order. According to our general 
strategy we shall calculate (3o and other /3's in the 0(a l ) prag , i = 0, 1, 2. 
The 0(a l ) prag expressions for /3q , i — 0, 1, 2 read[| 

^ r Wi,P2,gi,g 2 ) = (i + 4 r) ) (i + 4 r) ) i E y^r^ 2 ^ ^ 

k,l=l,2 

*f = 0, S? = \l, if = *f + j7 2 , if = if + ^7 S , (10) 
if = 0, 4" = \l„ if = if + ^?, if = if + ^7?, (ID 

1 It may look that we miss pure (a/7r) term in 5j F . The calculation shows || that such a non- 
logarithmic contribution is accidentally equal zero. 
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where 



#11 = Z(pi, gl), #12 = Z(pi, -<f 2 ), #21 = Z(-p 2 , Qi), ^22 = A~P2, -ql), (12) 

with all 3- vectors taken in the rest frame of the four- momentum X, that is in the frame 
XMS. 

Let us first explain the fact that instead of having a single da Born / dVL(d) with a single 
# we take an average over four In fact we could adopt one #, for example #o — 
Z(pi — P2,qi — $2) where all three-momenta are taken in XMS. The main reason for 
our apparently more complicated choice is related to the implementation of the first 
and higher-order real photon contributions in the next subsections. More precisely, it is 
well known [15|, 16| that the exact single photon ISR matrix element can be cast as a 



linear combination of the two da Born /dQ(i9k), k — 1,2 distributions. The same is true for 
FSR |16| . (Our implementation of the leading-logarithmic (LL) matrix element for 2 and 



3 real photons will also involve the linear combination of this type.) It is therefore logical 
and practical to use a similar solution already for (3q. One should also keep in mind that 
in the soft limit, when all photons are soft, then all four angles dki are identical and the 
averaging over them is a spurious operation anyway. 

The reader not familiar with exponentiation may have an even more elementary ques- 
tion: Why do we have a freedom of defining # in da Borri / dVL(d) in first place? Is this 
ambiguity dangerous? These questions are already discussed in refs. @,[l]J . The answer is 
the following: Strictly speaking the differential cross section da Born (s, d)/dVt and ^° arc 
defined within the two body phase-space. Later on they are used, however, in eq. (^) and 
in the definitions of (3^\i = 1,2, ... all over the phase-space with additional soft and/or 
hard photons. This requires some kind of extrapolation of (3q and da Born (s, #) / dfl beyond 
the two body phase-space. In ref. this extrapolation was done using manipulations 
on the four-momenta and in ref. [|lT[ it was done as an extrapolation in the Mandelstam 
variables s, t, u. Here we present another solution which is somewhere in between the 
previous two ones. What is really important, however, is that the effect due to change 
from one particular choice of extrapolation to another is always, for the entire calculation, 
a kind of "higher-order effect". For instance at 0{a l ) changing the type of extrapolation 
is an 0(a 2 ) effect! Of course, it is always wise to use some kind of "smooth" extrapolation 
which is able to minimize the higher-order effects. 

Another possible question is: Why we did not write down the second-order virtual 
correction factor in an additive way, like for instance (1 + 5^ + + 5^5^)7 We have 
opted for factorized form because it is generally known that the factorized form is closer 
to reality at higher perturbative orders. Another important reason is that the factorized 
form is easier for semi-analytical integrations over the phase-space in the next section. 

2.3 One real photon with virtual corrections 

The contributions are needed directly in eq. (|]) and the 0(a l ) prag version of /3i enters 
indirectly as a construction element in fa- They are constructed from QED distributions 
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with a single real photon emission and up to one virtual photon contribution. They are 
defined separately for initial- and final-state photons 

Pih X ^P^P2,(li,(l2,k j ) =D ( t ) I (X,p 1 ,p 2 ,qi,q 2 ,k j ) - S I (k j )^ % ~ 1) (X,p 1 ,p 2 , <?i, q 2 ), 
Pif(X,Pi,P2, 9i; 92) k'i) —Dil(X,pi,p 2} qi, q 2 , k[) — S F (k' l )f3Q~ 1 \x,p 1 ,p 2 , gi, 92), 

where 2 = 1,2. Let us define first all ingredients for the initial-state contribution. The 
single initial-state photon emission differential distribution at 0(a r ), r = 1,2, 3, with the 
eventual additional up to two-loop virtual correction from the initial- and/or final-state 
photon reads 

r=l,2 r=l,2 



where 



A^(^))(l + 4 r ~^ 



^ _ kj?2_ j p. = ^ = (1 - aj)(l - 4), 



(15) 



P1P2 P1P2 
Af(*) = 0, A«(*)^ 7 -~ 7 ln(*), 

Af\z) ee A?\z) + i 7 2 - i 7 2 ln(z) + ^ 7 2 ln 2 (z) 

We( °' 6) = 1 - Wl-a)' + (l-6)' U + a 

Again the question arises why the averaging over r in i?^ is introduced? In the case of 
just one ISR hard photon the averaging trivially disappears because d^i = $k2 and in this 
case our formula coincides with the exact 0{a x ) result, see [[I5|,[l7j, as it should. In the 
less trivial case of the presence of the additional hard photons there is an ambiguity in 
defining which is reflected in our "averaging" procedure; however, it is harmless i.e. 
the effect is of 0(a 2 )., 

It is necessary and interesting to check the soft limit. If in the presence of many 
additional photons (n > 1) we take the soft limit kj — > 0, keeping momenta of other 
photons constant, then §kr are in general all different. However, in eq. ([14]) the sums over 
da Born /dfl combine into a simple average over all four angles, as in eq. (|S|) - in fact the 
single photon distribution reduces to 

Dfi^iX.p^p^qi.q^kj) ~ S I {k j )(3^' 0) (X,pi,p 2 ,q 1 ,q 2 ) 

and therefore ^j(X,pi,p 2 , q%, q 2 , kj) is infrared-finite as required. The above argument 
shows that the extrapolations for /3 and J3\ have to be of the same type. If we have opted 
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for another extrapolation in eq. (|14D , for example without averaging, with a single angle 
$ kr — ¥ $ fc , then the extrapolation in eqs. ([|) would need to be changed appropriately. 

Another interesting limit is the collinear limit. If all (possibly hard) photons are 
collinear to initial- or final-fermions then all angles ^ sr , s, r — 1, 2 are identical and equal 
to the familiar leading-logarithmic effective scattering angle for the hard process in the 
"reduced frame" XMS. This will facilitate introduction of the higher-order LL corrections 
in the following. 

Another remark on eq. (|l4|) is in order: There are many equivalent ways, modulo 
term of 0(m 2 /s), of writing the single bremsstrahlung spin summed differential distribu- 
tion ||17| . Our choice follows the representation implemented in the Monte Carlo programs 
YFS2 0, KORALZ and MUSTRAAL [TJ, because it minimizes the machine round- 
ing errors (quite important due to the smallness of electron mass), and it is explicitly 
expressed in terms of Born differential cross sections - this feature facilitates introduction 
of electroweak corrections. 

The virtual correction (1 + Aj (a,, is taken in the leading logarithmic approxi- 
mation (sufficient for our 0(a 2 ) prag approach) and it agrees with the corresponding dis- 
tribution in ref. In the kj — > limit we have Aj(aj, j3j) — > 8j as expected, and as 
required for infrared finiteness of fifp- The other factor (1 + 5p) represents the contribu- 
tion from the simultaneous emission of the real initial and the virtual final photon. We 
again prefer the factorized form over an additive one (1 + Aj + 5 ft). 

The essential ingredients for the 0(a r ) final-state fl^p, r = 1, 2, is the single final-state 
photon emission matrix element with up to one-loop virtual initial- or final-state photon 
corrections 

r=l,2 r=l,2 J 



A 



(r-l)> 



where 



k'iQ2 > k' t qi « Vi 

vi = — , = — , vi - 



qiq2 <?1<?2 1 + vi + Q L + TJl + Q 

z, = (l -*),)(! "CO 

A ( §\z)^0, A«(z) = i 7/ +i 7/ ln(z), (1T) 

2 

rr 

W f (a,b) = 1 



m) (1 - a )(l -b) fa , b 



2qiq 2 (1 -a) 2 + (1 - b) 2 \b a 



All discussion on the ISR distribution of eq. ([14]) applies also to the above FSR distribu- 
tion. 
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2.4 Two real photons with virtual corrections 

— (2) ~(2) ~(2) 

The contributions P n ,P FF and (3 IF are related to emission of the real two initial, two final 
and one initial and one final photons correspondingly. They are genuine 0(a 2 ) objects 
because they appear in this order for the first time. For the same reason they do not 
include any virtual contributions. They are defined formally in the usual way 

Pu (X,p u p2,qi,q2, kj , Afc) = D% {X, p 1 ,p 2 ,qi,q 2 , kj , k k ) 

- Si(k j )^{ I ~ 1 '(X,p 1 ,p 2 ,q 1 ,q 2 ,k k ) - S I (k k )f3[ r I ~ 1) {X,p 1 ,p 2 ,q 1 ,q 2 ,k j ) (18) 

- S I (k j )S I (k k )^ r ~ 2) (X,p l ,p 2 ,q 1 ,q 2 ),r = 2,3, 
PP F (X , pi , p 2 , qi , q 2 , k[ , k' m ) = Dy F (X,p 1 ,p 2 ,q 1 ,q 2 ,k' l ,k' m ) 

~ S F (k' l )(3 < ff 1 \x,p l ,p 2 ,qi,q 2 ,k' m ) - S F {k' m )(3[ r F ~ 1) {X,p u p 2 ,q 1 ,q 2 ,k' l ) (19) 

- S F {kl)S F {k' m )f]i?S 2 \x,p u p 2 ,q l ,q 2 ),r = 2,3, 
Pif ( x , Pi , P2 , qi , 92 , kj, k[) = D { p F (X, pi,p 2 ,qi,q 2 , kj , k[) 

- S I {k j )(][ r f 1) (X,p 1 ,p 2 ,q 1 ,q 2 ,k' l ) - S F (k' l )^[ r f 1 \x,p 1 ,p 2 ,q 1 ,q 2 ,kj) (20) 

- S I (k j )S F (k' l )f3fl 2 (X,pi,p 2 ,q l ,q 2 ),r = 2,3. 

r 2 \ (2) 

The new objects in the above expressions are the differential distributions D IT ,D FF 

(2) 

and D IF for double bremsstrahlung. They are not taken directly from Feynman diagrams 
but they are constructed in such a way that: 

• If one photon is hard and one is soft then the single bremsstrahlung expression of 
eqs. (Ill JIB) are recovered 



• If both photons are hard and collinear then the proper LL limit, which we know 
from the double or triple convolution of the Altarelli-Parisi kernels, is also recovered. 

The resulting expressions are rather compact and the LL limit is manifest, this is not 
necessarily true for the exact double bremsstrahlung spin amplitudes (see next section). 
The method is similar to that of refs. [0,|12|- in the case of ISR we shall also include one- 
loop virtual corrections read from the triple convolution of the Altarelli-Parisi kernels, see 
ref . p . 
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Our construction in the case of the double real ISR reads as follows 



Du(X,p 1 ,p 2 ,qi,q2, h, fa 



Q 



a 



47r 2 (A;ipi)(A;ip 2 ) 4vr 2 (k 2 pi)(k 2 p 2 ) 

| e( Vl - v 2 ) (1 + Af- 1) ( Zl , z l? )) a + 6 {r - 1] 



da 



Born 



X2{ 



Born 



r=l,2 



r=l,2 



-(X 2 ,tf 



2r) 



(21) 



+ e(t; 2 - m) (1 + A^ 1) (z 2 , 221)) (1 + 5 



(r-l)i 



do- 



Born 



X2 (A 2 ;« / 1 ,/30 £ —-{x a ,^ lT ) + X302'A^) E 



d<7 



Born 



where 



a 



r=l,2 



(I.; 



r=l,2 



2r 



a 2 







02 



02 



l-a 2 I -ai 1-02 l-0i 

v i = a i + /3 i , Zi = (1 - &»)(! - Pi), Zij = (I- ati- atj)(l- fa- fa), 



X 2 (u;a,b) 



:i-^) 2 [(i- 



a 



(22) 



1 
2 



A // = °. A // 0*> z ij) = h ~ h Mzi) - h ln (%)- 



1 
6 



1 
6 



The variables &i, Pi for i-th photon are defined as in eq. (|T4T) . 

In order to understand our construction let us examine how the LL collinear limit 
is realized in the exact single bremsstrahlung matrix element of eq. ([14]) . If the photon 
carrying the fraction x\ of the beam energy is collinear, let us say, with p\ then &\ ~ x, 
Px ~ 0, all four angles are the same "$ sr — > $* and we recover immediately the correct LL 
formula 



r=1.2 



da 



Born 



r=l,2 



dn 



■(■& 



2rj 



I ^(T Born 



It is therefore natural to employ for the double emission the angular dependent Altarelli- 
Parisi (AP) factors of the kind 

^[(l-a 2 ) 2 + (l-^ 2 ) 2 ]^(l-«i) 2 + (l-A) 2 ]. 

The above formula is too simple, however, to reproduce correctly the result of the double 
convolution of the AP kernels in the case when both photons are collinear with the same 
fermion 

1 1 da Born 

2(1 + (1 - *i) 2 ) 2^ + - [^/(l - ^)] 2 ))^-(^)- 
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where x 2 = xi/(l — Xi) reflects the loss of energy in the emission cascade due to emission 
of k\. In order to match the above cascade limit we construct a better angular dependent 
AP factor as 



1 



[(1 - a,) 2 + (1 - hf] -[(1 - «' 2 ) 2 + (1 - fa) 2 }. 



The above fulfils both types of LL collinear limit, when two photons are collinear with 
a single beam or each of them follows different beam. Finally, let us reproduce the 
limit in which one photon, let us say the 1-st, is hard and the other, the 2-nd, is soft, 
V2 = «2 + $2 — > 0. In this case it is logical to split the above double bremsstrahlung 
angular dependent AP factor into two pieces 



) 2 ^[(i-«' 2 ) 
^[(i-« 2 ) 



P'2) 2 } 



and associate each one with the corresponding da Born /dfl, following eq. (|14]) . The order 
in the cascade does not matter. We simply symmetrize over the two orderings in the 
cascade - it is essentially Bose-Einstein symmetrization. 



The above construction clearly provides the correct limit D\j (ki, k 2 ) ■ 
for v\ = const and V2 — > 0. As a consequence (X, pi, p 2 , qi, q2, k%, ki 



S(k 2 )D^(k 2 ) 
is finite in the 
limit of one or both photon momenta tending to zero. 

The construction of eq. ( pT|) will be inadequate if both photons are hard and at least 
one has high transverse momentum. It reflects the fact that we do not control fully in 
EEX the second-order NLL, 0(a 2 L), contributions. However, we have known since a 
long time that the construction of the type of eq. (|21|) agrees rather well with the exact 
double bremsstrahlung matrix element calculated using spinor techniques, see [03]. For 
both photons having high transverse momenta there is only about 20% disagreement for 
the approximate and exact results (integrated over the double photon phase-space). This 
result is confirmed in the present work by the numerical comparisons of EEX and CEEX, 
where the double bremsstrahlung matrix element is exact. 

The double final-state bremsstrahlung distribution is defined/constructed in an anal- 
ogous way 

D i F ) F {X,p 1 ,p 2 ,qi,q 2 ,k' 1 ,k 2 ) = 
a a 2 9iP2 a 2gip 2 
V/ 4tt 2 (fcifcXfcij*) 4tt* (k' 2qi )(k' 2 p 2 ) fKVl -' 

^^Born 

r=l,2 



Ci)w f (fj 2 ,( 2 ) 

^lr)+X2(Ci;V2,Q 



E 



da 



Born 



=1,2 



dn 



2r) 



+ 0(v> 2 



=1,2 



dn 



J^-Born 

^) + X2(C 2 ;^,CDE-^( x ' 



=1,2 



dn 







2ri 



Ay 



(r-l). 



z 3. 



(23) 
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where 



Vi — t , ' V2 ~ -, , j Ci — 1 . > 5 C2 — -, . > • (24) 
1 + 772 1 + 771 1 + C2 1 + Ci 

The "primed" Sudakov variables are here defined differently than in the ISR case because 
the fermion momenta qi^ get affected by photon emission. Virtual corrections are absent 
because we restrict FSR to 0(a 2 )LL- The above expression is tagged with r = 2, 3 for 
0(a r ), however, FSR we implement essentially only in 0(a 2 ) and the only correction in 
0(a 3 ) is the ISR one-loop correction. 

The distribution for one photon from the initial-state and one photon from the final- 
state at 0(a r ) r = 1, 2 we construct as follows 



D ( jf(X j p 1 ,p 2 , qi,q 2 , kj, k[) = 

^Mtt 2 (k jPl )(k jP2 ) w * a *>W ^/4tt 2 (k' lPl )(k' lP2 ) 



1 - %) 2 (l - 77O 2 ^ Bor " 2 q , , (i-%) 2 (i-0) 2 ^ Bor \ v 



(™ + V o" " o W ,0 (25) 



2 2 dn y ' 2 2 dn 

1 - (i - f,,r ^ + (i-4-F(i-6)^ (xV22 , 



2 2 dfi v ' y 2 2 rffi 

(l + A^(^)) (l + Aj-^4) 



where the variables ctj, /3j,fji,^i and other components are defined as in eqs. ( JI1| , [TE| ). The 
above construction is in fact the easiest because two photons cannot be emitted in a 
cascade from one line and we fully exploit the four scattering angles in the Born differential 
cross sections. It is trivial to check that all soft and collinear limits are correct. 



2.5 Three real photons 

The differential distribution for of 3 real ISR photons is essentially obtained by the triple 
convolution of the AP kernel, for each beam separately and the the two results are com- 
bined with help of additional convolution. This exercise was done for the collinear sub- 



generator of BHLUMI |12| and we exploit here these results. Even though the collinear 
limit is of primary importance, we have to be very careful in construction of the fully 
differential triple photon distribution to preserve all soft limits: when all three photons 
are soft, when two of them are soft , and only one of them is soft. In these limits the 
three-photon differential distribution has to reproduce smoothly the previously defined 
Born, single and double bremsstrahlung distributions times the appropriate soft factor (s). 
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Otherwise we may have a problem with IR finiteness of 



f3?n( x ,Pt, Qj, fa, fa, fa) = Dfj'jiX^i, qj, ki,k 2 , k 



(r) 



Si{ki)Pu{X,pi, qj, fa, fa) - S I (k 2 )p[ 2 /(X,p i , qj, fa, fa) - S I (fa)f3[ 2 / (X,p h Qj, fa, fa) 
S I (k 1 )S I (fa)(3 < l 1 j(X,pi, qj, fa) - S^fafS^k^hx^qjM) 
~S I {fa)~S I {fa)P^{X,p i ,q J ,fa) - S I {k 1 )S I {k 2 )S I {k z )^\X, Pi ,^ 



(26) 



It is therefore not completely straightforward to turn the strictly collinear expression for 
three real photon distributions of ref. |12| into the fully differential (finite px) triple photon 
distribution which we need. As in the case of double real ISR the guiding principle is 
that (i) the hardest photon decides which of the angles is used in da Born / dQ(X 2 , fi^) and 
(ii) we have to perform Bose symmetrization, that is sum over all orderings in a cascade 
emission of several photons from one beam. For three real photons there are no virtual 
corrections. 

Our construction in the case of the triple real ISR reads as follows 



D ii (X,Pi,p 2 ,qi,q 2} k 1: k 2 ,fa) = ] [ Q 



a 



2piP2 



1=1,3 



47r 2 (kipi)(kip 2 ) 



®{vi - v 2 )Q(v 2 - v 3 ) 



da 



Born 



r=l,2 



dQ 



iX 2 ,Ar) 



X3 



0i,a' 2 ,p 2 ,al^) E 



da 



Born 



r=l,2 



dQ 



ix\$ 2r ) 



(27) 



+ remaining five permutations of (1,2,3) 



where 



X 3 («i; fl2, h, a 3 , b 3 ) = -(1 - u,) 2 [(1 - a 2 f + (1 - b 2 ) 2 } [(1 - a 3 ) 2 + (1 - b 3 ) 2 ] 



a 



«3 



03 



(28) 



" 1 - a x - a 2 ' 3 i-fa-fa 



In most cases such an approach should be enough; however, in some special cases with 
two hard photons explicitly tagged it may not be sufficient. We have programmed and 
run special tests (unpublished) relying on the up to 3 hard-photon ISR amplitudes P0| 
constructed with the methods similar to these in ref. ||21|| , in order to get additional 
confidence in the approximate real emission distrubutions presented in this Section. 
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3 Amplitudes for Coherent Exclusive Exponentiation 



The Coherent Exclusive Exponentiation (CEEX) was introduced for the first time in 
ref. @. It is deeply rooted in the Yennie-Frautschi-Suura (YFS) exponentiation ||. It 
applies in particular to processes with narrow resonances where it is related also to works 
of Greco et.al. [p2|,p3| . The exponentiation procedure, that is a reorganisation of the QED 
perturbative series such that infrared (IR) divergences are summed up to infinite-order 
is done at the spin-amplitude level for both real and virtual IR singularities. This is to 
be contrasted with traditional YFS exponentiation, on which our EEX is based, where 
isolating the real IR divergences is done for squared spin-summed spin amplitudes, that 
is for differential distributions and spin density matrices^. 

Our calculations of the spin amplitudes for fermion pair production in electron positron 
scattering is done with the help of the powerful Weyl spinor (WS) techniques. There are 
several variants of WS techniques. We have opted for the method of Kleiss and Stirling 
(KS) f24], p5|j, which we found the best suited for our CEEX. In ref. [0 the KS spinor 
technique for massless and massive fermions was reviewed and appended with the rules 
for controlling their complex phases, or equivalently, the fermion rest frame (all three 
axes) in which the fermion spin is quantised - this is a critical point if we want to control 
fully the spin density matrix of the fermions. This fermion rest frame we call the GPS 
frame and the rule for finding it we call the GPS rule. For the sake of completeness we 
include definitions of the KS spinors, photon polarization vectors, and our GPS rules in 
Appendix A. 

The very interesting feature of CEEX is that, although it is formulated entirely in 
terms of the spin-amplitudes, the IR cancellations in CEEX occur for the integrated cross 
sections (probabilities), as usual; in practice they are realised numerically. There is no 
contradiction in the above statement. In order to avoid any confusion on this point, we 
shall provide the new detailed proof of IR cancellations in CEEX scheme in one of the 
following subsections. 



3.1 Master formula 

Defining the Lorentz invariant phase-space as 

f f n n ^ 

J dUps n (P;p 1 ,p 2 ,...,p n ) = J (2?r) 4 (5(P - Y^Pi) Yl (27r) 3 2p! 



(29) 



we write the CEEX total cross section for the process 

e~{p a ) + e + { Ph ) -> f(p c ) + f( Pd ) + 7 (fc 1 ) + 7 (A; 2 ) + ... + -y(k n ),n = 0, 1, 2, oo (30) 



2 The realization of EEX for spin density matrices is an obvious generalisation of the EEX/ YFS 
exponentiation which, however, was never fully implemented in practice. 
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with polarized beams and decays of unstable final fermions being sensitive to fermion spin 
polarizations, following refs. ||, as follows: 

1 °° f 

V {r) = flux( -^ ^ J dU V S n+2(Pa+Pb]Pc,Pd, h, . . . , k n ) Pcl EX {p a , Pb, Pc Pd, h, ■ ■ ■ , K) 

(31) 

where, in the CMS flux(s) ~ 2s, 

tiLxiPaiPbiPcPd, K h,..., k n ) = L e nn-,Pa,., Pd ) Q {n) J- Yl 

0-«=±l Ai,Ai=±l 

3 

fc a fc f> a \ a \ a a \ b \ b JJl n U<ri<72 ■■■<Tn) l JJl n \\aia 2 ■ ■ ■ a n )\ X C X C A d A d c'^d > 

i ,j,l,m=0 

(32) 

and assuming domination of the s-channel exchanges, including resonances, the complete 
set of spin amplitudes for emission of n photons we define in 0(a r )cEEX r = 0,1,2 as 
follows: 

n 

^ es ■■.£)= e n s fr 1 ^ 0) (5; ^) > (33) 

pe{/,F}" i=l 

85 = E n 4f 4" G; + E w tT " J • (34) 

P e{/,F}« i=l [ j=l s \j] J 

m>(2) /pfci t„\ _ 
MV n \Xtn ■ ■ ■ a„J 

{ n m (p k i-x\ B {2) (v k ih-x\ 
Po (x, X p) + 2^ m + M~M 
j=l 5 [j] l<j<l<n 5 [j] 5 [l] 

(35) 

In the following subsections we shall explain all basic notation, then in the next section 
we shall discuss in detail the IR structure in CEEX, effectively deriving all the above 
formulas. At 0(a r ) we have to provide for functions = 0,1,.. .,r from Feynman 

diagrams, which are infrared-finite by construction ||. Their actual precise definitions 
will be given in the following. We shall define/calculate them explicitly up to 0(a 2 ). 

3.1.1 Spin notation 

In order to shorten our many formulas, we use a compact collective notations 

(P\ — (PaPbPcPd\ 

\\) - \x a x b x c xj 
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for fermion four-momenta pa, A = a,b,c,d (i.e., pi = p a ,P2 = Pb,Qi = Pc, Q2 = Pd) and 
helicities \a,A = a,b,c,d. For k = 1,2,3, a k are Pauli matrices and cr° M = 8\ tfJi is the 
unit matrix. The components £{, £%,j, k = 1, 2, 3 are the components of the conventional 
spin polarization vectors of e~ and e + respectively, defined in the so-called GPS fermion 
rest frames (see Appendix A and ref. @ for the exact definition of these frames). We 
define i° A = 1 in a non-standard way (i.e. pa -£a — m e, A = a, b). The polarimeter vectors 
he are similarly defined in the appropriate GPS rest frames of the final unstable fermions 
(pc - he = rrif,C = c,d). Note that, in general, he may depend in a non-trivial way on 



the momenta of all decay products, see refs. [g, 26-28 1 for details. We did not introduce 
polarimeter vectors for bremsstrahlung photons, i.e. we take advantage of the fact that 
luckily all high-energy experiments are completely blind to photon spin polarizations. 

3.1.2 IR regulators and YFS form- factor 

Here we introduce/explain our notation for IR integration limits for real photons in 
eqs. ([H]) and ( j32|) and in the following sections. In general, the factor G(f2) in eq. ( |3TD 
defines the infrared (IR) integration limits for all real photons. More precisely for a single 
photon, Q is the domain surrounding the IR divergence point k = 0, which is in fact 
excluded from the MC phase-space. In CEEX there is no real distinction among ISR and 
FSR photons, Q is therefore necessarily the same for all photons. We define a character- 
istic function 0(f2, k) of the IR domain Q as 0(f2, k) = 1 for k e Q and Q(Q, k) = for 
k ^ fl The characteristic function for the part of the phase-space included in the MC 
integration for a single real photon is G(f2, k) = 1 — Q(Q, k). The analogous characteristic 
function for all real photons is, of course, the following product 

n 

Q(0) = Y[Q(n,k). (36) 
i=i 

In the present calculation corresponding to the /C/C Monte Carlo program we opt for Q 
defined traditionally with the photon energy cut condition k° < E m i n . 

The YFS form factor |H1 for Q defined with the condition k° < E m - m reads 



Y(Q;p a , ...,p d ) = Q 2 e Y n (p a ,p b ) + Q 2 f Y n (p c ,p d ) 

+ QeQfY n (p a ,Pc) + QeQfYn(pb,Pd) ~ QeQ fY n (p a , p d ) - Q e QfY n (pb,Pc 

where 

Y n (p,q) =2aB(Q,p,q) + 2Q)a^B(p,q) 



(37) 



(3f 



2a5ft 



d A k i ( 2p — k 2q — k 



k 2 (2tt) 3 \2kp-k 2 2kq-k 2 



2 
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is given analytically in terms of logarithms and Spence functions. As we see, the above 
YFS form factor includes terms due to the initial- final state interference (IFI). The above 
form-factor will be derived in the following. The additional contribution to the YFS 
form-factor due to the narrow Z-resonance will be discussed in detail separately. 



3.1.3 Partitions and s-factors 

The coherent sum is taken over the set {p} = {I, F} n of all 2 n partitions - the single 
partition p is defined as a vector (p 1 , p 2 , . . . , p n ) where pi = I for an ISR photon and 
pF = F for a FSR photon, see the analogous construction in refs. |22,p3|. The set of all 
partitions is explicitly the following 

{p} = {(I, 1,1,... ,1), (F,I,I,... ,1), (I,F,I,... ,1), (F,F,I,... ,/),... {F,F,F,...,F)}. 

The s-channel four-momentum in the (possibly) resonant s-channel propagator is X p = 

Pa + Pb- J2 k i- 

pi=i 

The soft (eikonal) amplitude factors s& , uj = I, F, are complex numbers and they are 
defined as follows 



b a (k,p a ) | b a (k u p b ) 



^KPa 



{F} fF\n \ ba\ki,p c ) „ 

5 h = 5 U Kh) = +eQ f ^j— - eQ f 



2kip b 

bg(ki,p d ) 
2kip d 



2 e 2 Q 2 e ( Pa 



Pb 



2 V kiPa kiPf) 



e Qf ( Pc Pd 



2 \hp c kip d 



u- a (k)u a (() V 2(k 



(39) 

2 

) 

(40) 
(41) 



see also Appendix A for more details. As indicated above, the moduli squared of the 
CEEX soft factors coincide up to a normalization constant with the corresponding EEX 
real photon soft factors S(ki). 



3.1.4 Born 

The simplest IR-finite /3-function is just the Born spin amplitude times a certain 
kinematical factor (see the next subsection) 

The Born spin amplitude 03 ( p x ; X) is a basic building block in the construction of all of 
our spin amplitudes - let us define it already at this point. The many equivalent notations 
for 23 will be introduced for flexibility - in view of its role as a basic building block in the 
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calculation of the multi-bremsstrahlung amplitudes. Using Feynman rules and our basic 
massive spinors with definite GPS helicities of Appendix A, Born spin amplitudes foiQ 
the e-(p a )e + (p b ) -> f(p c )f(Pd) process are given by 

<B X) = <B (lllllZ; X) = <B R J] g£] (X) = B MM (X) = 

X] n iTW ( G ^)[H ( G />)[«fl ^ = $Z 



ze 2 



B=j,Z B=j,Z 

( G f M )[H = «(P6, A 6 )G^(p a , A a ), = «(pc, X c )Gf fl v(p d , X d ) 

G lv = 7m J] G L = 7/i w ^9x J , = ^(1 + A75), 



(43) 



11% (X) 



2 

A=± A=± 



X 2 - M B 2 + iT B X 2 /M B 



where gf are the usual chiral (A = +1, — 1 = R, L) coupling constants of the vector boson 
B = 7, Z to fermion / in units of the elementary charge e. If not specified otherwise, the 
"hook function" Hb is trivial H 1 = Hz = 1. It will be used to introduce special effects 
into Born spin amplitudes, like running coupling constants or an additional form-factor 
due to a narrow resonance. 

Spinor products are reorganized with the help of the Chisholm identity, see eq ( |219[) 
in the Appendix A, which applies assuming that electron spinors are massless, and the 
inner product of eq. fl216| ), also in the Appendix A: 

r [ B,e B,f rp rpi . B,e B,f tji tt 1 

fft£ (Y\ ^„2 0x *- x b[9\ a g-Ag i A e A a J-\ h \ d + g Xa 9x a A e Aj, U K\ d \ , ... 

® MM (*) = 2 * e X*-M B * + iT B XyM B > W 



where 



T XcXa =u(p c , X c )u(p a , A ) =S , (p c ,m c ,A c ,p o ,0,A a ), 

=«(p&, X b )v(p d , X d ) = S(p b ,0,-X b ,p d ,-m d ,-X d ), 
u 'x c x b =u(Pc,Xc)v(p b ,-X b ) = S(p c ,m c , X c ,p b ,0, X b ), 
U\ a \ d =u(p a ,-X a )v(p d , X d ) = S(p a ,0,-X a ,p d , -m d ,-X d ). 



(45) 



Note that the use of the Chisholm identity is a technical detail which should not obscure 
the generality of our approach. What we need in practice is any numerical method of 
evaluation of the Born spin amplitudes defined in eq. fl43|) , and Chisholm identity is just 
one possibility. 



3.1.5 Off-space extrapolation 

In eq. (^) Born spin amplitudes are obviously used for pi which do not necessarily obey 
the four-momentum conservation p a +p b = p c + Pd- in the exclusive exponentiation this is 

3 For the moment we require / ^ e. 
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natural and necessary because, in the presence of the bremsstrahlung photons, the relation 
X = Pa+Pb = Pc+Pd may not hold. In eq. (^) only fermion momenta enter as an argument 
of the Born spin amplitudes. Photon momenta play only an indirect role, they disturb 
fermion momenta through energy and momentum conservation (sometimes referred to as 
a "recoil effect"). The self-suggesting questions are: Is this acceptable? Is this dangerous? 
Can this be avoided? The clear answer is: It is unavoidable and natural feature of the 
exclusive exponentiation that certain scattering matrix elements originally defined within 
ra-body phase-space are in fact used in the phase-space with more particles. Let us call it 
off-space extrapolation, analogously to off-shell extrapolation^. It surely makes sense, and 
in principle is not dangerous, provided it is done with a little bit of care. 

A technical remark: In the actual calculations of the multiphoton spin amplitudes 
fermion momenta pi in eq. (fHf ) may be replaced, and occasionally will be replaced, by the 
momentum k of one of the photons. This will be due to purely technical reasons (specific 
to the method of calculating multiphoton spin amplitudes). In such a case, the spinor 
into which k enters as an argument is always understood to be massless. 



3.1.6 Pseudo-flux factor 

One demonstration of the "off-space extrapolation" is the presence of the auxiliary factor 
F = X p / (Pc+Pd) 2 - 111 the framework of CEEX, its presence is not really mandatory and it 
disappears in the "in-space" situation p a + Pb = Pc + Pd- In other words, the F-factor does 
not affect the soft limit; it really matters if at least one hard FSR photon is present. It 
is not related to narrow resonances, but rather to the leading-logarithmic (LL) structure 
of the higher-orders. Nevertheless, the F-factor is useful, because it is already implicitly 
present in the photon emission matrix element at 0(a l ) and in all higher-orders, as can 
be seen in the LL approximation. It is therefore natural to include it at the early stage, 
already in the O(a ) exponentiation. If we do not include it at the O(a ) then it will be 
included order by order anyway. However, in such a case, the convergence of perturbative 
expansion will be deteriorated. As we shall see below, the introduction of the F-factor will 
slightly complicate the higher-order exponentiation and construction of the (3 functions, 
but the gain is worth the effort. Furthermore, the F-factor has also been always present in 
the "crude distribution" in the YFS-type Monte Carlo generators, see for instance ref . [|7| , 
so it also improves the variance of the MC weight, especially for C(a°)cEEX- 



3.2 IR structure in CEEX 

Let us discuss in detail the origin of the C(a r )cEEX expressions eqs. d3T|- 32j) and the 
mechanism of the IR cancellations. Our real starting point is the infinite order pertur- 
bative expression for the total cross section given by the standard quantum-mechanical 

4 In the off-shell case particles do not obey p 2 = m 2 , here we also modify the dimension of the phase- 
space. 
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expression of the type "matrix element squared modulus times phase-space" 
00 1 f 1 

ff(Bfl) = E^/^» + w^^'-'^4 E I^KSS-'-fc)! 2 . ( 46 ) 

where dr n is the respective 727 + 2/ Lorentz invariant phase space, and M n are the cor- 
responding spin amplitudes. To simplify the discussion we take the unpolarized case, 
without narrow resonances. 



3.2.1 IR virtual factorization to infinite-order 

According to the Yennie-Frautschi-Suura fundamental factorization theorem |J , all virtual 
IR corrections can be re-located into an exponential form-factor| order by order and in 
infinite order 

Jy[(oo) _ e aB 4 {p a ,p b ,p c ,p d )^{oo) ^ 

As the convergence of the perturbative series is questionable, the above equation is in 
practice treated as a symbolic representation of the order-by-order relation which at 0(a r ) 
reads 

r—n , „ \ r —l 

M^El^^ (48) 



1=0 



where the index I is the number of loops in 2Jtn +n '. The above identity is quite powerful 
because 9Jtn +n ' are not only free of the virtual IR-divergences, but are also universal: they 
are the same in every perturbative order r - for example for one photon, the one-loop (IR- 
subtracted) component, SDT^ , is the same in the fifth-order and, let us say, in the second 
order, where it appears for the first time. The above identity can also be reformulated as 
follows 

r—n 

= y^ml +n] = \e- aB ^>P»^MP] , (49) 

' O(ot r ) 

1=0 y ' 

where, JYti^ has to be calculated from Feynman diagrams in at least[] 0(a r ). The above 
steps are exactly the same as in ||. 

The YFS form-factor £? 4 for e~{p a ) + e + {p b ) — > f(p c ) + f{Vd) + n l reads 



aB 4 (p a ,p b ,p c ,p d ) = J k2 _ d ^ 2+K j^y \Ji{k) - J F (k)\ 2 , 

2p f } + & 



(50) 



Ji = eQ e (J a (k) - J b (k)), J F = eQ f (J c (k) - J d (k)), Jf(k) 



k 2 + 2k -p f + ie 



5 In the LL approximation it is, of course, the doubly-logarithmic Sudakov form-factor. 

6 The use of %ft n r+m ' at 0{cS r+m ^), m > will yield the same result - this is another way of stating 



the universality property. 
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Using the identity (^2 k Z k J k ) 2 = - ^2 i>k ZiZ k (Ji - J fc ) 2 , valid for J2 z k = 0, where Z k 
is the charge or minus charge of the particle in the initial- or final-state respectively, we 
may cast (see ref. ||) B4 into a sum of the simpler dipole components 

B^PaiPbiPaPd) = Q 2 e B 2 (p a ,Pb) + Q 2 fB 2 (p c ,Pb) /r1 . 

(51) 

+ Q e QfB 2 (p a ,Pc) + QeQfB 2 (p b ,Pd) ~ Q eQ f B 2 (p a , p d ) - Q e Q fB 2 (p b , p c ) , 

B 2 ( Pi , Pj ) = J k2 _ d J 2+u j^f (j(Pi, k) - J( Pj , k)) 2 . (52) 

In the above we assume that IR singularities are regularized with a finite photon mass m 7 
which enters into all B 2 s and implicitly into s-factors (and in the real photon phase-space 
integrals, see the following discussion). 

3.2.2 IR real factorization to infinite-order 

The next step is isolation of the real IR singularities and it is worth to elaborate on 
this point because here the CEEX method differs in essential details from the original 
YFS method |J. We use again results of the basic analysis of real IR singularities of 
ref. ||, the essential difference is that we do not square the amplitudes immediately - it 
is done numerically at the later stage. The validity of the whole basic analysis of the IR 
cancellations in ref. || remains, however, useful because it is done in terms of the currents 

jf(k) = -^^,f = a,b,c,d. (53) 

The above currents are simply related to our s-factors: 

s { J } (k) = const x Q e (j a - j b ) ■ e CT (/3), 
si F} (k) = const x Q f (j c - j d ) ■ e a {(3). 

It is important to remember that the whole structure of the real IR divergences is entirely 
controlled by the squares of the currents \j(k) | 2 , for j = j a —jb or j = j c —jd, independently 
whether we prefer to work with the amplitudes or their squares, because only the squares 
|j(/c)| 2 are IR divergent and the other contractions do not matter (as was already stressed 
in ref. ||). Similarly, if we express spin amplitudes in terms of s-factors, only the squares 
|s(/;;)| 2 are IR divergent and not the interference terms like dt{s(k)(. . . )*}. 

Having the above in mind we may proceed using results of ref. || and we see that for 
instance the most IR divergent part of M n is proportional to the products of n s-factors 

m n . . .£) ~ A (5; x) Sai (h) Sa2 (k 2 ) . .. 5an (k n ) (55) 

where the function (3q is not IR divergent any more, and we assumed for the moment the 
absence of the narrow resonances, using the sum of ISR and FSR s-factors[] 

sjjf) =sf>(fc)+s«(fc). (56) 

7 In the non-resonant case we may set X = p a + Pb, for example. 
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However, there are also non-leading IR singularities. Suppressing inessential spin 
indices the whole real IR structure is revealed in the following decomposition || : 

n n 

Wl^ih, k 2 , h, k n ) = A, J[*(fe a ) + U<k s ) 

s=l j=l s^j 

+ J2$2(k jl ,k j2 ) Y[ s(k s ) + M k h,k h ,k j3 ) \\ s(k s ) + ... ^ 

h>h s^ji,j 2 ji>32>h s^ji ,32,33 

n 

+ y^^ ra _i(fei, ...kj-i, k j+1 , fc n )s(%) + Pn{h, h, h, k n ) 
3=1 

where functions /3j are IR free and include finite loop corrections to infinite-order. Let us 
stress that these functions $ are genuinely new objects. They were not used and even not 
considered in ref. ||. 

3.2.3 Finite-order /?'s 

The decomposition of eq. (|57|) has also its order-by-order representation, which at 0(a r ), 
r = n + /, reads as follows: 

n n 

wtt? +i \h, h, k 3 , k n ) = &® i[ 5 (k s ) +J2ft 1+l) ( k j) n*) 

8=1 i =1 

jl<j2 s^n,32 h<32<33 s^j 1 ,j 2 ,33 

n 

+ E - Vl» -> *n>(*j) + P { n +l \h, k 2 , k 3 , k n ) 

= TT sffe 1 \ ft® 1 @i + 1 ^ + , + ^(fejn fcj2) 

11 l S yo + 0(%) S (% 2 ) h ^ <j3 B{k h )B{k h <k h ) 

>p -fcj-i, fcj+i, fen) pt +l \k 1 ,k 2 ,h,...,k n ) 

U n*(k s ) + n*(k s ) 

(58) 

The new functions fin (hi, k 2 , k$, fe„) contain up to /-loop corrections, and are not 
only completely IR-finite, but are also universal: for instance the (5\ (k), which appears 
for the first time in decomposition of Wt[ (k), is functionally the same when decomposing 
Wl 2 3 \ki, k 2 ) or any higher-order . This feature is essential for reversing the relations 

of eq. fl5"8"D, that is for practical order- by-order calculations of (3^ +l ^ from DJln \ obtained 
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directly from the Feynman rules: 

fi? +l \h, k 2 ) = Wlf +l \h, k 2 ) - p[ 1+l \h) 5 (k 2 ) - p[ l+l \k 2 )s(h) - ffisfaXfo), 
fii 3+l \h,k 2 ,k 3 ) = wi 3+l) (h,k 2 ,k 3 ) 

- fi? +l \h, k 2 )s(k 3 ) - $ +l \h, k 3 )s(k 2 ) - P { 2 +l \k 2) k 3 )s(h) 

- ^(hXhXh) - fi[ 1+l \k 2 X k iXh) - fi^ihXhXh) 

n 

fin + ^(^1) • ) k n) = y^ri + ■'•■> k n) ~ ^ ] fin— 1 (^1) •••%—!> •••> k n )s(kj) 



ii<ja s&uh i =1 s=1 

The above set of equations is a recursive rule, i.e., higher-order /3's are constructed in 
terms of lower-order ones. In practical calculations we do not go to infinite-order but 
we stop at some and the above set of equations is truncated for {3 n n+l ' ) by the 

requirement n + I < r. The above truncation is harmless from the point of view of IR 
cancellations because we omit higher-order fi's which are IR-finite. As a consequence of 
the above fixed-order truncation eq. (|57D takes the following form: 

Tt^(h,k 2 ,k 3: ... : k n ) = 

fi[\kj) fi 2 \f £ j l ,kj 2 ) ^ fii i k ju k jn kj 3 ) 



n*) K , +E^+E3z%^+ E 



j=i S ^ ~^ j2 s ( k hM k n) h ^ <h 3 { k h) s ( k h s ( k h) (60) 

+ E 



fir ( k jl> k hl ■■■■> k 3r) 



... ; h < k hH k h)~*{ k i 

where, contrary to eq. (|58|), we now allow only for r < n; in such a case the sum has r + 1 
terms instead of n. 

The above formula represents the general finite-order 0(a r ) exp case while for r = 
only the first term survives, and in our 0(a 2 ) case there are three terms. The CEEX spin 
amplitudes in our master formula eq. ( |3"3"| ) represent the cases of r = 0, 1, 2. 

Just to give an explicit example, in the recursive calculation of /3's in 0(a 3 ) we would 
need to calculate j3$\l = 0,1,2,3, = 0,1,2, fi { 2 +l \l = 0,1 and /^ 3) . In 

the present work, at 0(a r ), r = 0,1,2, we shall employ the following set of recursive 
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definitions based on eqs. ([59] 



A? } (a) = awS (5), 1 = 0,1,2, 

$i+0 g S) = ^(i+O gfa) _^(0 (P),^), / = 0,1, 

(61) 



where the 9Jt-amplitude is given by eq. (£49]). Here we restored spin indices but we still 
specialize to the non-resonant case, and our /3's do not have the partition dependent X p 
argument as in /3's of eqs. ( |33| - |35|) . We shall provide a definition for /3's in the resonant 
case in the following section |3.3.4| . 

3.2.4 IR cancellations in CEEX 

At fixed-order 0(cn r )cEEX, and remembering that | exp(_B 4 )| 2 = exp(29?_B 4 ), we have ob- 
tained 



a 



(r) 



oo - 

Y,-\j dr ^ + **>■■■> e™^--^- £ iJOTf) (h, k 2 , . . . fc 

n=0 ' spin 



2 

(62) 



where DJln^ is defined in eq. fl60"D and we factorize out the s-factors 

1 2 ™ 

- ^ |3rtjp(fci, A: 2 , fc 3 , K) | = d n (k 1 ,k 2 ,k 3 , ...,k n ) ]\\s(k s 



)\ 5 



spin 



d n (ki, k 2 , k%, fc re ) 



A(r-) A (%) $2 ^ii' Aj 2 ) | ^-v /3 2 fcj 2 , 



J3 ' 



/3r ^(A;^, kj 2 , /c Jr ) 



; s(%)s(A; i2 )...s(%J 



(63) 



In the above the function d n (ki, k 2 , k^, k n ) is IR-finite and we are allowed set m 7 — > 
in it. Apart from 2a9?£> 4 the IR regulator m 7 still remains in all s{ki) -factors and in the 
lower phase-space boundary of all real photons in J d 3 k/2k°. 

The above total cross section is perfectly IR-finite, as can be checked with a little bit 
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of effort by analytical partial differentiation^ with respect the photon mass 
d 



dm 1 



a {r) 



n=0 ^ spin 

1 n f 

+ E ~1 E / rfr «-l( P ; ^3,P4, fcl, • • • , fc,-l, fc s +l, • • • , fcn) e 2 ° RB4 
„=1 n - s =l ^ 

/ ^ |s(A:s)|2 } n |s( ^ )|2 fcs > ^ 



,2 



9 r /■ d j jfe, 

X - — 



(64) 



It is now necessary to notice that 

d f f d 3 k. 



<9m 7 1 ./ 2k2 



is a 5-like function concentrated at k s = and we may therefore use the limit 

dn{ki, • k s , /c n ) > d n (ki, /c2, k s — i, 0, = d n —i[ki, k^-, k s —±, k s -\-i, fcn) 

The above helps us to notice that all terms in X^=i are identical and we may sum them 
up, (after formally renaming the photon integration variables in the second integral) and 
rewrite eq. (|64]) as follows 



oo 

Y,-\J P3,P4, fcl, • • • , fen) e 2Q ^ 4 - J] W (h, h,... k n ) 

n=0 spin 



(65) 



where the independence on m 7 of the sum of the 1-photon real and virtual integrals is 
due to the usual cancellation of the IR- divergences in the YFS scheme, shown explicitly 
many times. 

The integral of eqs. ( [4"B] ) and ( p2] ) are perfectly implementable in the Monte Carlo 
form, with small m 7 being the IR regulator, using a method very similar to that in ref. J7| . 
Traditionally, however, the lower boundary on the real soft photons is defined using the 
energy cut condition k° > ey/s/2 in the laboratory frame. The practical advantage of 
such a cut is the lower photon multiplicity in the MC simulation, and consequently a 
faster computer program^. If the above energy cut on the photon energy is adopted, then 

8 This method of validating IR-finiteness was noticed by G. Burgers |2£j. The classical method of 
ref. H relies on the techniques of the Melin transform, which could be also used here. 

9 The disadvantage of the cut k° > Ey/s/2 is that in the MC it has to be implemented in different 
reference frames for ISR and for FSR - this costs the additional delicate procedure of bringing these two 
boundaries together, see ref. El and/or discussion in the analogous i-channcl case in ref. [jflj. 
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the real soft-photon integral between the lower LIPS boundary defined by m 7 and that 
defined by e can be evaluated by hand and summed up rigorously (the only approximation 
is m 7 /m e — > 0) in the following. 



3.2.5 Explicit IR boundary for real photons 

A general notation for the IR domain Q was already introduced, see eq. (|36|). Let us 
now exclude the Q domain from the real photon phase space (integrate out analytically). 
Splitting the real photon integration phase space we rewrite the eq. (^2[) as follows 

n=0 ' j=l J 3 J 

( 66 ) 

/ dT 2 (P-J2 k 3-, Ps,p 4 ) e 2a ^d n (h,k 2 ,...,k n ). 

After expanding the binomial product into 2 n terms let us consider for instance the sum 
of all (™) = n terms in which one photon is in Q and the other ones are not: 



s=l s j^s 3 

/n 
dr 2 {P - y^fcj-; p 3 ,Pi) e 2a ^ B4 d n (k 1 ,k 2 ,...,k s _ 1 ,0,k s+ i,...,k n 

7=1 



(67) 

k 's(k)\ 2 e(n,k) 



in 



WlJ 2k° 



/n—l 
dT n+1 (P; p 3 ,P4, h, h, ■-, K-i) Yl ®(^' k j)\ B ( k j)\ 2 d n -i{ki, h, • -, K-i) 

3=1 

A similar summation is performed for the ( n ) terms where s photons are in Q giving rise 
to 



a 



oo ^ n 

(r) 



£^£C)(/§wW)y 

n=0 s=0 w 7 

/n—s 
dr 2+n _ s (P; p 3 ,p A , k u k 2l k s ) JJ { Xfcj) j 2 0(Q, kj)} e 2a?liBi d n _ s (k u k 2 , k n _ s ) 

J dr 2+n (P; p 3 ,p A ,h,k 2 ,...,k n )ex V ^\ 5 (k 3 )\ 2 Q{^k 3 )^ e ™( Pl ,.., P4 ) 

X II i l S (^)| 2 @( fi ' k i)} d n(ku k 2, k n ). 
3=1 

(68) 
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The additional overall exponential factor contains the well known function 



/d 3 k- 
-^■Hkj^Q^kj) = Q 2 e B 2 ( Pl ,p 2 ) + Q 2 f B 2 (p 3 ,p A ) 

+ QeQfB 2 (pi,P3) + QeQfB 2 (p 2 ,p 4 ) - QeQfB 2 (pi,P4) ~ QeQ fB 2 (p 2 , p 3 ) , 

b iM s -f*L m , kj) (*(*)-*<*))' , / § m , kj) <g> (|L _ £) 2 , 

(69) 

which forms together with 2o$tBjj)i, ...,P4) the conventional YFS form-factor 

Y{n-p u ...,p 4 ) = 2aB A ( Pl , ...,p 4 ) + 2aKB 4 (pi, ...,p 4 ) (70) 



in our master eqs. ( pl| , |32|) . The dependence on m 7 in Y cancels out. Photon mass gets 
effectively replaced by the size of Q in its role of the IR regulator. The YFS form-factor 
Y can be decomposed into six dipole components, see eq. (j37D and can be calculated 
analytically in terms of logs and Spence functions, see refs. [0O| [32|1 keeping all fermion 
masses exactly. 

As already indicated, in the MC with the YFS exponentiation it would be possible 
to do without Q (declare it as empty) and rely uniquely on the IR regularization with 
a small photon mass m 7 only ||. In such a case the formulas (|38D for YFS form factor 
would include only the second virtual photon integral part. 

For the sake of the completeness of the discussion it is necessary to examine once again 
the IR cancellations in the total cross section with Q as the new IR-regulator 

oo .. „ n 

a(r) = E- / dT ^(P; p 3 ,P4,k 1 ,k 2 ,...,k n )l[{\ 5 (k J )\ 2 e(n,k J )} 

n=0 ^ J j=l (71) 

x e B4(n; P1 ,..., P4 )+2a5RB4(pi,..,P4) ^ ^ fc ^_ 

IR-finiteness of the total cross section now simply translates into independence on the Q 
domain, (assuming, as usual, that the size of Q is very small) 

The proof can be done along the same lines as the previous one for the photon mass. Let 
us assume that we want to vary Q — > Q' = Q + 5Q, that is Cl' = Cl — 5Q. Note that Q' can 
be much bigger or smaller than of Q, the only requirement is that both are very smallj^] 

10 <5O does not need to be infinitesimal with respect to f2. 
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We proceed as follows 

n=0 ' j=l <J j J 3 > 

x J dr 2 (P - kj\ V^Va) e S4(n ; pi 1 ...,P4)+2a^ 4 (p ll ... )P4 ) d ^ ^ ? ^ 

= EAEQ { / Sl*(*)l a e(«U)} / dr 2+n . s (p ; ps^ku-M 

x n{K%)| 2 ©(^%)} e ^( n ^--^ +2aRB *^'-**) d n - s (ki,k 2 , ...,k n - s ) 

3=1 

00 -IP 3 

X II i |S(^)| 2 ©(^, fcj)} dn(fcl, fc 2 , fcn)- 

(73) 

recovering the same expression as ([711), but with f2' instead of Q. 
3.3 Narrow neutral resonance in CEEX 

The main new feature of CEEX in comparison with EEX is that the separation of the 
IR real singularities is done at the spin amplitude level and after squaring and spin 
summing them (numerically) the higher order terms are retained while in CEEX they 
are truncated. For more detailed discussion of see section |4.3| , where we show explicitly 
the relations among (3's of EEX and /3's of EEX. Keeping the above in mind, we still 
have at least three possible versions of CEEX. In the following we shall describe them, 
concentrating mostly on the third one which is designed for the neutral resonances^ 
and which is the principal version implemented in the /C/C Monte Carlo. Let us stress 
immediately that the resonance may be arbitrarily narrow. However, our approach works 
without any modification for any value of the resonance width. 

3.3.1 General discussion 

We believe that CEEX is the only workable technique for treatment of narrow resonances 
in the exclusive MC. To understand the essential difference among three possible formu- 
lations of CEEX it is enough to limit the discussion to the simplest case of the O(a ). 
The three possible options are: 

11 The case of exponentiation for charged resonances like W ± resonances is not yet covered in the 
literature. 
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(A) Version for the non-resonant Born without partitions: 

n 

^ ••■£) = II « (**) + < (**)) ® HH (74) 

i=l 

(B) Version for the non-resonant Born with partitions: 

n 

anS°g£5 •••*:) = E II »hhW (75) 

pe{/,F}« i=l 

(C) Version for the resonant Born: 

^ g££ • • •» = E f[ 4f } (^) E «hW 

(76) 

Let us immediately define the additional form-factor for the Z resonance (case (C)) 

°mv - J w ^ T -_ v*) w (*))• - ') • (") 

where M 2 = M| — iMz^z, the currents J M are defined in fl50|) , while for the non- resonant 
part we have Bj(X) = 0. The Bf(X) form-factor sums up to infinite order the virtual 
aln(r z I 'Mz) contributions - we postpone discussion of its origin and importance to the 
latter part of this section. 

Coming back to the more elementary level we see that the case (B) becomes (A) if we 
can neglect the partition dependence of the four momentum in the Born: Q5[6 a ][ C£ q(X p ) — > 
%$[ba][cd\(P)i where P = p a + pi, or P = p c + p d or any other choice which does not depend 
on momenta of the individual photons. This is thanks to the identity: 

n n 
1=1 P6{/,F} i=l 

Only case (C) is efficient for the resonant process, so obviously (A) and (B) are limited 
to non-resonant processes. The immediate question is: which of them is better? If (A) 
is not summing higher order much better than (B), then it has the clear advantage of 
being simpler - summation over partitions makes the computer code more complicated 
and adds heavily to the consumption of CPU time. The answer is that, although we 
did not investigate quantitatively the differences between (A) and (B), we think that (B) 
sums up the LL higher orders more efficiently than (A) and is therefore better, even if 
there is no resonance. In our case, since we want to cover the resonant process anyway, 
it is a natural choice to use (B) for the non-resonant background component of the spin 
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amplitudes (off-shell 7 exchange) even if it is not vital. Once summation over partitions is 
in place, it is the easiest to use it for the non-resonant background as well. The additional 
bonus of better higher order convergence provides an extra justification. Summarizing, if 
(C) is implemented then (B) comes for free. 

Having discussed the differences among the three options let us now concentrate on 
option (C) for the resonant process, remembering that for the non-resonant background 
component it becomes automatically (B). First of all, for the narrow neutral resonance (the 
Z boson in our case) the emission of the photons in the production and the decay processes 
are well separated by a long time interval, and are therefore completely independent and 
uncorrelated. In the perturbative QED this simple physical fact is reflected in a certain 
specific class of cancellations among the ISR and FSR photons on one hand and among 
virtual and real corrections on the other hand. For the inclusive observables like the total 
cross section or charge asymmetry the effects of the ISR-FSR interference (IFI) in the non- 
resonant case are of order a/ir, typically up to 1%, as can be seen from many example of 
the explicit 0{a 1 ) calculations. The IFI effect will be of order (a/ii)(E mSlX /Ebeam), when 
the experimental cut on photon energy is E max . Note that the IFI effect is not directly 
enhanced by the big mass-logarithms like Ins jm\ ~ 20. For the resonant process the 
IFI effects in the inclusive observables are of order (a/ir)(Y /M) and are therefore often 
negligible on the scale of the experimental error. One has to remember, however, that the 
additional suppression factor Y/M disappears if the experimental cut on photon energy 
is of order of the resonance width, E mauX /Eb eam ~ Y/M, and for an even stronger cut 
-Emax < T the IFI effect becomes of order (a/7i)(E max /Y). 

If Y/M is extremely small, like for the r lepton, the IFI cancellation can be taken 
for granted and the photon emission interference between production and decay can be 
neglected whatsoever. In the case of the Z boson close to the Z resonance (LEP1) the 
IFI effect is detectable experimentally but it is small enough that it can be omitted in the 
Monte Carlo programs used for correcting for the detector acceptance only. In this case 
KORALZ/YFS3 |10| with the EEX matrix element was the acceptable solution. 

The most convenient solution is the universal Monte Carlo in which IFI is included, 
which can evaluate IFI effects near the resonance, far from the resonance, for inclusive 
quantities and for strong energy cuts E max ~ Y. This is exactly what our CEEX offers. 



3.3.2 Derivation of the resonance formfactor 



As we have already pointed out (following refs. |22|,|23|]), in the presence of narrow res 



onances it is not enough to sum up coherently the real emissions, taking properly into 
account energy shift in the resonance propagator (only due to ISR photons). It is also 
necessary to do the same for the virtual emission, and also sum them up to infinite-order 
- this is why the resonance form factor exp(i?^) has to be included, see eqs. (|76|) and 
(|77|) . In the following we shall derive eq. fl77|) for Bf and show analytically that the IFI 
cancellations do really work, as expected, to infinite order. 

Let us write again the formula for standard YFS function in eq. (15(1) in a slightly 
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modified notation 

S(k) = S I (k) + S F (k) + S Int (k), (79) 
S I (k) = \J I (k)\ 2 , S F (k) = \J F (k)\ 2 , S Int (k) = -23?(J/(fc) ■ J* F {k)) 

In the presence of the narrow resonance, the YFS factorization of the virtual IR con- 
tributions has to take into account the dependence of the scalar part of the resonance 
propagator on photon energies of order V (the numerator treated in soft photon approxi- 
mation as usual). The relevant integrals with n virtual photons look as follows: 

^_^)^ E n/^^ SpA) * (80) 

n=0 p£V n i=0 J V ' 1 T P 

where M 2 = M 2 —iMT, and V n is set of all 3 n partitions (p±, p 2 , p n ) with pi = I, F, Int, 
and P p = P — ^2 ki includes only momenta of photons in Si nt and not of photons in Si 

Pi=Int 

or Sp. The (P 2 — M 2 ) factor is conventional, to make the integral dimensionless. We shall 
show that the above integral factorizes into the conventional YFS formfactor (dependent 
on the photon mass m 7 ) and the additional non-IR factor due to the resonance R = Z 

I = exp(5f (m 7 , s, M)) = exp( J B 4 (m 7 , s) + A#f (s, M)). (81) 

Our aim is to find the analytical form of the the additional function ABf. In the current 



calculation we use the following approximate formula, also used by Greco et.al. f22|,|23 

aAflf M = -2Q eQf a - in g) .„ (^) = In . (82) 

In the following: 

• We shall derive the above approximate result and 

• show explicitly that the above approximate virtual interference part of the formfac- 
tor cancels exactly with the corresponding real interference contributions. 



Since soft virtual photons entering into Sj and Sp in eq. (|80| ) do not enter the resonance 
propagator, we may therefore factorize and sum up the contributions with Si and Sp: 

°° 1 ni f i d^k- °° 1 ™ 2 f i d^k- 

1 = e — { n j {2 l)3 k 2 - ii m 2 s ^ e — x n j ( 2^)3 k 2 - 2 m 2 SF ^ 

ni =0 1 i 1= J \ J H 7 n 2 =0 z i 2 =0 K ' 12 T 



x J? ^ S/ ( 2 -) 3 ^-%^ ( " t3) (p - e;:!^) 2 - m 2 (83) 

00 _ n „ ,4 



^ n! I! / (2 I)3 *W*0 (P _ E / fc 



j; 2 — M 2 
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Now we neglect the quadratic terms in photon energies O(kikj) 

1 111 



\n 2Pfc, 



(P - YJ 1 . kA 2 — M 2 P 2 - 2P V n , kj — M 2 P 2 - M 2 i - y n _ 
J ^ 1 1 " P 2 — M 2 

_ M 2 11 i 2P~ - p 2 _ A>2 11 



J4) 



P 2 



P 2 -M 2 



P 2 - M 2 AA ( P _ fc,) 2 - 



M 2 



and this leads to 



J = e 



APf(r) 



aBj+aBp 



exp 



(2tt) 3 /c 2 - m 2 



ze 



'Sji n (k) 



P 2 -M 2 



(P - A;) 2 - M 2 



3 S 4 (m 7 )+ABf (r) 



(2^)3 1 2 " 



Slin(k) 



P 2 -M 2 



(P - A;) 2 - M 2 



- 1 



55) 



How solid is the above "derivation"? Strictly speaking it is justified in the limit where we 
follow Yennie, Frautschi and Suura in ref. || and express the k — > emission amplitude 

as 

M—>t(ei + 0(k/M) + ^-(e 2 + 0{k/M)) 
k\ T z 

where e\ t i are constants independent of k, so that 

\lPkjj{P 2 - A? 2 ) | < 1, 

that is if photon energy is below the resonance width. This restriction is thus entirely 
analogous to the usual YFS expansion into IR-singular part and the rest. We note that 
Greco et.al. in refs. [p2,|23| have also pointed out that the result for AP^(r) in eq. ( |35| ) 



follows from the YFS expansion; here shall show how this happens in detail. 

The best situation would be to have a more precise evaluation of the integral of eq. fl8C?D 
(the integral is probably calculable analytically). For the moment, however, following 
refs. |22|,[23|]: we choose an easier "pragmatic" approach based on the fact that the virtual 
and real contributions from IFI for photons E 1 > T do cancel, as a consequence of the time 
separation between production and decay, and we shall check that the above cancellation 
really works. In this way we trade analytical evaluation of the more difficult multiphoton 
virtual integral for an easier evaluation of the multiphoton real integral. 



3.3.3 Cancellation of the virtual formfactor with the real emissions 

Let us therefore examine analytically the real multi-photon emission contribution from 
the IFI0- The starting point is the integral in which the total photon energy K = X]j=i kj 

12 In the practical CEEX calculation the contribution from IFI is evaluated numerically, inside the MC 
program. 
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is kept below E max = iwv^, where T < E max « y/s: 

-Eyfl|E E n«wta) 3^ W„ - E *,) 



n=0 u i=l 1 <n...(r n p e{I,F}"j=l P 

°° 1 f n d 3 k- n 

/ 11^ Yl J2 H s IMpj} s \j]{ P r} 

n=0 Vol /r i=1 1 ci-ovi p,p'e{/,F}" j=i 



.V 2 - M 2 I X 2 , - M 2 



ftr Q < „ A /3 — »l-»np,p'e{I,F}nj=l P \ P 



oo 



e^t / n^sr e n mm n ^ 

e aBf{P^K I -K IF ) / ^BfiP-Kt-Kpj) 



Slnt(kj) S Int {kj 



(P - Kl - Kjp) 2 -M 2 \(P- A, - Kpi) 2 - M 2 ^ 

(86) 

where the Born amplitude we have simplified to the level of the scalar part of the resonance 
propagator and we denote 

&(*i) = £i«wi a > ^*i) = Ei*S } i 2 > 

Sint{kj) = E s fj] } (%]«)* = ( 5 b']{°})*' 

(87) 



A /F = £ = E fc J' K = K i 2 + K F 2 + K IF + K FI 

Pj =IF Pj =FI 

As we see, the product of two sums, each over 2 n partitions p, p' G {/, F} n , is now replaced 
by the single sum over 4 n partitions p, G {I 2 , F 2 , IF, FI} n , where IF,FI represent the 
interference terms. 

Keeping track of the dependence of the propagators on Kp, K IF and K FI , the sum- 
mation over the number of photons can be reorganised, leading us back to the following 
factorized formula 

ni J 3 b oo n 2 ,3, 

<*w) = > J -t / TT ^Sr 2S,(k„) V-r/ TTTrnr 2S F (h„ 



E^/n^^..)E^/n^^ 

m=0 1 J ii=l *i n 2 =0 ^ ^ j 2 =l * 2 

1 f H 3 k- ~ paBfiP-Kj.. 

~ 1 /• J^t f/ 3 fr / p aB?(P-K l2 -K FI ) 

E^/n^",^ 



2 — M 2 



(88) 



ri4=0 ^4 = 1 l 4 

6(P max - X /2 - K F 2 - Kj F - Kpj 
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(P - A> - K FI ) 2 - M 2 



where Kp = ^ k h , K F % = Y, i2 h 2 , K IF = Y^ ia h 3 and K FI = Yl u h 4 - The sums over 
the pure initial and final state contributions, and the interference contributions are now 
well factorized and can be performed analytically. As a first step we integrate and sum 
up contributions from the very soft photons below e^/s, similarly to what was shown in 
ref. 0: 

I* -Em ax /* ax 

d(tj max ) = / dE S(E — Ei — Ep — E Int ) / dEi dE F dE IF dE FI 



o 



e A ft / ^ste^w^sfr - y: k) 

E A ft / ^ F (k i2 ) f***W+**B* S{EF _ k 

~ 1 ^± f d 3 ki ~ e aAB?(P-K l2 -K IF ) 

» 1 "1 f (fki ~ ( e a^{P-K l2 -K FI ) 

e 2QSAB ^(£M-^^ 3 ) 5 



«3 



(89) 

where E = The integration over photon momenta can be performed without any 



2 

approximation leading to the following result 



O-Omax) = / dv 8{V-V 1 -Vp-V IF -V FI ) 

JO 

J dvj F^jrv?- 1 ^Bm+^tBi J dv F F^pvf- 1 e ^F(E)+2^B F 

r ,y T t 1 1 i / p oABj( a (l-«j)(l-wi F )) \ . (90) 

J dVlF b [ 2 j 2 T/n ^ U(i - - v IF ) -m) e 

J FI [ 2 } 2 llntVpi \s(l - Vl )(l - v FI ) -m) 6 

in which is explicitly free of any IR divergences. 

The essential question is whether we have perfect cancellations of the ln(T/Mz) terms 
in the interference subintegral 

fV m ax-Vi-v F -v FI /^ Tr ,\ 1 l -, p aABf(s'(l-v IF )) 



l Int 



M / dv JF F(^) L llFV f ^~ l _£ (91) 



We omit from consideration the constant IR-finite factor e aB i™t(E)+a$iB Int Decause ft does 
not depend on resonance parameters. The bulk of the integral comes from the neighbour- 
hood of vj F = and the integrand is ~ 1/v 2 at large v due to the resonance; we can 
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therefore extend the integration limit to / °° dv[ n t at the expense of an error of 
One possible evaluation method is to use the standard techniques of the complex func- 
tions. First, we reformulate the integral as an integral over the discontinuity C\ along the 
real axisQ 

hm = FOf)e^ \ \ I dz - ± (92) 

V 2 / i sm(7r|7/ nt ) J Cl 2 s' — M 2 — s'z 

Since the contour can be closed in a standard way with the big circle, the integral is given 
by the value of the residue at z — 1 — M 2 /s'. 

1 / A Tin* -1 „ 

1 M 2 -s'\ 1 



\ 2 / sm (ix-^ij. 



l lint 



M 2 — s' \ 2 ) sin(7ri 77rii ) 6 { s> 

1 



2 

7 (1 + 0{ llnt )) 



M 2 - s 
The above is true because 

M = - 2Q < !n g) !n (^) = -i 7M >n (^) (94, 

We have therefore proven the full cancellation of the dependence on the resonance pa- 
rameters for the integrated cross section. 

3.3.4 Definitions of /3's with partitions 

The 0(a r ), r = 0, 1, 2, /3-functions for the variant of the CEEX with summation over the 
partitions, as in eqs. (|33| - |35D , are derived with the recursive relations of eqs. fl59|) (similar 
to those of eqs. floT])). The only additional complication is that we must keep track of the 
indices which say whether an external real photon is of ISR or FSR type and of the total 
photon momentum after emission of the ISR photons (the one which enters resonance 
propagator, if such a resonance is present) 

$g>&P)=mg>&P), / = 0,1,2, 

0$$ (x^ p - = Kn ~ ki) - /f g; p - h) i = o, l, 
( p xa\;P) = an£J? g£; p ) - $° g; 4f i = o, 1, 

o(2) (Pkihi. Y \ — m?( 2) (Pkik 2 . v \ 

- /d es; *j - #&, } g& *j - ^ 0) g; 4r } (*4)4: 2} (u 

(95) 

13 We have also pulled out of the integral the e aAB * factor, because the most of integral comes from 
the neighbourhood of the singularity at vjf = 0. 
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where X w = P — Yl P = Pa + Vb- Introduction of the partition index Ui defin- 

U)i=I 

ing whether a photon belongs to ISR or FSR is in a sense not such a deep and great 
complication - it is now just another (third) attribute of the photon like its helicity. 

Let us look closer into the structure of the term like /^r^i (^|; J^js^ 2 ^^)- For 
example io\ = F and U2 — I it reads P[^ F ^ (^; P — ^2) sl?^), that is the total shift in 

X in ft 1 * depends not only on the type uj\ of "its own photon" but also on the type cu 2 of 
the photon in s^ 2 ^ factor which multiplies it! 

The 971-amplitude in eq. (|95|) is given essentially by eq. (|49"D with the formfactor 
including the resonance part (if present) 

STtfWfi (Pki kn. v \ — L-aB4-oBf(X u )'v r (r)i? (phi k n . X \] ( qfl \ 



0(a r ) 



As we see the type R = 7, Z of the "resonance" formfactor Sf 1 has to be adjusted to the 
type of the component in (we have temporarily introduced an explicit index R into 

M and QJt and 7 is essentially a "resonance" with the zero width). 



3.4 Virtual corrections, no real photons 

We now start to accumulate the actual formulas for the /3-functions entering the CEEX 
amplitudes of in eqs. (|33| - p5[ ) with the case of no real photons and up to two virtual 
photons. The "raw material" are the M-amplitudes from Feynman diagrams which are 
turned into /3-functions using the recursive relations of eqs. (pop. 



3.4.1 Photonic corrections 



a c a c a c a c 




Figure 3: First order diagrams. 



Let us start with the simple case of 0{a 1 ) spin amplitudes with one virtual and zero 
real photon coming directly from Feynman diagrams, which will be used to obtain the 
first order (3^\ The relevant spin amplitudes are 

M« (I; X) = 25 (5; X) [1 + QlF^s, m e , m 7 )] [l + Q 2 f F x {s, m f , m 7 )] + M« & X) , 

(97) 

where F\ is the standard electric form-factor regularized with a photon mass, see fig. |3]. 
We omit, for the moment, the magnetic form-factor F 2 ; this is justified for light final 
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fermions. It will be restored in the future. In F\ we keep the exact final fermion mass. If 
not stated otherwise, the four-momentum conservation p a + Pb = Pc + Pd holds. 

In the present work we use spin amplitudes for 7-7 and 7-Z-boxes in the small mass 
approximation m 2 e / s — »■ 0, rrij/s — »■ 0, see fig. [| following refs. f33|.|31"|. 



B,e B,f rp rpi 1 B,e B,f tji tt 

X 2 - M B + iY b X 2 /M b ^ 

Oi — — 

-QeQf [h a ,x c /bdp(M|, m 7 , s, t, u) - 5\ a -\ c /bdp(M|, m 7 , s, u, t) 

71 

where 

/bdp(Mb, ™^ *, «, t) = In (-) In ( j^-) - 2 In ( -) In ( 



u) \(t u y/ 2 J \uj \ M 2 



+ Li2| M_Li/^ 



M 2 i V s J V M| / \ M 2 ) \ M 2 

ms \ M 2 Z ) u \M 2 Z 



(99) 



M| = Mf - zM z r z , M 2 = m 2 , and the function / B dp is that of eq. (11) of ref. |4|. 
The standard Mandelstam variables s,t and u are defined as usual: s = (p a + Pb) 2 , t = 
(Pa — Pc) 2 , t = (pa — Pd) 2 - Since in the rest of our calculation we do not use m 2 /s — > 0, we 
therefore intend to replace the above box spin amplitudes with the finite-mass results^ 
that were given in ref. fl35| .) 



Now using eq. fl96| ) we determine 

(5; X) = <B (I; X) (l + 8^ t {s)) (l + 5^ t {s)) + ft« x ( P A ; X) (100) 



where 



$Virt( s ) = Q 2 e Fi{s,m e ,m^) - Q 2 e aB 2 (p a ,Pb,m 1 ) = Q 2 e --L e , 



tt2 

S vlrti s ) = Q 2 f Fi{s,m f ,m 1 ) - Q 2 f aB 2 (p c ,p d ,m 1 ) = Q 2 --L f , (101) 



L e = ln ( -^r ) + z'7r - 1, L f = ln ( ) + m — 1. 

V.m 2 y V m // 

Note that we departed in eq. ( |100| ) from the strict 0(a l ) by retaining the ^v\rti s )^vlrti s ) 
term, i.e., replacing the "additive" form 1 + 5y^ t (s) + Syj^ t (s) with the "factorized" form 

1 For the 7-7 box we use the spin amplitudes with the exact final fermion mass. It seems, however, 
that the 7-Z box for the heavy fermion is missing in the literature. 
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(1 + 8yJr t (s))(l + SyJ t (s)). The above does not need really much justification - it is 
obviously closer to the reality of the higher-orders, so the "factorized" form is preferable. 
The only question is whether the above method does not disturb IR-cancellations. It does 



not, as it is seen from the definitions of 5y\l t (s) and S v ^ t (s, 




Figure 4: Second order vertex diagrams. 

The IR-subtraction in using eq. fl96[) at 0(a l ) leads to the IR-finite CRbox- The 

above subtraction is equivalent to the following substitution 



/bdp(M|, m 7 , s, t, u) -> /bdp(M|, m 7 , s, t, u) - /m(m 7 , t, u) 



(102) 



where 

f m (m y ,t,u) = ^-B 2 (p a ,p c ,m J ) - ^-B 2 (p a ,Pd,m^) = In 



7T 7T \M/ WtuJ 2 V M 



(103) 



and the additional resonance factor exp ( — aBf (s)) in eq. fl96[) induces the additional 
subtraction in the 7-Z box part: 



/bdp(s, t, u) /bdp(s, t, u) - aBf{s) 



(104) 



see eq. ( [32]) for the definition of a£>f . 

Our 0(« 2 ) expressions for /3q are still incomplete. We base them on the graphs 
depicted in fig. [|. (In fig. ^ we omitted some trivial transpositions of the diagrams.) 
Following again the eq. (p6|), we obtain 



A® (5; X) = 25 (5; X) ( 1 + m e ) ) ( 1 + m 



(2)e 



<T)( 2 ) /p. 
-""Box U> 



X) (105) 



42 








+ 



AAA 




Figure 5: Missing second order diagrams. 



(2) 

In the present calculation we neglect the two-loop double-box contributions in ^g^, 
depicted in the first row in fig. |5| and vertex-box type of diagrams, see examples of diagrams 
in the second row in fig. |5] Pi. In fact we keep only the first order box contribution 3^bq X 
in our incomplete 0(a 2 ) type matrix element. 

Two remarks: in spite of the temporary lack of the above contribution we are not 
stuck because what we neglect is IR-finite! This statement is not so trivial as it may look 
because in the calculation without exponentiation neglecting such contributions would 
violate IR cancellations, and correcting for such a violation would be rather complicated 
and physically dangerous. Secondly, what we neglect is expected to be numerically small, 
of 0(a 2 L 1 ) and therefore it does not make much harm to our overall physical precision. 

Coming back to the 0(a 2 ) corrections to the electric form factor from the diagrams 
in fig. f|, they are well known since they were calculated in refs. [18, [37|-|39| and they 
contribute as follows 



t (2)e 
'Virt 



«N 2 

a 

7l 



4 2 T 



Li 



Li 



3 

32 
3 

32 




(106) 



In the above we kept terms of 0(a 2 L 2 ) and 0(a 2 L 1 ), and neglected the known p9|.[40 
negligible terms of O(a 2 L ). 

15 In fact the two- loop double-box contributions became known recently |3(|, so there is a chance to 
include it in the future. 
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3.4.2 Electroweak corrections 

In the not so interesting case of the absence of the electroweak (EW) corrections the 
couplings of two neutral bosons 7 and Z are defined in a conventional way: 

G ZJ = g ZJ - Xg Z /(7), Gf = g ZJ , \ = +,- = R,L, 
9v £ = Qe = h 9yj = Qf, 9a = °» 9aj = °» 



z>e _ 2T e 3 -4Q e sin 2 9 W ZJ _ 2T| - 4Q f sin 2 9 W 



z,/ _ f^/ ~ ^/ M " (107) 
16 sin 2 6V cos 2 6*vk ' yy 16 sin 2 9 W cos 2 0^ 



j A -1 ^ • 2/i 9/1 1 9 A 



16 sin ^cos 2 ^ 16 sin ^cos 2 ^ 

where T 3 is the isospin of the left-handed component of the fermion (Tj = —1/2, T e 3 = 
-1/2). 

The actual implementation of EW corrections is practically the same as in KORALZ 



T0|| • It goes as follows: The 7 and Z-propagators are multiplied by the corresponding 
hook-functions (scalar form-factors) due to vacuum polarizations 

1 G M 2 

# 7 -> # 7 x - — — , H z ^H z x 16 sin 2 9 W cos 2 9 W — z p E w (108) 

2 - ILy a QED 87rV2 

In addition the vector couplings of the Z get multiplied by extra form factors. First of all 
we replace 

Zj6 2T 3 -4Q e sin 2 # M / 2T e 3 - 4Q e sin 2 9 w F e EW {s) 



■ 

,/ ( 109 ) 



16 sin 2 9w cos 2 9w 16 sin 2 9w cos 2 9w 

2T 3 - AQ f sin 2 9 W _ 2Tf - AQ f sin 2 9 w F f EW (s) 



9V = ir„:„2n o „ => 



16 sin 2 9\y cos 2 9\y 16 sin 2 6 W cos 2 9\y 

where F EW (s) and F EW (s) are electroweak form factors provided by the DIZET pack- 
age |TJ]], which is part the ZFITTER semianalytical code || and correspond to electroweak 
vertex corrections. 

The electroweak box diagrams require more complicated treatment. In the Born spin 
amplitudes we have essentially two products of the coupling constants 

Z.e ZJ , Z.e x Z,e\ ( ZJ . > ZJ\ Z,e ZJ \ Z,e ZJ . x Z,e ZJ Z.e ZJ 

9x 9—x = (9 V - X 9 A )\3v + X 9 A ) = 9 V 9 V ~ X 9a 9y + X 9 V 9a ~ 9a 9a » 

Z.e ZJ 1 Z,e , Z,e\/ ZJ x ZJ\ Z,e ZJ , Z.e ZJ x Z,e ZJ . Z.e ZJ 

9x 9x = \9v - X 9a XflV - X 9 A ) = 9v 9v ~ X 9 A 9v ~ X 9v 9 a + 9 a 9 a ■ 

(110) 

In the above the following modification is done for the doubly-vector component 

z, zj _ > ^Tf-ST^QfF^-STfQ.F^ + WQfQfF^t) 
9v 9v (16 sin 2 9 W cos 2 9 W ) 2 ' 1 j 

where the new form factor F EW (s,t) corresponds to electroweak boxes and is angle de- 
pendent. The Born spin amplitudes modified in the above way are used also in the case 
of the presence of the single and multiple real photons, see next sections. 
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3.5 One real photon 



The discussion of the f3\ tensors corresponding to emission of a single real photon we start 
with the tree level case (zero virtual photons). The starting point is the well known 0(a l ) 
split amplitude for the single bremsstrahlung which we shall reconsider separately first in 
the case of the emission from the initial state beams (ISR) and later for emission from 
the final state fermions (FSR). This will be the "raw material" for obtaining /3i using 
eqs. (H). 



a 



1 



a 




b 




d 



d b 

Figure 6: ISR diagrams 



The first-order, 1-photon, ISR matrix element from the Feynman diagrams depicted 
in fig. reads 



M 1{/} g£) =eQ e v(p b , X b ) M 1 ia + ™ _ ~ & (h) u(p a , A.) 



where 



+eQ e v(p b ,X b ) f ai (ki) ^_2^ ^ M m u(Pa,K), 



B=f,Z 



;ii2) 



;ii3) 



is the annihilation scattering spinor matrix, including final-state spinors. The above ex- 
pression we split into soft IR parts^] proportional to (^±m) and non-IR parts proportional 
to fa. Employing the completeness relations of eq. ( |222j ) in the Appendix A to those parts 
we obtain: 

/ v -° htPa] [cd\ U IpaVlXai + 2k X p b ^ ^A^ipJ "° \_p h \ a J M 



/ , -° Y\ b p J Lp oiaJ 2k 1 p b ^ '- Af,CTlp ^ ^ Xa * [ ] ' 



;ii4) 



16 This kind of separation was already exploited in ref. |EM. We thank E. Richter-Was for attracting 
our attention to this method. 
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The summation in the first two terms gets eliminated due to the diagonality property of 
U and V, see eq. ( |229| ) in the Appendix A, and leads to 

M 1 <'>g£)=4?(*i)<BK]+' W (<*), 

r <" 8S) = E » K 1 ] m" RSH » Kl m. (115) 

The soft part is now clearly separated and the remaining non-IR part, necessary for the 
CEEX, is obtained. 

a 



c 



a 





b 



d b 

Figure 7: FSR diagrams. 



d 



The case of final-state one real photon emission (FSR), see fig. [7], can be analysed in 
a similar way. The first-order FSR, 1-photon, matrix element is 



M l{F} gS)=eQ / fi(p e) A c )^(fci; 

+ eQ f u(p c , A c ) M {F} - 



2kip c 
4 d + m- fa 
2kip d 



M v(p d , \ d ) 



;ii6) 



where 

M{F} = *e 2 n^(X) (G£j M Gf u , (117) 

B=j,Z 

is spinor matrix for annihilation scattering, including initial spinors. Similarly, the expan- 
sion into soft and non-IR parts for the FSR spin amplitudes is done in the way completely 
analogous to the ISR case 

gS)=*Jf } (*i)®(S)+w gy, 

w 85) = |^ E " KW] » M Kl - ^ £ » M K:?] v psa , (ng) 

4f } (fe 1 )-eQ / ^ ( f 1 ' Pc) -eQ / -^ (fcl ' Prf) 
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For the purpose of the following discussion of the remaining non-IR terms it is useful 
to introduce an even more compact tensor notation: 



U 



Ulf,, jh ® [ P x b b \ a a ] [ P x c P \ d ] = ©MM> 



;ii9) 



etc. For the "primed" indices we understand contractions, for instance 



0-'=± 



V 



(120) 



Using the above notation, the complete 0(a l ) spin amplitude for 1-photon ISR+FSR, 
coming directly from Feynman diagrams, with the explicit split into IR and non-IR parts, 
ISR and FSR parts, reads 



(lt\)(P-h) + Tl^ F} g£)(P) 



*g> © g; P - h) + r {1} g£; P - fcj + 4f 23 & P) + r {F} g£; P) 



2kp a ~^ 1 '^"^^^ 2k Pb 
r { F } (ll\;X) = ^LUiciv] X^X) - |^©Mi'](X) ^ 



(X)U { 



Via] 



[1] 

V[blV] ©ri'aceflPO 



(121) 



In the lowest-order the Born spin amplitudes 03 are defined in eq. fl4"4]), and we show 
explicitly as an argument the four-momentum X which enters the propagator of the s- 
channel exchange. Note that the formulas here differ by an overall sign from those of 
ref. H 



3.5.1 First- and second-order fii 

Now we employ the tree level, 0{a 1 ) variant of eqs. (|95|) getting the following results: 

3(1) (Ph. p\ _ rr i (Pkl . pW [(Pc+Pd + h? Ak(P-Y\ (122) 

Pi{f} Uoi ' = r (n U«n > + I — ( Pc + Pd ) 2 1 J ^ U) x ) > 

The "context dependent" reduced total momentum X (the total four-momentum in the 
resonance propagator, if present) is in the above definition uniquely defined as X = P — k\ 
in the case of ISR, and X = P in the case of FSR. In the general context of the CEEX 



amplitude of eqs. (|33| - p^ ), that is in presence of the additional "spectator" ISR photons 
in a given term, X is also defined quite unambiguously: X p includes not only k\ but 
also all additional ISR momenta in the process. For the pseudo-flux factor there is some 
ambiguity, however. In the presence of the additional "spectator" ISR photons it can be 
defined either as (p a + Pb — k\) 2 j (p a + Pb) 2 or (p c + p d + ki) 2 / (p c + Pd) 2 - We are free to 
choose any of them and we opted for the second choice (it seems to lead to more stable 
MC weights). 
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The one-loop level, 0(a 2 ) case of (3^ is quite interesting because this is the first 
time that we deal with the nontrivial case of the simultaneous emission of virtual and 
real photons. It is therefore instructive to write the formal definitions of f3\ following 
eqs. ( |96| ) and (|9^) in this particular case: 




Figure 8: One-loop corrections to single bremsstrahlung. 



What is presently available from the Feynman diagrams? For the moment we have at 
our disposal the amplitudes corresponding to vertex-like diagrams in fig. and we miss 
diagrams of the "5-box" type shown in the third (bottom) row in fig. More precisely, 
after applying the IR virtual subtraction of eq. (|123|) we expand in the number of loops, 
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keeping track of the initial/final state attachment of the virtual photon: 

+ aQ e Q f ^l hBoxFj g£;X). 

In the above expression the first term describes the already discussed tree level single 
bremsstrahlung, the next two correspond to vertex-like diagrams in fig. |8|, and the last 
one represents the "5-box" type diagrams in the third row of fig. p. In the present version 
we temporarily omit from the calculation the contribution to p[ from the last, "5-box" 
term which looks as follows: 

Pl{u},Box5 (\n> X ) = a QeQf m i{uj},Box5 (w X ) ~ S \l] ^BcL (ft X ) ~ S \l] ^Box gl X ) ■ 

(126) 

As we see the trivial IR-part, which we remove, is proportional to the ordinary box con- 
tributions already discussed before. We expect the above to contribute in the integrated 
cross section to be at most of 0(a 2 L 1 ), and in the resonance scattering it will be sup- 
pressed by the additional T/M factor. 

Limiting ourselves to the pure "vertex-like" diagrams of fig. ||, for one real ISR (u> = I) 
photon we obtain from the Feynman rules the following 0(Q 2 a 2 ) result 



+ ®g;X) 4?(ki)p$: t (s,aJ) 



(127) 



where 



p$: t (s,aj) = -Ql ~(V(s,a,P) + V(s,(3,a)), 

71 I 

V(s, a, $) = ln(a) ln(l - /3) (128) 
+ Li 2 (a)-^ln 2 (l-a) + |ln(l-a) 



2"" v " ' 2"" v " ' 2(l + (l-a) 2 ) 
and we use Sudakov variables 

2p a p b 2p a p b 
Let us make a number of observations concerning eq. ([12 

• The terms of 0{a i ) like \si Pvln\ 2 * n ^ ne cross section, although beyond 0(a 2 ), 
are not rejected, as it would be the case in the ordinary 0(a 2 ) calculation without 
exponentiation. They are included in the process of numerical evaluation of the 
differential cross sections out of spin amplitudes. (It is essential that they are IR- 
finite.) 
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• The term r^jySyj^. contributes 0(a 2 L 2 ) to the integrated cross section - one L 1 is 
explicit (from the virtual photon) and another L 1 is from the integration over the 
angle of the real photon. 

• The term ~ ln(a) ln(l — (3) contributes a correction of 0(a 2 L 2 ) to the integrated 
cross section with the double logarithm L 2 resulting directly from the integration 
over the angle of the real photon: 

/ ^ nP { vt(k){^ { J } (k)V)} ~ QW [ f ln(5) ~ QW In 2 -L 

J K Jm\js a m e 

• The other terms in ^1 contribute at most 0(a 2 L l ). 

• The FSR virtual corrections are included multiplicatively through the factor (1 + 
8ylrt(s)) and not additively like (l + ^y l > f (s) + 5y i ' ) / t (s)). This is our deliberate choice. 

• The subleading term a(l — a)/(l + (1 — a) 2 ) has in fact a more complicated spin 
structure than that of the Born amplitude (it should be restored in future). The 
unpolarized integrated cross section is however correct in 0(a 2 L l ). 

The analogous 0(Q 2 a 2 ) contribution for one real FSR (u = 0) photon is 
$® n g*; X) ee r {F} £*; X) (l + (l + 4»{s) + (a, ft) 

+ (i + (i + $£w) 

(130) 



where 



£') = ^/(ln(l - fi") + ln(l - /?')), 

„, 2kp d ~ 2kp c „„ &' ~ ? (131) 

a = , p = , a = — , p = — . 

2p c p d 2p c p d i + af + p/ 1 + & + P' 

In the above FSR amplitudes we keep only the LL part averaged over the photon an- 
gles, similarly as in EEX. This corresponds to present status of our CEEX amplitudes 
implemented in /C/C MC version 4.13, and we expect this to be improved in the future. 



3.6 2-real photons 

In the 0(a 2 ) contributions from two real photons are completely at the tree level, without 
virtual corrections (in the future 0(a 3 ) version we shall include virtual corrections to dou- 
ble bremsstrahlung in the LL approximation). The double bremsstrahlung is considered 
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in three separate cases, two ISR photons, two FSR photons and one ISR plus one FSR 
photon. The corresponding spin amplitudes will be given without any approximation, in 
particular we do not use the small mass approximation rrif/^fs « 1. The main problems 
to be solved will be 

(a) to write all spin amplitudes in a form easy for numerical evaluation, that is in terms 
of U and V matrices, 

(b) to extract (3 2 functions by means of removing IR-singular parts. 



3.6.1 2-real ISR photons 





b d b 2 b 2 

Figure 9: Feynman diagrams of the ISR double bremsstrahlung. 

The second-order, two-photon, ISR matrix element from the Feynman rules, see Fig. 
reads as follows 



TVf( 2 ) (PaPbk l k 2 . p jL u \ 
M 2{II} {\ a X b *iV2i^ ~ Ki- k 2 ) 



le 



£ ng(P -fa- k 2 ) (Gf„ ) M (eQ e ) 2 v( Pb , X b ) | 



G„„ — - faAfa, 



-2fap a -2k 2 p a + 2fak 2 



{j>a + m) - jjg 
-2k 2 p a 



Xk2 



, , (-^ b + m)+ fa 



{-j> b + m) + fa+ fa qb 
-2k lPb - 2k 2 p b + 2fak 2 e ^ 



(-A + m)+ fa B j^g + m) - fa 
Cih) 7TTZ G e,« 7TTZ faS k 2) 



-ZfaPb 



e '» ~2k 2Pa 



2)}u(p a ,\ a ). 



(132) 



We shall use eq (p5|) which in this case reads 

P% } KSSi P~h-k 2 )= gj£; P -fa- k 2 ) - g£; P -fa- k 2 ) *g>(fc) 

- (&; p - fc i - **) - $ 0) & p - h - k 2 ) 4?(ki)4?(k2). 

(133) 
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We shall proceed similarly as in 1-photon case, we shall isolate from the above expression 
the group of terms containing two factors (p'+m), then the group containing single factor 
(p' + m) and finally the rest. Such a split represents almost exactly the split in eq. ( |60D 
into contribution with two s-factors (double IR singularity), with single s-factor (single 
IR sing ularity) and the IR-finite remnant # 2 > which is our primary goal. In other words, 
we decompose 9R 2 {ii} severa l terms/parts, as described above, and we apply the 



IR-subtraction of eq. ( |133|) term-by-term. 

Let us discuss first the doubly IR-singular part proportional to two factors (p'+m). To 
simplify maximally the discussion let us neglect for the moment 2k%k 2 in the propagator. 
Using the completeness relations of eq. ( |222p and the diagonality property of eq. ( ^29| ) in 
Appendix A, we can factorize soft factors exactly and completely 



(Ja + m) (ji a + m) 

^ 2k lPa + 2k 2Pa ^ 1J 2k 2Pa **' 



2, 



* (, s (-A + m) (-fa + m) B 



(*0 "-^^ Gl + (1 - 2) \u( Pa , \ a ) 



{- j>b + m) qB (fa + m ) 
2k lPb ^ 2k 2Pa 

(ki,Pa) b a2 (k 2 , Pa ) b ai (ki,p b ) b a2 (k 2 , Pb ) 
1 2k lPa + 2k 2Pa 2k 2Pa 2k lPb 2k lPb + 2k 2Pb 



(ki,Pb) b a2 (k 2 , Pa ) 
2k lPb 2k 2Pa 



'1<->21 



(GfJ M «™(k)«£?(k), 



where the identity 

1 1 1 111 



2k lPa + 2k 2Pa 2k lPa 2k lPa + 2k 2Pa 2k 2Pa 2k lPa 2k 2Pa 



(134) 



(135) 



was instrumental. 

If we restore the terms 2k\k 2 in the propagator the corresponding analog of ( |134| ) 
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M DoubleIR leadg tQ 



A(2)Double rpk 1 k 2 ] _ TvfDoublelR rpfcifc 2 l M} 1 1. \ M}(h \ m IP 
- I S [l] S [2] A « + S [l] S [2] A & ) * b 



El, 



•«(*,) - < = -eQe^H ««(*,) - «g> = +e Qe b -^fl, (136) 

^ _ 2fcip/ + 2fc 2 P/ 1 ±2k x k 2 „ . , 

and the upper sign should be taken for f = a,b. Obviously / g( 2 ) Doublc i s IR-finite because 
of the A f factor. In the above we have introduced a more compact notation for s-factors. 
In addition from now on we shall use the following shorthand notation 

r if = 2ki-pf, rij = 2ki-kj, f = a,b,c,d, i,j,= l,2,...n. (137) 

The next class of terms which we are going to consider carefully is the one in which 
we sum terms with a single (tf+m), more precisely, let us include terms, which may lead 
to a single IR singularity (if k\ « k% or k 2 << fci), that is with (tf+m) next to a spinor, 
at the end of the fermion line: 

K$? R KiH = ^ E u b (x) {g% ) M 

B=-y,Z 



— Ha — ^2a + ^12 _ r 2a 



21 



j, x (-jfr + m) . s #i+ ^2 ^ 

+ ^(«i) C 2 {h) : G 

-ru -ru - r 2 b + r 12 

+ + (1 ~ 2)L(p 8> A.). 

— 7" lb —J" 2a 

(138) 

Using the compact notation, already introduced when (re) calculating single bremsstrahlung, 
we express ^ 2 Tfn in a form friendly for numerical evaluation, that is in terms of U and 
V matrices, 

rSinglc 
1 2{II} 



IVfSmglelR (pk 1 k 2 \ 
JVL 2{//| KXcriai) 



n -^[bl'][cd]U[l'la] ~ %>[b2>][cd]U[2>la] (a) n (ft) ^622'] ® [2'a] [cd\ + V[m']%$ [V a][cd] 
d\c^e [21 m%e *>m 



-ria - r 2a + r 12 1 J -ri a - r 2a + r X2 

r\ ( b ) ca U[ 2 '2a] . ~ V[blV] „ (a) . N 

eQ e Sft 1 + eQe 58[l'a][afl *[ 2 ] + U ^ 2 ) 

^2a ^16 



(139) 
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On the other hand the single-IR part to be eliminated is 

^l(l)[in2] ^l(l) [2] [1] — '[1] [2] ^ [2] 

e<5e23[bl'][ C d] — 1 ~ eQe ®[l'a][cd] S[ 2 ] + (1^2). 

v ''la ?~la / 

Altogether we get 

^Single (pkik 2 \ _ T^rSinglelR (pkik 2 \ _ o(l) C U} _ C U} 
^2{7/} UiTKrsJ ~~ JVL 2{//} KXcricri) ^l(l)[in2] ^l(l) [2] ^[1] 

/O a> ^[2'la] (a) . ^ (6) V[h21'] „ 

= -eQe _ na _ r2a+ri2 *[2] + * W _ ria _ r2a + ri2 



V[622'] V[b22'j ' 

+ eQ e si [ — [ — Q5 [2 'a]H + (1 <-> 2). 

11 v -ru - r 2a + H2 ~r<ib 1 



It is rather straightforward to see that the above is IR-finite. 

Finally, we have to include all remaining terms from eq. (|134|) which have not yet 
included in ^2{ii}- They are IR-finite (in the case of only soft photon energy) and they 
read 

${n } gS5) = ^ Yl U b(X) {G% ) m (eQ e ) 2 v(pb, A b 



g. (A + m)-„- w £ (fcjZ » < (fc 
-ria - r 2a + r 12 -r 2a 



— ^16 — rib — r 2b + r 12 



(142) 



+ IT" ^ 4( fc 2) + (1 - 2) 

'16 '2a 

Using tensor notation in the fermion helicity indices the above can be expressed in terms 
of U and V matrices as follows 

^Rest /Pfc!fc 2 \ / q y ^[ba'][cd]U[a'12"] ~ *%[bV][cd]U[Vl2"] ~ ^[b2'][cd]U[2'l2"} ~U[2"2a] 

+ (eQ ) 2 ^ W1 "^ ~~ ^[l"26']55[6'a][cd] + V[i"21']®[l'a][cd| + V[i"22']®[2'a][crf| (143) 

-ru -ru ~ r 2 b + r u 

■ ( rn \2 V [bll'] ~ ~t/[2'2a] n 9 x 

+ (eQ e ) — — 2?[i'2']M — : — + (1 ^ 2). 

— ^ 16 — ^ 2a 

The total ISR /32{n} is the sum of the three 

n (phik 2 \ _ ^Double (pk 1 k 2 \ , ^Single (pk±k 2 \ , SRest (pkjk 2 \ f-i aa\ 
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3.6.2 2-real FSR photons 

The case of final-state double real photon emission can be analysed in a similar way. The 
second-order FSR, two-photon, matrix element is 

<W } (K;^) = *'e 2 E n^p) (Ow (eG/) 2 *(Pc, Ac 

(^ c + m) + ffi j^c + m) + fa + fa c b 

%1 2k lPc ^ 2k lPc + 2£; 2 p c + 2fcifc 2 ^ (145) 
c b j-^d + m)- fa- fa * (-^d + m)- fa j, 



+ m) + fa (-^d + m)- # 2 
2£;ip c 2A; 2 p rf 



fa nrz /.^ ^rr- fa + (i - 2) ^(p d , A,) 



Similarly, the expansion into soft and non-IR parts for FSR spin amplitudes is done in 
the way completely analogous to the ISR case. The subtraction formula is now 

o(2) (Pk 1 k 2 . p\ _ m>(2) (Pkik 2 . p\ oW fpfci. p\ R {F} , i \ 
P2{FF} \\cn<72i r ) ~ JJl 2{FF} KX<ri<r 2 ^) \\crii ^ ) S a 2 \ K2 > 



ffL ( P x^P) 4f - $ 0) G; P) 4f } (*04f } (fe). 



(146) 



'1{F} 

First we obtain the contribution from terms with two (ft — m) factors 

A(2)Double rpfcifcal _ TyrDoublelR Tpfci^l A F K^(n (Pc + Pd + fa + fa) 2 

P2{FF} [Xcna 2 \ - Jn 2{FF} \_Xaia 2 \ S [l] S [2] *HH (pT+PdP 

- I AjSpjfipj +^S[i]5 [2 ] I^MH 

'"ic ''id 

«if>(**)=*S > (**)+4?(*i)=^ ) +<, 

and is explicitly IR-finite. The second group of terms with only one (jf — m) factor at the 
end of the fermion line is 

K^fT (Itltl) = ^ £ V%{X) (G^) [ba] (eQtf u( Pc ,X c 

B=j,Z 

(ft c + m) fa + fa nB B - fa- fa (-fa + m) 

m ~ m _ , _ , _ ^f,v + u f,» - , _ , - ?[i] r~ ?p] 

r lc r lc + r 2c + r 12 r ld + r 2d + r 12 ^2d 

r lc ^ 7,7 ^ + ^ n c f ' u — Wd f[2] ( } J ( } ' 

(148) 
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and it translates in the matrix notation (in fermion spin indices) into 

T^SinglelR (pk 1 k 2 \ _ 
JVL 2{FF} U o-i(T2/ ~~ 

~ (c) U[ c2 V] (A U[ c2 2>] 

=e v/ s fn ; : ^[ba][i'd] + ei^f Sm ; : -o[ba][2'd] 



-Vjl'ld] (d) ,. n m -V[2'ld] Id) 

■ ■ S [21 + eLJf U5 bo c2' ; ; S, 2] 

rid + r 2 d + r 12 ' 2| r ld + r 2d + r 12 [2J 



, pn Jc) n, ~V[2'2d] n U[clV] ~ Id) , r,N 

+ e V/ S [i] -0[H[c2'] + e V/ -O [60] [I'd] S [ 2 ] + I 1 ^ ^J; 

On the other hand the single-IR part to be eliminated is 
fl« «J F > 1 3(1) jn _ r mJF} , injn 

^l(0)[in2] ~T ^l(0)[2]*[l] ~~ '[1] *[2] "I" '[2] ^[1] 

= (+eQ e »[6a][l'«fl - e<5e^«[6a][cl'] ) s{ 2 ^ + (1 <- 2) 

Altogether we get 

ASingle (vk 1 k 2 \ _ M SingleIR ( P k 1 k 2 \ _ nO) (pki\ JF} [fc 2 ] _ A(l) /pfc 2 \ JF} [fell 
^2{FF} U<7i<7 2 ; — JVL 2{FF} U0-10-2J ^1(0) UoJ * L0-2J ^1(0) U<7 2 J * M 

_ pn Jc) J / U [c22>] U[c22'] \ ~ t/[e21'] ~ 
— e V/ S [i] ^ ( — — ] -0[6a][2'd] + — ^[ba][Vd] 



r 2c + n c + r 12 r 2c J r 2c + r lc + r 12 

- V WM -Vji'irf] ^ -V[2'id\ 

+eQf < "J5[ba][ci'} : : H : : »[6o][c2'l 

\ r id + r 2d + r i2 r M / ri d + r 2rf + ri2 

MM V (Pc+Pd) 2 ~ ) [1] [2] +(1 2) ' 



Finally we include the remaining terms in eq. 



M^V } g££) = *e 2 £ n^(X) (G^) H (eQ;) 2 fi(p C) A c ) 

F= 7 ,Z 

* fa * (P'c + m) + fa + fa B 

nil mi ; ; Lr /,i' 

, n B (-tfd + m)- fa- fa - fa 

+ ^f,v ; ; mi 

J rid + r 2d + r 12 11 r M 
fa riB ~ fa 



fa — G% fa + (1 «- 2) W A,), 
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which in the programmable matrix notation looks as follows 

o Rest (pk\k 2 \ Tv,fRest (pk\k 2 \ 

H2{FF} \Xa x a 2 ) — JVL 2{FF} \\<tkt 2 ) ~ 

I r\ \2 U[cll>>] U[i»2c']%$[ba][c'd\ + U[y>2V]?B[ba}[l'd] + ^ [l"22']93[6a][2'(f] 

—{eQf) 

ric r le + r 2c + r 12 

, / ^ \2 ~~ ®[6a][cd']^[d'12"] — ®[6a][cl']^[l'12"] ~ 55 [6a] [c2'] ^[2' 12"] ~ V[2"2d] (153) 

rid + r 2 d + ri2 r 2 d 

/ ~ no ^[cll'l „ —V[2'2d] ,-. „s 

+(eQ / ) 2 <8 H[1 ^ — ^ + (1 <- 2) 

The total contribution from double FSR real photon emission is 

n (Pk 1 k 2 \ _ ^Double (pk 1 k 2 \ , ^Single (pk 1 k 2 \ , SRest (Pk-Lk 2 \ (A^A\ 

P2{FF} \x<n* 2 ) ~ P2{FF} \X<na 2 ) + h>2{FF} {x*i<t 2 J + P2{FF} {x*ia 2 J ■ l i04 J 

3.6.3 1-real ISR and 1-real FSR photon 

As we have seen in the previous cases of double real emission most complications are due 
to simultaneous emission from one fermion "leg" . The case of one real ISR and one real 
FSR photon is easier because there is at most one photon on one leg: 

KS£2S&; p ~^ = ^ 2 £ u b (p - h) 

B=-y,Z 

eQ e v( Pb , A*) (Of /**™-/ 1 ffo+ fa h G*^ u( Pa , X a ) (155) 



eQ f u(p c , Ac) ( Gj„ — fa + ^ ^ G f v ) v( Pdl X 



and the subtraction formula is now 



(i) 

L{F} 



(156) 



The simplicity of this contribution is manifest in the fact that $2{if} is obtained by 
simple subtraction (omission) of all terms proportional to one or two — m) factors 

km( p xt\ll;x)= ie 2 Yl n%(x) 

B=-y,Z 

eQ e v( Pb , X b ) [cl^ fa + fa A- Gl) u( Pa , X a ) (157) 

V '"la ?"lb / 

eQ f u( Pc , X c ) (g% ^ fa + fa h. G b\ v(pii Xd) . 
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In the computation- friendly matrix notation it reads 
/Wl(K 2 2 ;*)=*e 2 £ IL%(X) eQ e eQ f 

B=-y,Z 

■<B \ —U[Vla] V[blV] (n B \ \ ( fr<B \ ~~V[2'2d] U[ c22 '] , B 



/ mB \ i 1 xu \ i i ul1 J snB \ i inB \ \* ■** , ~ i^b \ 

X l G e ,/J[M'] 1 (k e ,/J[l'a] 0/,J[c2'] 1 0/,J[2'd] 

-ri a -rib / \ r 2d r 2c J ' 

eg e eg/ «[6r][ C 2'](Aj + — — K [bV][2 , d] {X) 

—ria r 2 d r 2c —ri a 

+ — — »[l'a][c2'](Aj— + — rg[Va][2>d]{X) 

—r ib r 2d —rib r 2c 



(15* 



4 Relations between CEEX and EEX 

Having shown the CEEX and EEX schemes in a detail, we would like to compare certain 
important /interesting fearures of both schemes in a more detail. In particular we would 
like to show how the two examples of the EEX scheme can be obtained as a limiting case 
of CEEX, and to show the exact relation among /3's of EEX and /3's of CEEX. From 
these considerations it will be clear that the CEEX scheme is more general than the EEX 
scheme. 



4.1 Neglecting partition dependence 

Let us first examine the interesting limit of CEEX in which which we drop the dependence 
on the partition index X p — > P, where P = p a + Pb, for example. Note that it is not in 
the EEX class. In this limit, in the simplest case of the O(a ) exponentiation we have: 

•y-2 n n 

j^e^Y *&*JII4r } =► ^ n (4? } +4i } )' ( 159 ) 

per cd i=i i=i 

because of the identity 

n n 

En4r , -n(4r+4i > )- 

peV i=l i=l 

The relevance, advantages and disadvanteges of this scenario were already discussed in 
Section ??? Is this realized in /C/CMC????? Note that in the above transition we keep the 
ISR®FSR interference contribution. 



4.2 Neglecting IFI 

The second important case we would like to discuss is the case of the very narrow res- 
onances, when the ISRcgESR interference contribution to any physicaal observable is so 
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small that it can be neglected whatsoever. This corresponds to a well defined limit in the 
CEEX scheme. In this limit, in the simplest case of the O(a ) exponentiation we have: 
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What we did in the above transition, we have neglected the ISR®FSR interferences en- 
tirely by dropping non-diagonal terms p ^ p' in the double sum over partitions, and we 
have replaced the resonance-type form-factor by the sum of the traditional YFS formfac- 
tors for ISR and FSR (no interferenence) . In this way we have got the O(a )EEX which 
at this order is identical to O(a )cEEX- At the C(a r )cEEX r = 1,2, in order to get from 
0{a r ) ceex to 0(a r )-EEX we have in addition to truncate /3's down to (3\ as will be shown 
in the next subsection. 

The 0(a r )EEX t — 1, 2 neglecting the ISR®FSR interferences was used in YFS2/3 |7],[| 
of KORALZ [|K]] and it is well justified close to Z resonance position at LEP1, see also 
relevant numerical results in the next Section. At LEP2 the above approximation cannot 
be justified any more. 



4.3 Relation among (3's for EEX and /3's of CEEX 

For the sake of completeness of the discussion, it is necessary to find out the relation of 
the /3's defined at the amplitude level to the older EEX/YFS /?'s defined at the level of the 
differential distributions. Let us suppress all spin indices, understanding that for every 
term like |...| 2 or !R[AB*] the appropriate spin sum/average is done. The traditional /?'s 
of the YFS scheme at the 0(a 2 ) level are 
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where subscript |( a r) means truncation to 0(a r ). Now for each 97tn n+ ^we substitute its 
expansion in terms of /3's according to eq. ( p8|) getting the following relation 



/^° = IOT(rf), ' = 0,1,2, 
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/^(fc) = |/^)| 2 + 2R[$ , W(A0)*] (a H.i )) ' = 0, 1 
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(162) 
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As we see, the relation is not completely trivial; there are some extra terms on the RHS, 
which are all IR-finite. From the above exercise it is obvious that /3's are generally more 
complicated objects than the /3's and that for example the inclusion of the spin density 
matrix formalism into the /3's would be a quite nontrivial exercise - the great advantage 
of the CEEX scheme is that this is done numerically. It is also seen that in the (3q and f3\ 
some higher-order virtual terms are unnecessarily truncated, which probably is worsening 
perturbative convergence of the EEX/YFS scheme in comparison with CEEX. The above 
formula shows in a most clear and clean way the difference between the EEX and CEEX 
exponentiation schemes. 
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5 Semi-analytical approach 



In this section we shall present results of semi-analytical calculations which reproduce 
in the 0(a 2 ) prag , or even in the 0(a 3 ) prag , certain selected results, mainly integrated 
cross-sections, of the Monte Carlo calculation. 

The semi-analytical approach in which one integrates over the phase space analyti- 
cally, often leaving the last one- or two-dimensional integrations for numerical treatment 
(usually non Monte Carlo), is the oldest one. Four decades ago there were no computers 
powerful enough even to dream about the numerical integration over the complete mul- 
tiphoton phase space. Even now, in spite of proliferation of the MC programs, the non 
Monte Carlo semianalytical programs are still very popular and useful, especially pro- 
grams used to calculate the total cross section and charge asymmetry for the fermion-pair 
final state, near the Z resonance, like ZFITTER and TOPAZO || 41]. Semi-analytical 



programs have certain advantages over the MC programs - they are generally faster in 
terms of computer CPU time and are therefore better suited for fitting input parameters 
of the Standard Model, like the Higgs massQ. Nevertheless, semi-analytical calculations 
have also two long-standing important disadvantages: 

(a) They are able to provide predictions for rather very primitive or absent experimental 
cut-offs. In practice they always have to be used in combination with the MC event 
generators anyway. The MC is used to remove or "straighten" the real experimental 
cuts to be closer to these which are practically implementable in the semianalytical 
calculation. Obviously this introduces additional systematic errors. 

(b) They are prohibitively complicated beyond the three-body final state, that is they 
are relatively easy up to 0(a x ), (single photon emission in the fermion pair produc- 
tion). In the collinear approximation (structure function method) one is able to add 
the effects due to emission of the second photon and further photons; however, this 
can improve the precision only within the leading logarithmic scope and makes the 
introduction of the realistic cuts even harder. 

The other important role of the semianalytical calculations is to provide the numerical 
tests of the Monte Carlo programs. They are typically used to check technical correctness 
of the phase space integration, the so-called technical precision, and also correctness of 
the implementation of the SM matrix element. In the following we shall see examples 
of both kinds of tests. In particular, we shall see the test of the technical precision of 
/C/CMC at the 2 • 10~ 4 level based of the semianalytical formula obtained in this section, 
in the case of a single kinematical cut on the total energy of all photons. 

The role of semianalytical calculations as a test of the Monte Carlo programs is im- 
portant but limited. One can easily imagine the situation in which the numerical problem 
shows up not for simple kinematical cuts, for which the "calibration" with help of the 
semianalytical program has been done, but for more realistic and complicated cuts. It 



17 It is definitely possible to fit input SM parameters with help of the Monte Carlo event generators, 
as it is currently done in the VK-mass measurement in LEP2. 
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is therefore always true that the ultimate test of the MC calculation is always the com- 
parison of two MC programs, because it can be done for arbitrary cut-offsQ This recipe 
may look often prohibitively expensive in the effort required to realize it; however, at 
the sub-permille precision level the amount of work required for the realization of the 
semianalytical formulas and testing the semianalytical program is probably comparable. 
So altogether, for the sub-permille precision prediction the approach with two indepen- 
dent MC programs seems to be the most economical one. In any case the semi-analytical 
calculations will be always very useful especially when the precision requirements are not 
excessive and we do not deal with the observables involving complicated experimental 
cuts. 



5.1 Inclusive exponentiation, IEX 

An important ingredient in many semi-analytical calculations is the "exponentiation". 
The meaning and the technique of exponentiation in the context of the semi-analytical 
calculations is however different from the exclusive exponentiation discussed in most of 
this paper. Let us elaborate more on what the exponentiation really means in the semi- 
analytical approach. 

As already discussed in ref. ||, in the typical semianalytical approach one is practicing 
what we call an "ad-hoc exponentiation" or "naive exponentiation" procedure in which 
one takes the QED finite-order, let us say 0(a) or 0(a 2 ), analytical result for a certain 
one- or two-dimensional inclusive distribution and this result is "improved" by hand, in 
such a way that the soft limit (the limit in which hard photons are eliminated) in the 
resulting distribution conforms to the Yennie-Frautschi-Suura work M. The well known 



examples of the ad- hoc exponentiation are presented in refs. ||42|,[43|1 and later in ref. [44 
for the initial-state bremsstrahlung in e + e~ annihilation. There are also many other 
examples of the ad-hoc exponentiation, including calculations for the deep inelastic and 
Bhabha scattering processes. The ad-hoc exponentiation procedure may get improved 
gradually by taking into account higher-order effects. For example, the 0(a) procedure 
of ref. [^2| was extended to 0(a 2 ) in ref. [43|] and later to 0(a 3 ) in ref. ||44j[ . The problem 



is that the procedure of ad-hoc exponentiation is essentially rather art than science - 
one may regard it at best as a "by hand interpolation" between two kinds of analytical 
formulas - one valid in the soft photon limit and another in the hard photon limit, see 
ref. j|] for more detailed discussion on this interpretation. Since this approach is not 
systematic, it is therefore difficult to estimate the uncertainty of the results and it has to 
be "reinvented" for each perturbative order and for each inclusive observable again and 
again. 

The self-suggesting question is therefore: Is there any better and more systematic way 
of re-formulating the ad-hoc exponentiation for any inclusive distributions in the analytical 
form? I would be also desirable to find a direct connection to the exclusive exponentiation 
YFS exponentiation (of the EEX or CEEX type) which is discussed and implemented in 

18 This kind of test was for instance done for the first modern 0(a 1 ) Monte Carlo event generator 
MUSTRAAL of ref. (HQ, witn tlic VC1 T m S h precision at that time of 1%. 
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this work. There is an obvious hint in which direction to go. If all soft photons are soft, 
then we know exactly the analytical formula for the multi-photon phase space integral, 




) 



already obtained in the original YFS paper ||. So why not include hard photons in 
the game? Let us therefore define the "YFS inclusive exponentiation" as a result of the 
analytical phase space integration of the distributions of the YFS exclusive exponentiation: 



In this way we have got a clear and clean connection between the YFS exponentiation (as 
implemented in the Monte Carlo) and the "YFS inclusive exponentiation" - the connection 
is simply the analytical phase-space integration! As a result, we do not need any obscure 
"recipes" any more and what we only need to know is how to integrate (analytically) the 
phase-space. 

The above looks promising but the YFS exclusive exponentiation involves non-trivial 
integrals over the multiphoton phase-space, this is why it is implemented in the form of 
the Monte Carlo integration/simulation numerical program. The relevant integrals over 
n real photon phase space do not seem at first sight to be treatable analytically at all. 

The situation is not so hopeless, however, as it may look at first sight, and in the 
following we shall present the solution. We start, as promised, from the full YFS ex- 
pression, the same as in the Monte Carlo and are able to do "the impossible" — that is 
to perform the phase-space integral analytically. We calculate the phase-space integral 
approximately. We shall do it, however, in such a way that in the approximate method 
becomes exact in the soft photon limit, The soft photon contributions are there integrated 
exactly and only the remaining "non- infrared" contribution will be calculate using approx- 
imate methods, typically the leading-logarithmic (LL) collinear approximation. The LL 
approximations in non-IR parts may concern both the phase-space parametrization and 
the matrix element^. 

The profit from the above approach is two-fold: Contrary to the traditional ad-hoc 
exponentiation we gain, for a given exponentiated inclusive distribution, a clear and clean 
connection between the YFS exponentiation (as implemented in the Monte Carlo) and 
the "YFS inclusive exponentiation" , we do not need any obscure "recipes" any more. It 
means that the resulting inclusive exponentiation (IEX) is now a systematic order-by- 
order procedure, this feature is simply inherited from the exclusive YFS exponentiation. 

As already stressed, the inclusive YFS exponentiation will never fully replace the YFS 
Monte Carlo because it is possible to deal analytically with only a very limited number 

19 Let us note that the LL evaluation of the phase-space integral was already employed to some extent 
in the original YFS work ||. At that time, because of lack of fast computers, it was the only accessible 
method. 



YFS inclusive exponentiation 



Analytical integration of 
YFS multiphoton integrals. 
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of the distributions, without cuts or with very simple cuts, while in the Monte Carlo one 
may calculate an arbitrary distribution in the presence of the most complicated cuts. 

In the following subsections we shall show explicitly the analytical integrations leading 
to 0(a 2 ) prag IEX results. We shall typically compare the Monte Carlo with EEX matrix 
element and IEX formulas, both in the 0(a 2 ) prag class. Their difference will be then nec- 
essarily of 0(a 3 ) prag , i.e. up to factor 10 smaller - quite a strong test of both calculations. 
On one occasion, we shall go to a more difficult level of the 0(a 3 ) prag , in which case the 
difference between MC and IEX is of order 0(a 4 ) prag . 

Finally, let us note that the set of IEX formulas presented in this section was used over 
many years as a basic test of the precision of the YFS2 [j7| and KORALZ/YFS3 P,[l0 



programs. Some of the IEX results were already shown in ref. |7j and ||. Most of them 
are, however, shown here for the first time. In the mean time the analogous set of IEX 
results was obtained and published for the t-channel dominated process |13|]. In fact 
ref. |13| shows an even more sophisticated case in the 0(a 3 ) prag class than the 0(a 2 ) prag 
results presented here. Using the experience of ref. ]TB|] it would be definitely possible, 
for the s-channel process of this paper, to upgrade systematically the calculation of this 
section to 0(a 3 ) prag , both for ISR and FSR. 



5.2 Semi-analytical formulas for ISR 

We shall start the construction in IEX expressions with the ISR case, first showing the 
basic techniques working out the example with the O(a ) EEX matrix element. In this 
case the multiphoton differential distribution is just the Born cross section times real soft- 
factors. While for the other IEX formulas the phase space will be integrated basically in 
the 0(a 2 ) prag , the case of the of C(a°)EEX we shall make more effort and do it in the 
0(a 3 ) prag , like in ref. [O]. Let us attract attention of the reader to the fact that we have 
the matrix element in the C(«°)eex and the phase space integration is in 0(a 2 ) prag or 
0(a 3 ) prag . There is no contradiction in this, as we shall see in the following. 

5.2.1 Baseline high precision results for 0(q;^)eex 

The complete 0(a 2 ) prag calculation/exponentiation according to the rules layed down 
in the beginning of this Section will be rather involved, let us therefore illustrate our 
calculational methods with the simplest possible example. Even this simple example 
features some nontrivial technical features and we shall therefore present two versions of 
the calculation. 

The basic example discussed in the following is the O(a ) initial-state YFS inclusive 
exponentiation. In the master equation (^) we set the charge of the final fermion to zero, 
Qf = and we replace the sum of /?'s with the O(a ) version of /3q proportional to the 
Born differential cross section 

$ 0) (fc,(h) = J- -^rt> + fe) 2 ,tf), Pf = (l-4mj/(g 1 + g 2 ) 2 ) 1/2 , (163) 
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where the normalization is such that 
d 3 q 1 d 3 q 2 



[ d -f d -f P\X - qi - g 2 ) flfVfc) = aWfa + q 2 f). 



The initial-state C?(a°) YFS formula reads 
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Integration over the final-state fermion 2-body phase space is done trivially leading to 
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where the essential multiphoton integral 
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is the main object of our interest. Note that we have split Yr(e) = 7 hie + ^y^s- 

The QED matrix element, beyond the soft photon integral, is in this simplified case 
totally absent. The inclusive YFS exponentiation, as defined above, amounts to calcu- 
lating analytically the multiphoton phase-space integral for po(v). As explained above, 
we shall do it in 0(a 2 ) prag , but we shall keep the proper soft limit undestroyed. Let us 
note first that in the soft limit v — > the function po(v) coincides with the soft integral 
of eq. ( |5.1| ), i.e. p(v) — > f(j,v). Since the most singular part in this limit is known we 

prag 



isolated it and we expect the 0(a 2 ) nrn n result to be in form 



Po(v) = f{l,v)(l + wyf 1 (v)) 



(168) 



where fi(v) is nonsingular. How does one find the function fi(v) ? Let us inspect the 
difference 



do(v) = p(v) - /(7,v) 
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5(l-v-\p 1 +p 2 -Y,k j ) 2 ) -5(v-- s ( Pl +p 2 )- (j^ k i)) 
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This new object has rather interesting properties. First of all, the Ola 1 ) integrals cancel 
exactly and the first nontrivial integral is of 0(a 2 ). This second-order integral is not, 
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however, infrared divergent! According to our rules we are therefore allowed, without any 
danger of spoiling the soft limit, to calculate it in the leading- logarithmic approximation. 

Let us present now our first of two methods of calculating po(v). In the LL approxi- 
mation we replace collinear singularities in the photon angle ?9 7 = 0, ir by 5-like peaks 
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I 7 5( Ci _ i) + I 75 ( Q + i) 



where 

Pe = (1 - 4mg/s) 1/2 , q = cos 6^, Si 
and using the above LL substitution we get 
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Two immediate remarks are in order: Out of four terms in the product [S(c\ — 1) + 8(c\ + 
l)][5(c2 — 1)6 (c2 + 1)} only two contribute, these with two anticollinear photons, c\ = 
1, c 2 = —1 and Ci = — 1, c 2 = 1. The result of integration depends critically on the careful 
regularization and for this reason we keep explicitly the e infrared regulator. A quick 
calculation gives a zero value for the integral. The very similar phenomenon is present 
in the calculation of /(7, u) where a naive calculation up to second-order gives the t? 7-1 
factor only. The remaining F( 7 ) = 1 - f27 2 + - factor comes from careful consideration 
of the k° > e^/s/2 condition for two photons. With our proper regularization we obtain 



dn(v) = 



7 2 ln(l — v) 



(172) 



which is finite in the v — > limit. 

Now comes the second calculation method which will be often employed in the fol- 
lowing. In this variant we take into account the influence of additional soft photons (in 
addition to the two hard ones). They do not change the second-order result but provide 
the proper infrared regulation replacing the former e regulator. The LL treatment of the 
phase-space will be a little bit different. Starting from eq. ( |5.1|) we split (in CMS frame) 
the photon integration into its forward and backward hemisphere parts 
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and after changing the summation order we get 
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where P = pi + P2, K + = kf , if - = ^ kj . The above sum of integrals factorizes 

i=l j=l 

into two sums. Each of the sums can be evaluated exactly leading to the following simple 
convolution 

/(%«) = / dv + dv-6(v-v--v + )f(l,v + )f(l,vj, (174) 

This identity holds for the integration result anyway, but we have also obtained it through 
the direct phase-space integration. So far all calculations were exact and we only reorga- 
nized the phase-space integration which will be useful in the next step. Let us consider 
the do(v) difference again 
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As before, the whole integral is finite in v — > limit and it gets the first non-zero con- 
tribution in the second-order. From the previous exercises we know that the essential 
second-order leading-logarithmic contribution comes from two anticollinear photons - this 
is why we divided photon phase-space into two hemispheres. Now, the LL approximation 
is realized by substituting in the first 5 

K ±fl = (K ±0 ,0,0,±\K ±0 \ ). 

Note that, contrary to the previous calculation, we did not modify the S factors, we did 
not introduce collinear <5's in the photon angular distribution and we keep infinite numbers 
of photons. In spite of the apparent increase of the complication level, the integral reduces 
to a nice form 



dv+dv. 
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which is calculable analytically! Neglecting terms 0(7 3 ) we obtain 

d (v) = -ff^-f^' 1 \ 7ln(l - v). (177) 

Note that since in the present variant of the calculation we have treated soft photons 
more friendly we recovered the proper soft factor f(j, v) as a factor in the solution. 
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Table 1: The comparison between ZC/CMC and IEX 0(a 2 ) prag formula of eq. (181) for the constant 



Born cross section. (200GeV?) The difference between KXMC in EEX mode and semianalytical 
formula divided by Born is plotted with dotted line, as a function of the w max cutoff on the total 
energy of all ISR and FSR photons. We also include the integrated cross section divided by Born, 
and multiplied by factor 10~ 2 , as dots for IEX and a line for MC. 



Summarizing, the 0(a 2 ) prag phase-space integration result is 
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and the corresponding cross section reads 
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The above integration methods provide us with the 0(a 2 ) prag phase-space integration 
result for any of the /3 contributions as listed in eqs. (|). For example the contribution 
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from (3^ reads 





= F(-f)e 8YFS 7w 7 - 1 (l + 5\ 1] + Sf ] ) (l - J 7 m(l 
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From now one we shall not restrict ourselves to O(a ) prag EEX matrix elements, but 
rather consider the complete EEX-class 0(a 2 ) prag EEX matrix elements as defined in 
Section 
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Figure 10: The comparison between OCMC and IEX 0(a ) prag formula of eq. ( |181|) for the s- 
dependent Born cross section at 189GeV. The difference between OCMC in EEX mode and IEX 
formula divided by IEX is plotted as a function of the v m & x cutoff on the total energy of all ISR and 
FSR photons. 



The practical significance of IEX formula of eq. |180| is rather important. The biggest 
terms neglected in it are of 0(7 3 ) and 0{a r y) and we expect them to stay below 0.1%. 
(This will be true provided they are no extra enhancement factors, see discussion below.) 

— (2) 

In other words we expect for the contribution in the EEX matrix element in Section 



FJ the result of the Monte Carlo phase-space integration will agree with the formula ( |180| ) 
to within about 0.1% for an arbitrary cut v max . 

Let us check the above conjecture with the numerical exercise. In the numerical test 
we shall already include at this moment not only ISR f3^ contribution of eq. (|180Q , see 
also Table FJ, but also the analogous FSR (3^ contribution which will be calculated^ in 

20 We could present results of the numerical tests (which we have done) for ISR alone. They look 
however very much the same like as simultaneous ISR and FSR so wc decided not to present their figures. 
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the next sub-section, see eq. (|204|) and Table |||. We consider the total cross section with 



the cut on the total photon energy defined by v max as follows 

/Unax — "max) 

dv j( du a f Born (s(l - «)(1 - v)) p®p(v) pf-, 2) (u) (181) 

In order to get a clearer picture about the magnitude of the discrepancy between EEX MC 
and IEX formula we use the artificially flat Born cross-section 0B Orn (s(l — u)(l — v)) — > 
<J Born( s ) i R b°th- Results of the comparison are presented in Figure |l|. Following our 
expectation the difference is well below 0.1% for the entire range of the photon energy 
cutoff w max . 

The situation does not look as good when we switch-on the s-dependence in Born cross- 
section. In Figure |T0|(a) we see the relevant comparison. At the CMS energy of 189GeV 
the position of the Z radiative return is at v = 0.75 and we clearly see a worsening there 
with respect to the previous case in Figure [I] where the discrepancy is now almost 0.2% 
(0.4% in terms of 0B Orn ). The situation is even more dramatic in the last bin which 
corresponds to f max = 1 — 4m^/s and here the discrepancy among 0(a 2 ) prag IEX and 
0(a 2 ) prag MC EEX is -2% of the total cross section, that is -7% in terms of the Born 
cross-section! This is, of course due to the Z resonance and 1/s behaviour of the Born 



cross-section at very low s (especially for the case of the p channel shown in Figure 1C). 
In order to be sure that the above effect is not due to some technical problem in the MC 
integration we had to improve our IEX formula and upgrade the analytical phase space 
integration for ISR to the level of 0(a 3 ) prag . The comparison with the 0(a 3 ) prag IEX 
for the same EEX 0(a 2 ) prag MC we show in Figure |10|(b). Now the difference is reduced 
to less than 0.1% everywhere, an in the last bin it is reduced from -2% to +0.2%, as 
expected. The additional terms of the 0(L l a 2 ) and 0(L 3 a 3 ) are shown in Table ^ at 
the end of this Section. We do not show the details on the phase space integration which 
provide these two additional terms. The method is generally rather similar to the one 



used in this Section and in ref. 13 
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5.2.2 Beta-bar-one, (3\ 

In the following step our aim is to calculate analytically the ISR contribution to the total 



— (2) 

cross section from j3\f as given by 



n=0 U - J <ll <l2 J j=1 K j 

e y ( n,) ^2ff){X )Pl}P2) q 1} q 2 , k^/S^k^ ( Pl+P2 - qi -q 2 -J2 k 3 ) 
j= i V j=1 / 

n=0 "-7 ?i <? 2 7 i=1 ^ 

1 - 9(fi/; /c))<5 (4) ( pi +p 2 - q\ - q2 - k - J^fcj )j3$(qi + q 2 ,Pi,P2,qi,q2,k). 
V j=i / 



;i82) 



We start again from the EEX 0(a 2 ) prag matrix element for the initial-state bremsstrahlung 
and we shall perform the phase-space integration also in 0(a 2 ) prag . We integrate first over 
final-state fermion four-momenta 

i ft 9 2 1 } (183) 

= B[ 2 \ Pl ,p 2 ,k) a Bom ((g 1 + g 2 ) 2 ) 

where 

^ 2) (Pi,p 2 , fc) =^7T^rr-,w e (aj)(i + A5 1) (d,/5))i{(i - «) 2 + (i - M 

4vr 2 {k Pl )(kp 2 ) 2 1 J (1 8 4) 

-5 / (p 1 ,p 2 ,fc)(l + ( 5i 1) ) 



and obtain 



B?\ Pl ,p 2 , k) a Bo ™(s(l - v ))5{v-- s (P-J2 ^ - kf) (185) 



dvpf\v)a Born (s(l -v)) 



(2) ~ 

In the calculation of p\ we could follow the first of the methods employed for [3 . Let 
us describe it briefly without going into details of the calculation. We calculate the first 
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two non-zero integrals 0(a) and 0(a 2 ). The first one 0{a) has to be calculated keeping 
both the leading O(La) and the subleading term O(L a). This can be done following the 
well known 0(a) analytical calculations [15]. The 0(a 2 ) integral with two real photons 



can be treated in LL approximation, i.e. keeping only 0(L 2 a 2 ) terms. This can be done 
introducing collinear peaks in the photon angles as demonstrated in the case of (3q. Both 
integrals are connected due to infrared regulation with e. The first one is proportional to 
e 7lne ~ 1 + 7 hie and the term 7lne from the first one cancels the infrared divergence in 
the second (independently of the LL approximation). As in the case of j3o one has to pay 
attention to the subtle "edge effects" in the e regularizationF]. 

Let us describe in detail the second method in which soft photons provide the conve- 
nient infrared regulation. The main 0(a) contribution comes from the configuration in 
which we have k° ~ Vyfs/2 and one or more soft photons. This part has to be calculated 
exactly in 0(a). We split, as before, 



P 



( 2 ), 



ft 



(2), 



4 



(2), 



;i86) 



(2) 

in such a way that d\ (v) is vanishing in 0(a) - it can be therefore calculated in second- 
order LL while fx(v) is simple enough to be calculated exactly in the 0(a). We define 



fp"-±h I rif^ 



T d% - 



5 { v ~ -/ ■ E k i + k))B[ 1 \p 1 ,p 2 , k) = e s ^s j ^/( 7 ,« - - g P ■ k)B[ 1) (p 1 ,p 2 , k) 

(187) 



where 



B[ 1) (p 1 ,p 2 ,k) 



a 



^(H^ W ' {&J)1 2{ (1 " &? + (1 - « 2 ^ - Sl(fM (188) 



(2) 

The remarkable feature of is that we could integrate over spectator photons exactly. 
Note that the e regulator has disappeared from the k— integral. In the next step we 
integrate exactly over photon angles following the old O{o) calculations and we are left 
with the single integral over the photon energy x = 2k / y/s, with the strongest singularity 
(v — x) 7_1 being regularized nicely by soft photons 



/f 2 V) 



3 Syfs 



F( 7 ) / dx 7 (i> — x) 1 1 7 



1+ r. 



d Syfs 



-P( 7 ) 7f 7 



1 1 1 

1 H — v 7f 

2 2 . 



;i89) 



0( 7 3 ). 



21 Generally, the calculation for /3i is more difficult than for f3o and /?2 because this is the only case in 
0{a 2 ) where we deal with simultaneous real and virtual photon emission. 
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Figure 11: The comparison between ZC/CMC and IEX 0(a 2 ) pra g formulas for the integrated cross 
section as a function of the cut-off parameter u max on ISR and FSR photons. Presented are /C/CMC 
(solid line) and IEX (line with open circles) results for: (a) ISR (3$ and the FSR (3$, (b) FSR (3$, 
and the ISR j3 0I , multiplied by factor 0.1 (in order to fit into the scale). The difference between 
/C/CMC and IEX is shown as a dotted curve. Center of the mass energy is 189GeV. Final state 
fermion is muon. 



(2) (2) 

Now we shall calculate the remaining part d\ of p\ . Since it vanishes at 0(a) we 
may calculate it in LL approximation. Although strictly speaking it is not necessary, 
we treat photons gently (as in (3 Q example) such that we do not use the crude collinear 
approximation. As before, we split the photon angular integration into forward and 
backward hemispheres and we integrate immediately over the final fermion momenta 



T) r> 1 



6<-k/2 



n 



i=l 



d 3 kj 



MK)Q IK 



n 



e 4 >7r/2 



j ~ 1 OjKir/2 1 



5 v 



-5(v 

+ 5(v-l + -{P 

\ s 



(P - k - K + - K~f 



2P ■ (jfe + K + + K- 



(190) 



B?(pi,P2,k) 



K~ 



i yB? ) ip 1 ,p 2 ,k)-Bi i) ip hP2 ,k)}^. 
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Using the collinear replacement K ± = (K ±0 , 0, 0, ±|if ±0 | ) in S's allows us to integrate 
over spectator multiple photons 



df\v) = I dv+ I dv- I /( ±v+ )f[J.,v 



d 3 k 



k° 



7 



0<n/2 



(191) 



5^v — 1 + (1 — x — u+)(l — — (jfu — a; — u+ — f-j 

+ <j(v-l + (l-a;- v + )(l - [sffri,^, fc) - BPfrupz, k)}\. 

where x = 2k / y/s and the other notation is the same as in /3 case. Integration over 
photon angles leads to 



(192) 



df\v) = J dv + j dv^j dxe 5 ^ f(l, v+ jf(l, v _ 



6\v — 1 + (1 — x — u+)(l — v-) \ — 6\v — x — v + — vS\ 761(2;) 
+ 6(v-l + (l-x- v+){l - u_)) 7 2 ^(a;) J, 

fei(x) = -1 + \x, b 2 (x) = -l + l x - i[l + (1 - x) 2 ] ln(1 ~ X) 

Let us show very briefly the calculation of the part proportional to the difference 
of #'s which is somewhat more tricky. We convolute first b\ with photons in the same 
hemisphere and next with photons from opposite hemisphere 



df A (v) = I dVdv. 



J dVdv_ s(v-l + (l- V)(l - -6(v-V- u_) 
A j dxdv + s(y - x-v^f^vj/jjb^x) 



e SYFS f(^,v. 



e &YFS p2( 1 ) 7U 7-1 



dj/2/» 7 (l-y)a 7 ^(l-vj/) ^ 7 -l 







2 V 2 



~ 7 ln(l-i;) ) +0{^ 



(193) 
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The remaining part of d± is easier to calculate because it is explicitly of 0(Y 



dfl{v) = e dYFS F(j) 7^7 





-1 - 




{Si 


- \) 









i(l + (l-) 2 )^^ 



(9( 7 3 



(194) 



The contribution from the initial-state (3\ with 0(a 2 ) prag QED matrix element and with 



0(a 2 



I prag 



analytical integration over the multiphoton phase-space reads 



pf(v) = e s ^s F(y) jv^l \ (-1 + \ 



+ 7 



(195) 



0( 7 3 ). 



The contribution with 0{a) pTag QED matrix element and with analytical 0(a 2 ) prag mul- 

(2) 

tiphoton phase-space integration is obtained by retaining only d\^ and it reads 



+ 7 



1 2 1 
■-V 

2 2 



■1 + l -v\ ln(l 



(196) 



0( 7 3 ). 



5.2.3 Beta-bar-two, (3 2 

In the following step our aim is to calculate analytically the contribution to the total cross 
section from f3\, as given by 



^2 



oo 

E 



1 fd%d% ^«is ;(piiP2 . %)(1 _ e(S j ;;fcj))e y, n ,. 



Ill 



1? <& 



E 

n>j>k>l 



<5 (4) { pi + p 2 - qi - q 2 ~ ^2 k i 

3=1 



(197) 



This contribution is in a sense more trivial than the previous two because it is pure 0(a 2 ), 
it does not have any infrared singularity in the two-photon phase-space integral. 

We can calculate the contribution from f3 2 with the same methods as in the case of (3q 
or f3\. The integral is reorganized easily such that the integration over photon momenta in 

— (2) 

the P2JJ is isolated and we are able to integrate over final-state fermion momenta bringing 
the integral to the standard form 



i 



/ dv p 2 2) (v)a BoTn (s(l-v)). 
o 



(198) 
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Figure 12: The comparison between /C/CMC and IEX 0(a 2 ) prag formulas for the integrated cross 
section as a function of the cut-off parameter u max on ISR and FSR photons. Presented are /C/CMC 
(solid line) and IEX (line with open circles) results for: (a) ISR 0® and the FSR pj®, (b) FSR 0® 

and the ISR (3$ (c) ISR p[f and the FSR 0$. The difference between /C/CMC and IEX is also 
included (dots). Center of the mass energy is 189GeV. Final state fermion is muon. 



(2) 

The function p 2 (v) can be calculated in the LL approximation with either of the 
two methods (keeping an additional spectator photon or not) and after integration over 
photon angles the integral boils down to the following integral over longitudinal photon 
momenta, separately for the case with two collinear and two anticollinear photons 



P?\v) 



dv-dv + — 5{v — v + — vJ) 



+ 



2w_i v. 



-x(w-) X 



2v+v 



I -v. 



-X(v+) x 
1 



— uj(v-) w(v+) 



+ J dv-dv+^- 8(v - 1 + (1 - t>+)(l - »_)) 



(199) 



—^—X(v + )x{v~) - — v(vJ) - — w(v+) — 



2 1 

=7 Z v, 



where x( x ) = (1 + (1 ~ x ) 2 )/2 and u(x) = — 1 + x/2. 

By eventually keeping additional soft photons in the calculation we get our final result 
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Table 2: Contributions to the function pi(v) = d s + A#(w) from fa,k = 0, 1,2. The ISR 
matrix element is at 0(a r ) prag with YFS/EEX exponentiation; r = 0,1,2 is marked in first 
column. Phase-space integration is done analytically always within 0(a 2 ) prag , except of the 
O(a ) prag case in the first row, where the phase-space integration is done in 0(a 3 ) prag . 



for the initial-state 0(a ) prag contribution from /3 2 in a more elegant form 

- Jyfs ey„a ™,7-i J I w „2 l , m/W 3 > 



p^(t,) = e dYFS F(j) -fiH' 1 ^-jv 2 1 + C( 7 d ) (200) 

We have compared numerically the above formula with the /C/CMC in the case of 
FSR switched off and found an agreement better than 0.1%. In Figure |12|(a) we present 
the comparison in which, as in the case of the previous fas, FSR is switched on. In 
Figure |T2](a) we compare the convolution of the ISR p^j an d the FSR (3^ : 

/Umax rv/(l—U max ) 
dv J du 0^(8(1 -u)(l-v))pf^(v)pf^(u) (201) 

The above IEX result is compared with the /C/CMC results, and they agree within 0.2%. 
In Figure |l2|(b) we show the analogous comparison for the convolution of the FSR ft^p and 
the ISR (anticipating IEX results for FSR f3^p in the next section) and we find the 
similar agreement. Finally, there is another more trivial contribution in the fa 2 ^ family 
which correspond to the case with one real photon emitted in the initial state and one in 
the final state. This case does not require a separate analytical phase space integration 
effort because the relevant IEX formula involves the convoluton of the already known 
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expression for the ISR f3{/ and the FSR (3\p. The corresponding numerical comparison 
of the IEX and EEX MC is shown in Figure Figure |12|(c). In fact the IEX matrix element 
was deliberatery constructed in such a way (factorizing virtual corrections) that the above 
convolution-type IEX formula results. 



5.2.4 Summary on IEX for ISR 

The entire initial-state 0(a 2 ) prag integrated cross section is obtained by combining con- 
tributions from all three /5's and it reads 

i 

r (2) _ f J„, s Born/ 



a y> = / dv pf{v)a iiom {s{l-v)), 







-1 + 5 ] 



2 8 

v 1 + 3(1 - vf_ 
2 8 



ln(l - v) 




This above ISR formula has been obtained as a result of ad-hoc exponentiation (interpo- 
lation) in ref. |7j and was used there as a numerical parametrization/testing of the cross 
section from the Monte Carlo program YFS2. It is is now derived starting from YFS 
exclusive exponentiation by means of direct phase-space integration^! 

Summarizing our IEX calculations for ISR, we have obtained through the analytical 
integration over the ISR multiphoton phase space the inclusive exponentiated cross section 
for the IEX matrix elements in the O(a ) prag , C , (a 1 ) pra9 and 0(a 2 ) prag for each i = 
0,1,2 separately. The phase-space integration was always done analytically within the 
0(a 2 ) prag - All results from the above extensive study are summarized in Table ^| where 
we list the two functions ds and A H (v) in the formula 



<7l 



dv Pl (v)a Born (s(l - v)), pi(v) = e &YFS F( 7 ) 7^ (d s + A* («)) 







(203) 



7 ol ( 1 tt 2 \ n a/. s A , e c ^ 

Syfs = i + - -n+T L 7 = 2- In— -1, F(j) 



4 7T V 2 3 J ' ' tx \ m 2 e /' y " r(l + 7 )' 
All notation is recalled for the convenience of the reader. 



5.3 Semi-analytical formulas for FSR 

The calculation of the 0(a 2 ) prag IEX formula for the FSR, with the u max cutoff, that is 
u — 1 — s'/s < n max is quite similar to one in the ISR case and we do not enter into 

22 Ad- hoc exponentiation is of course easier to do and in ref. [Q even the 0(a 3 ) pra g formula for the 
initial-state brcmsstrahlung was given but the derivation method presented here is much better founded 
and the result does not depend on any kind of interpolation or guesswork. 
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Table 3: Contributions to function p F {u) = d' s + A' H (u) from fa, k = 0, 1, 2. The FSR matrix 
element is C(a r ) pras with YFS/EEX exponentiation, r = 0, 1,2 is marked in the first column. 
Phase-space integration is done analytically always in 0(a 2 ) prag . 



details. We only discuss the basic differences between the ISR and FSR cases and present 
the final result. 

If we switch off the ISR completely then the FSR integrated cross section for the 
0(a r ) prag r = 0, 1, 2 EEX matrix element reads 

0>(Mmax) = ^Born J du p F {u), p F {u) = e &YFS F(nff) (d' s + A' H (tt)) 

(204) 

*™ = I - ~ 2 vHi -«) + - (-5 + ¥ j , 7/ = 2- (t^ - ij , 

where the functions <f 5 and A^(-u) obtained with analytical integration of the phase space 
using the 0(a 2 ) prag approximation, are listed in Table [| 

The main difference and a complication in the phase space analytical integration with 
respect to the case of ISR is that the YFS formfactor S' YFS depends in the case of FSR on 
the integration variable u. This is why terms of 0(L 2 a 2 ) are different in the two cases. 
In Table [I] we show separately the contributions from each j3. Note that in the case of 
FSR we did not integrate analytically the phase space for (3 at the 0(a 3 ) prag , like in the 
case of ISR. (It was not necessary in order to reach the precision level of 0.2%.) We have 
checked numerically the agreement of the /C/CMC with the eq. ( |204| ) separately for each 
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Table 4: The comparison between /C/CMC and IEX 0(a 2 ) pra g formula. The difference between 
/C/CMC in EEX mode and the IEX formula divided by Born is plotted with the dotted line, as a function 
of the cutoff u max on the total energy ISR and FSR photons. Include is also the 10~ 2 xcr(-y max )/(TBorn, 
as dotted line for IEX and solid line for MC. 

type of the /3, with the ISR switched off (plots are not shown). We have already presented 
the complete set of numerical results in the case of the ISR switched off in this section, 
for each combination of the ISR and FSR /3's. 

5.4 Semi-analytical IEX for ISR and FSR 

The last numerical test which we show in Figure ^ is the case in which we switch on all 
ISR and FSR /3's listed in both tables | and |. 



It is done for the constant Born cross section, the case with the variable cross section 
will be shown in the next section. We use the IEX formula of the pure 0(a 2 ) prag type 
(without 0(a 3 ) prag improvements for ISR). The everall agreements between IEX and 
/C/CMC is within the advertised 0.2%. By looking into all previous figures in this and the 
previous subsection it is interesting to note that this difference does not come from one 
particular combination of the ISR and FSR /5's, but from several ones. 

The reader may wonder why we elaborate so much in this section for the IEX semi- 
analytical formula which are related to the EEX type of the matrix element in /C/CMC 
if in fact the main matrix element in /C/CMC is now CEEX. The main reason is that 




u/(l— it, 



(205) 



o 







NO 



historically the EEX was the first available example of the exclusive exponentiation, and 
the IEX semianalytical formula were developed in parallel, providing the valuable cross- 
check of the MC. At this stage, as we see in the next section, both IEX and EEX provide 
the reference calculation and valuable test for the CEEX. The precision of the present 
0(a 2 ) prag IEX is limited, but it could be improved to the full 0(a 3 ) prag if necessary. More 
important limitation in the present 0(a 2 ) prag IEX as a test of the CEEX model is the lack 
of the ISR®FSR interference. We believe that this effect can be included in the semian- 
alytical IEX if necessary. The ad-hoc variant of the 0(a l ) exponentiation including the 
ISR®FSR interference is already available in refs. |f22|,p3|1. 
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6 Numerical results and tests 



In this section we shall mainly present the numerical results from the /C/CMC in which 
the Standard Model amplitudes for the process e~e + — > f f + n r y of the previous Section |2| 
(EEX) and Section ^ (CEEX) are implemented. The analytical results of the Section |5| 
will be also exploited to obtain numerical results from the semianalytical program /C/Csem. 
These results are mainly for the fi~fi + final state. For more results on the quark final states 
and other interesting numerical results from /C/CMC we refer the reader to forthcoming 
proceedings of the LEP2 Monte Carlo Workshop p5 |. 
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Figure 13: General structure of the /C/CMonte Carlo program. 



The general structure of the /C/CMC code is depicted in the Figure [L3|. The program 
is divided in the two distinct parts (levels): 

(a) Phase-space Monte-Carlo integration engine with the common importance sampling 
for the entire family of QED distributions (EEX and CEEX) 

(a) Collection (library) of programs for the SM/QED spin amplitudes and differential 
distributions, at various orders, with various styles of exponentiation. 

In this work we do not enter into the description of the MC integration algorithm in the 
universal MC integration engine. The Monte Carlo method of the phase space integration 
is fully documented (for the first time) in ref. Q and some aspects of the phase-space 



82 



parametrization are documented in the forthcoming ref. Here we regard this low- 

level MC program as a black box capable to integrate the phase space exactly (up to a 
statistical error). 

The life is however not that simple and the numerical program which "in principle" 
is doing something "exactly /rigorously" may still give imprecise results due to program- 
ming bugs and numerical instabilities, especially in a program as complicated as /C/CMC 
is. This is why we always introduce a concept of the technical precision of the given pro- 
gram/calculation, see below. The basic aim of our numerical exercises presented in this 
Section is the determination of the total theoretical precision associated with our calcula- 
tion of Standard Model predictions for experimental observables (although we limit the 
discussion to QED part of SM for most of our discussion). As for observables we shall 
concentrate mainly on the total cross section and charge asymmetry at LEP1, LEP2 and 
Linear Collider energies. 

What are the technical and physical precisions? The technical precision we define as 
all uncertainties related to pure numerical problems like programming bugs, numerical 
instabilities, numerical approximation, etc. In our case the question of the technical 
precision will mainly concern the MC integration engine. It is important to determine it 
at the early stage of the work and it should be generally much better than the physical 
precision. On the other hand, the physical precision is the total uncertainty related to 
neglected higher orders in coupling constant a or in other expansion parameters like the 
inverse of the big-log 1/L, or the ratio of the width to mass T/M for a narrow resonance. 
For physical precision we understand that the above truncations are done in the spin 
amplitudes and/or differential cross section. If some of them are done in the phase space 
integration then we tend to associate them with the technical precision (as phase space 
integration is a technical problem). 

We shall start this section with the basic discussion of the technical precision, then 
we proceed to a subsection elaborating on the physical precision for the EEX matrix 
element, based on comparisons of the /C/CMC and semi- analytical results, and later we 
discuss the physical precision for the case of the full CEEX matrix element. In this Section 
we also present numerical results and a rather complete discussion of the effects due to 
the ISR-FSR interference in the fermion pair production process. 

We note that it would be good to include also more numerical tests and lower energies 
~ lOGeV, and for very high energies ~ ITeV, and some more tests specific to spin 
effects. However, the basic pattern of the spin correlations in double r decay was already 
cross-checked in ref. [Q]. 

6.1 Basic test of technical precision 

The best way to determine technical precision is to compare results of the two or even 
more independent calculations which implement the same physics model but differ in 
technical details of the actual implementation like the method of phase space integration, 
independent coding, etc. The best two possible methods are: (a) to compare two indepen- 
dent Monte Carlo calculations or (b) to compare Monte Carlo results with the results of a 
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Figure 14: Evaluation of the technical precision of the /C/C Monte Carlo using simplified QED 
multiphoton distribution. The difference of the /C/C MC result and semianalytical result divided by 
semianalytical is plotted as a function of v max = 1 — s' min /s. results are shown for p + p~ final state 
at y 7 ! = 189GeV. In the case (a) the phase-space limit v mSLX = 1 — 4m^/s is taken; the last bin 
represent the entire phase phase space. In the case (b) v max = 0.999. 



semi-analytical calculation. The method (a) is generally better because it can be done for 
arbitrary kinematical selections (cuts) and for the simplified QED matrix element, while 
method (b) is limited to a simple or absent kinematical selections. In the following we 
shall use method (b). 

For our basic test of the technical precision we use the simplest possible variant of 
the QED model, that is of the type C(a°)EEX defined in Section §. For this type of 
QED matrix element we were able to integrate analytically the total cross section in the 
Section |5|. The relevant formula can be read from the first row in Tables ^| and [3]. For the 
sake of completeness we write down the complete expression explicitly: 

/■1'max 

4vn = / dv a f Born (s(l - u)(l - «)) pf\v) P P(u), 
Jo 

pf\v) = F( 7e ) e^(^) le v^ x [\ - \ le ln(l -v)-~ hi 2 (l -v) + 7e 2 ) , 
p P(u) = F{ lf ) e h/-h/Mi-«)+f(*+£) 7/U 7/-i _ ^ 7/ i n (i _ u ^ 

(206) 

As we remember the coefficient in front of the 0(L 3 a 3 ) term is zero, as marked explicitly. 
It was essential to calculate analytically and introduce the ISR term of ^(L 1 ^ 2 ) because 
it amounts numerically to several per cent for the cross section located close to v — 1. 
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In fig. [T4| we present the comparison of the /C/C MC with the semianalytical formula of 
eq. ( P06|) . The difference between the MC result and the semianalytical result is divided 
by the semianalytical result and as we see it is remarkably small! The comparison is done 
for the final state at y/s = 189 GeV, as a function of v mSuX . In the last point (bin) 

the entire phase space is covered, w max = 1 — 4m^/ s. 

The conclusion from the above exercise is that we control the phase-space integration 
at the level of 2 x 10~ 4 for t> max < 0.999, including the Z radiative return, and at the level 
of 3 x 10~ 3 for no cuts at all. 

The possible loophole in this estimate of precision is that it may break down when we 
cut the transverse momenta of the real photons, or switch to a more sophisticated QED 
model. The second is very unlikely as the phase space and the actual SM model matrix 
element are separated into completely separate modules in the program. The question of 
the cut transverse momenta of the real photons requires further discussion. Here, it has 
to be stressed that in our MC the so-called big-logarithm 

L = In (^) - 1 (207) 

is the result of the phase space integration and if this integration were not correct then 
we would witness the breakdown of the infrared (IR) cancellation and the fermion mass 
cancellation for FSR. We do not see anything like that at the 0.02% precision level. In 
addition there is a wealth of comparisons with many independent codes of the phase 
space integration for n 7 = 1,2,3 real photons, with and without cuts on photon px- 
It should be remembered that the multiphoton phase space integration module/code in 
/C/CMC is unchanged since last 10 years. For ISR it is based on YFS2 algorithm of 
ref. [[F] and for FSR on YFS3 algorithm of ref. 0, these modules/codes were part of 
the KORALZ Jl0| multiphoton MC from the very beginning, already at the time of the 
LEP1 1989 workshop f40fl , and they were continuously tested since then. The phase space 
integration for n 7 = 1 was tested very early by the authors of YFS2/YFS3 against the 



older MC programs MUSTRAAL [|16j and KORALB JJIJ and with analytical calculations, 
at the precision level < 0.1%, with and without cuts on photon px- The phase space 
integration for n 1 = 2, 3 with cuts on photon px was tested very many times over the 
years by the authors of the YFS2/YFS3/KORALZ and independently by all four LEP 
collaborations, using other integration programs like COMPHEP, GRACE and other ones, 
in the context of the search of the anomalous 27 and 37 events. Another important series 
of tests was done in ref. |T7| for ISR n 7 = 1,2 photons (with cuts sensitive to pt of 
photons), comparing KORALZ/YFS2 with the other independent MC's for the 1^7(7) 
final states. Typically, these tests, in which QED matrix element was programmed in 
several independent ways, showed agreement at the level of 10% for the cross section for 
n 7 = 2 which was of order 0.1% of the Born, or 0.2-0.5% for n 7 = 1 which was of order 1% 
of the Born, so they never invalidated our present technical precision of 0.02% in terms 
of Born cross section (or total cross section in terms of Z-inclusive cut). 

We conclude therefore that the technical precision of /C/CMC due to phase space inte- 
gration is 0.02% of the integrated cross section, for any cuts on photon energies Z-inclusive 
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and Z-exclusive, stronger thanf^] M inv (ff) > 0.1-y/i and any mild cut on the transverse 
photon energies due to any typical realistic experimental cuts. For the cross sections with 
a single photon tagged it is about 0.2-0.5% and with two photons tagged it is ~ 10% of the 
corresponding integrated cross section. These conclusions are based on the comparisons 
with at least six other independent codes. 

6.2 Physical precision, the case of EEX 
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We now start the presentation of the numerical results from /C/CMC run in the EEX 
mode with semianalytical calculations based of the results in Section |^. Note that the 
EEX matrix element of Section ^| is very similar (basically the same) to that implemented 
since many years in KORALZ program ||10|| . We do for two reasons: (a) these tests were 
historically the first, (they existed unpublished for many years giving us confidence that 
the KORALZ /YFS3 program provides correct results) and (b) they are now still useful as 
a reference calculation for the newer CEEX scheme. They will also allow us to introduce 
some notations and gradually introduce the reader to the subject of the discussion of the 
theoretical precision of our results. Of course, we shall remember that in the case of EEX 
we do not include the ISR-FSR interferences (IFI). 

In Figure [15] we show the dependence of the total cross section on the cut on the total 
photon energy (ISR+FSR). The comparison is done for the fi + fi~ final state at y's =189 
GeV, as a function of v max - In the last point (bin) the entire phase space is covered, i.e. 

23 It downgrades to 0.5% for M inv (fj,p) < 2m^, i.e. full phase space. 
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fmax = 1 — 4m^/s. The very striking (and well known) phenomenon is that the total 
cross section due to huge ISR correction is almost 3 times the Born cross section, in the 
absence of any kinematical cuts. Part of this ISR contribution is located close to v = 1, 
s' ~ 4m^/s, let us call it the 77* process, it amounts to as much as the Born cross section 
itself, a 77 . ~ cr Bo rn, while dominant part of the cross section ctzrr ~ 2o"Born is concentrated 
close to v = 1 — Mg/s ~ 0.75, and is associated with the so called U Z radiative return" 
(ZRR) process, that is the resonant production of Z, after emission of rather hard ISR 
photon, usually lost in the beam pipe. In the experiment the 77* process is almost always 
eliminated from the data, and the ZRR process is also not very often included in the data 
sample. The typical experimental cut is situated somewhere in the range 0.1 < v max < 0.3. 
As we see in Figure [15] (a), the total QED corrections (cr(v max ) — o"Born)/ cr Born is in this 
case quite close to zero, in fact slightly negative. 

In Figure |T5|(b) we compare the AT/CMC with the semianalytical expression based on 
the phase-space integration in Section |5|. In the MC calculation we use the second order 
EEX type of the QED model EEX2=C(1, a, aL, a 2 L 2 ) EEX , defined in Section @. The 
semianalytical formula used in Figure 0(b) is also in the class EEX2. It is defined as 
follows 

/•Umax 

"san = / dv "L>(1 - «)(1 - «)) P?>) (208) 
J 

where the distributions pp and py are from the Tables |2| and |3] in Section |5|. 

What kind of lesson can we draw from Figure |T5|(b)? We treat the result in Fig- 
ure fL"5Kb) as an indication that, the contribution from QED (non-IFI) photonic corrections 
to combined physical and technical precision in the EEX2-class integrated cross section 
for the standard cut v max ~ 0.2 is about 0.2%, for the ZRR process it is 0.7% and for 
the 77* process it is 3%. We are here talking about the technical precision of the coding 
of the EEX2 matrix element, not associated with the phase space integration (covered in 
the previous section). 

In the next Figure [16| we make an attempt to estimate the physical precision of the 
QED model in the EEX class. Specifically, we look into difference between EEX2 (as de- 
fined above) and EEX1, with the EEX1 being the C^c^eex of Section g EEX1=C(1, a, aL) 
It is plotted in Figure |l6l(a) both from /C/CMC and semianalytical formula. Taking con- 
servatively (see the discussion below) half of the difference between EEX2 and EEX1 as 
an estimate of the physical precision of EEX2 we arrive to a similar estimate of about 
0.2% for the standard cut v max ~ 0.2, 0.7% for the ZRR process and up to 3% for the 77* 
process. 

The other useful piece of information comes from Figure ^6|(b), where we plot the 
difference EEX3-EEX2, with EEX3=C?(1, a, aL, a 2 L 2 , a 3 L 3 ) EEX , provides direct insight 
into the neglected third order LL contributions. As we see it is always below 3 • 10 -4 , 
(This estimate will be also useful for the case of CEEX.) If the 0{L 3 a 3 ) corrections is of 
this size, then necessarily the main contribution to the above estimate of the theoretical 
error is coming from the 0(L l a 2 ) corrections! 

In fact the lack of the 0(a 2 L 1 ) corrections in both EEX2 and EEX1 is the main 
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Figure 16: An attempt to estimate physical precision for EEX: 0(a 2 ) and 0{a 3 ) Process, energy, 
definition of cuts the same as in Fig. [L5| . 



deficiency of the above tests, so they cannot pin down directly the size of this contribution. 
Keeping this limitation in mind, from the above test we nevertheless estimate tentatively 
the combined physical and technical precision in the integrated EEX3-class cross section 
of the /C/CMC to be 0.2% for the standard cut t> max ~ 0.2, 0.7% for the ZRR process 
and about to 1.5% for the 77* process. The caveat of this exercise is that we know 
retrospectively the QED non-IFI component of the precision on the KORALZ/YFS3 
Monte Carlo at LEP2 energies because the EEX of KORALZ and the EEX of /C/CMC 
are practically the sarneQ. The above does help indeed, in spite of the fact that the 
neglected IFI contribution to integrated cross section is of order 1%, because KORALZ 
in the non-exponentiated 0(a) mode can calculate IFI separately, see discussion in the 
following subsections. 

Let us finally make an ultimate effort to estimate the total precision, staying all 
the time within the EEX model. As we have already noted the most important miss- 
ing contribution seems to be the 0(L l a 2 ), most probably the ISR part of it. In the 
the semianalytical formula for the total cross section we are able to add it, since it is 
known from ref. [18]. The (D(L 3 a 3 ) corrections we may add as well and in this way 
we replace py by the p^f 1 of ref. which is the true 0(a 3 ) prag for ISR (according 
to the terminology in the Introduction) and 0(a 2 ) prag for FSR (no IFI). Let us call 
it EEX3best=C(l, a, aL, a 2 L 2 , a 2 L l , ck 3 L 3 )eex- The difference between semianalytical 



EEX3best and EEX3 from /C/CMC is plotted in Figure [17[ As we see this final test 
confirms the previous estimate of the physical precision of the EEX type of the matrix 
element. 



4 KORALZ version 4.02 and earlier have EEX implemented differently from /C/CMC. 
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6.3 Physical precision, the case of CEEX 

Quantitative determination of the physical precision should be based on the comparison 
of the calculations in two consecutive orders in the expansion parameters, for instance by 
comparing results from 0(a r ) and 0(ct r ~ l ) calculation, or 0(L r a n ) versus 0{L r ~ l a n ), 
etc. For example when only the Born and 0{a 1 ) results are available one should take 
the difference between the two (or some fraction of it) as an estimate of the physical 
precision. The above conservative recipe gives solid estimate of the physical precision and 
we shall employ it as our basic method in the following. In most cases in the literature, 
however, authors try to estimate the uncalculated higher order effects with some "rule of 
thumb". For instance in the case when Born and 0(a 1 ) results are known they take \L- 
as an estimate of the missing/uncalculated 0(a> 2 ) corrections. This has to be done with 
care because one may easily overlook some "enhancement factor" . For example the cross 
section close to a resonance can be modified by additional powers of the big logarithm 
In jr. In most cases these "enhancement factors" are already seen in the 0(a 1 ) calculation 
so it is not difficult to trace them. 

We are in rather comfortable situation because for QED "photonic" corrections we 
have at our disposal O(a ), 0{a 1 ) and 0(a 2 ) calculations (at least for ISR, where they 
are the biggest). We can therefore afford to take half of the difference between the 
0{a 1 ) and 0(a 2 ) calculations as a conservative estimate of the physical precision due to 
QED "photonic" corrections. We also profit from the fact that the exponentiation speeds 
up considerably the convergence of the perturbative series by "advanced summation" of 
certain class of corrections to infinite order, and by not introducing additional spurious 
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cut-off parameters dividing real emissions into soft and hard ones which are typical for 
the calculations without exponentiation (see discussion on the famous ko parameter in 
the 1989 LEP workshop gj). 

Let us mention that in our estimates of the physical precision we omit from the dis- 
cussion the 0(a 2 ) effects due to an additional fermion pair, either real or virtual. We 
do it because: (a) there are many MC programs which implement production of the four 
fermion final states (often with additional ISR) and (b) in the experiment this contribu- 
tion can be eliminated at the early stage from the data in the experimental data analysis 
aimed at single fermion pair production, see for example ref. ||48|| . In fact this point is 



still under debate, see forthcoming proceedings of the LEP2 Monte Carlo workshop 
It was proposed that in the final combined LEP2 data certain the so called non-singlet 
initial and final state secondary pair contribution will be kept in the data, as done by 
OPAL, see refs. ||49|-pT|. We have recently included the virtual corrections of the "vac- 
uum polarization" type with the fermionic bubble in the 0(a 2 ) photonic contributions 
like vertex corrections in yet unpublished version of 4.14 of /C/CMC. This is done having 
in mind combining results of /C/CMC with the other MC program for four-fermion pro- 
duction process like KORALW [52 ] . The tandem of /C/CMC and KORALW programs will 
be able to realize any possible scenario of the treatment of the soft/light pair corrections 
in the LEP2 data. 

In Fig. [18| we present the numerical results on which we base our quantitative esti- 
mate of the physical precision due to photonic QED corrections. In this figure we plot 
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the difference between 0(a 2 )cEEX and C(a 1 )cEEX for the total cross section and charge 
asymmetry at 189GeV as a function of the cut on the total energy emitted by all ISR 
and FSR photons for final state. The cut is formulated with the s' > s' min or equiv- 

alently v < w max condition, where s' is the effective mass squared of the pair and 

v = 1 — s'/s, as usual. One should remember that the actual experimental cut is around 
% M ~ 0.2 (eliminating Z radiative return) in the case of the standard data analysis, and 
sometimes around i> max ~ 0.9 in the case when Z radiative return is admitted in the data. 
The "kink" around v max ~ 0.75 is at the position of the Z radiative return. In either 
case, whether we admit or eliminate the Z radiative return, that is for f max ~ 0.9, the 
difference between C(q! 2 )ceex and 0(a 1 )cEEX m total cross section is below 0.4% and for 
the charge asymmetry it is below 0.002. 
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Figure 19: Evaluation of the physical precision for total cross section and charge asymmetry, at 
■y/s = 500GeV. The difference between 0(« 2 )ceex and 0(o 1 )ceex is plotted as a function of 



s' min /s for final state. 



Taking conservatively half of this difference among 0(« 2 )ceex and 0(« 1 )ceex as an 
estimate of the neglected 0(a 3 )cEEX and higher orders we conclude that the physical 
precision due to photonic QED corrections of our C(q! 2 )ceex calculation for all possible 
cutoffs within < u max < 0.9 range is 0.2% in the total cross reaction and 0.001 in the 
charge asymmetry. This estimate would be even a factor of two better, if we restricted 
ourselves to the most typical cut-off range 0.1 < t> max < 0.3. The above estimate will be 
confirmed by more auxiliary tests in the following. 

As we see we have improved on the physical precision estimate with respect to the 
previous estimates for the EEX model - in addition we do include IFI all the time. 
For the ZRR process we now quote for the integrated cross-section 0.2% instead of the 
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previous 0.7% and for 77* we have something like 0.3% instead of the previous 1.5%. 
This we interpret as a result of inclusion of the 0(La 2 ) ISR correction in our CEEX spin 
amplitudes. 

We have to stress very strongly that the estimate of the physical precision depends on 
the type of the observable (we took a and A FB ), the type of the final state (we took \i pair 
final state; for the quark final state the QED FSR effects are smaller due to the smaller 
electric charges of quarks) and on many other input parameters, for example the total 
CMS energy. The great thing about the Monte Carlo is that the type of the evaluation 
we proposed and implemented in this Section (half of difference 0(a 2 )— 0(a 1 )) can be 
repeated for any observable, any final state and any energy. For example in Figure [19 
we repeat our evaluation of the physical precision for a and A F b at the Linear Collider 
energy 500GeV. As we see the resulting precision is worse, the worsening is negligible for 
a mild cut of order v max < 0.5 and almost factor two for the Z radiative return, which is 
now placed close the v = 0.95. 



6.4 Absolute predictions, more on physical/technical precision 

In this Section we shall present the SM absolute predictions for the total cross section and 
charge asymmetry at LEP2 (189GeV) and at the Linear Collider (500GeV). We compare 
them with our own semianalytical program /C/Csem, with KORALZ [[M]] and in some 
cases with ZFITTER ||. They may not improve our basic estimates of the technical 
and physical precision from the previous sections, but they can confirm them (or disprove 
them!). 

In Table || we show numerical results for the total cross-section cr(f max ) and charge 
asymmetry A FB (v max ) as a function of the cut v ma _ x on the total photon energy (the 
cut-off parameter v mSuX is defined as in the previous subsection). Generally, in Table |5] 
we show results with the ISR-FSR interference (IFI) switched on and off. The /C/Csem 
semianalytical program (part of the /C/CMC package) provides "reference results" for cr 
and Afb, see the first column in Table [|, which are without IFI, obtained from using 
EEX3best formula defined in the previous Section |6.2[ For the charge asymmetry we use 
the convolution- type semianalytical formula like that of eq. (|208Q . (In fact we use this 
formula separately for the cross section in the forward and backward hemispheres and 
then we calculate A f -q from these partial integrals.) Results from the /C/CMC in Table |5] 
are shown for two types of QED matrix element: 0(« 2 )ceex with and without IFI. In 
addition we include results from KORALZ are for the 0(a l ) matrix element with and 
without IFI which will be discussed in the next Section. 

As tables with list of numbers are difficult to comprehend, we present the essential 
results of the Table |5| in Figure |20|, where they are all plotted as a difference with the 
reference results of our semianalytical program /C/Csem. (In other words the results from 
/C/Csem are exactly on the x-axis.) 

In the case of IFI switched on /C/Csem cannot be used as a cross-check of the /C/CMC. 
Remembering that IFI in KORALZ in the 0{a 1 ) mode (without exponentiation) is very 
well tested we combine the 0(a l ) IFI contribution with the CEEX result without IFI. Such 
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.5852 ± .0008 


.5947 ±.0009 


.5972 ± .0004 


.0078 ± 


.0004 


.90 


.3105 ± .0000 


.3115 ±.0004 


.3096 ± .0004 


.3170 ± .0005 


.3260 ± .0002 


.0037 ± 


.0002 


.99 


.2851 ± .0000 


.2867 ±.0004 


.2843 ± .0004 


.2912 ± .0004 


.3039 ± .0002 


.0024 ± 


.0002 



Table 5: Absolute prediction for total cross section an charge asymmetry. They are for final 
state at yfs = 189GeV. Results are plotted as a function of the cut-off on the total photon energy 
Vmax = 1 — s min/ s - The "reference" a and Afb in first column are from /C/Csem semi-analytical 
program. We have used Higgs boson mass lOOGeV and top mass 175GeV as input parameters. 



a hybrid solution denoted in Fig. ^0] as "CEEX2+IFI at (9(a )" us used as our primary 



test of the full CEEX matrix element with IFI switched on. The above procedure is 
done separately for cross sections in the forward and backward hemispheres such that the 
prediction for charge asymetry is also available. 

It is worth to mention that the above hybrid solution was already successfully used in 
refs. [[)3|,[)4]] for the study of the IFI contribution at Z peak, imposing strong acollinearity 
cut. It is also implemented in a semianalytical form in ZFITTER 6.x. On general ground 
we expect this recipe to be rather good, because IFI correction itself at 0(a l ) does not 
contain any large mass logarithm and is relatively small and can be handled additively. 

In Figure 20 we also show the numerical results from /C/CMC in the EEX3 mode (no 
IFI), from KORALZ in the EEX2 mode (no IFI), which are not included in Table 

Let us now comment on the results in Figure The EEX3 from /C/CMC differs from 
EEX3best of /C/Csem (no IFI in both) by about 0.7% for the ZRR process, as we have 
already seen, and we interpret this difference as the result of the missing Oil^o?^. The 
EEX2 of KORALZ 4.03 is closer to the EEX3best of /C/Csem for ZRR process - we do not 
see any contradiction in this since the implementation of EEX in KORALZ and /C/CMC 
differs in the details (causing a difference of 0{L l a 2 ) in the integrated cross section.) 

In the case of IFI switched off, the CEEX2 result, corresponding exactly to C(q; 2 )ceex, 
defined in Section [|, as implemented in /C/CMC 4.13, agrees very well with the EEX3best 
of /C/Csem. This result is compatible with the total theoretical precision of 0.2% for the 
integrated cross-section, even including the ZRR process. 
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Figure 20: Absolute prediction for total cross section an charge asymmetry. They are for final 
state at y/s = 189GeV. Results are plotted as a function of the cut-off on the total photon energy 
Vmax = 1 — s min/ s - The "reference" a and Afb in the first column are from /C/Csem semi-analytical 
program. 



In the case of IFI switched on, the hybrid solution "CEEX2+IFI at (^(a 1 )" also 
agrees with the full CEEX2 result confirming the total theoretical precision of 0.2% for 
the integrated cross-section, including the ZRR process. 



For the charge asymmetry in Figure 20 situation is quite similar. The IFI effect is 



up to 4% for strong cuts. In the case of IFI switched off, the CEEX2 result agrees with 
the EEX3best of /C/Csem to within 0.2%. For IFI included, the CEEX2 agrees with the 
hybrid solution rather well, to within 0.4%. Note that in the above Monte Carlo exercise 
we have used the symmetric definition of the scattering angle 9* of ref. [|55| (which is close 
to what is used in the LEP experiments). 

Summarising, the numerical results in Figure |2l] establish our basic estimate of the 
theoretical precision of the /C/CMC, due to QED effects, at LEP2 energies of about 0.2% 
for total cross section and 0.2-0.4% (depending on cut-offs) for charge asymmetry. 



Finally, we examine the analogous results from the /C/CMC at 500GeV in Figure |21|. In 
this case we include only results from the /C/CMC and /C/Csem. The pattern of agreement 
is up to a factor two the same as at 189GeV. 

6.5 Initial-final state interference 

The control of the initial-final state interference correction down to the precision of 0.2% 
in the integrated cross-section and in the charge asymmetry is rather important - this 
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Figure 21: Total cross section an charge asymmetry for fi + fi final state at ^fs = 500GeV. Results 



analogous as in Fig. 20 



is why we dedicate this section to a more detailed study of this QED correction. In 
particular we would like to answer the following questions: 

• How big is the ISR®FSR interference in a to t, A-fb^- 

• Do we know ISR®FSR at C^a 1 )? 

• Do we know ISR®FSR beyond 0(^)1 

• How sensitive is ISR®FSR to cut-off changes? 

KORALZ is the best starting point and reference for the problem of calculating the 



ISR®FSR. In Figure p2| we show results from the 0{a 1 ) KORALZ (no exponentiation) for 
the fi + fi~ final state at i/i=189GeV. Angular distributions from KORALZ, pure 0(a l ) 
(without exponentiation), were verified very precisely at the level of ~ 0.01% using a 



special analytical calculation, see ref. 0], so we know ISR®FSR at Oi^a 1 ) very precisely. 
As we see the ISR®FSR contribution to the integrated cross-section is about 3% and 
about 0.03 to A-p^. This is definitely above the ultimate experimental error tag for the 
combined LEP2 data at the end of LEP2 operation. The energy cut on the total photon 
energy is fixed in the results of Figure ^ to just one value v < v ma , x = 0.1 (where 
u max = 1 — s'/s is defined as usual). This is close to the usual value in the experimental 
LEP2 data analysis. We introduce also the angular cut | cos6*| < cos6* max and vary the 
value of cos^max, see Figure ^2](b), the value used in the experimental LEP2 data analysis 
is around cos# max = 0.9; this corresponds to two bins before the last one in Figure p2j(b) 
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Figure 22: Results from 0(a 1 ) KORALZ (no exponentiation) for /i + fi final state at A/s=189GeV. 
The energy cut is on s'/s, where s' = m 2 j- The angular cut is | cos9\ < cos# max . The scattering 



angle is 9 = 9* of ref. [55] 



(the last point in the plot is for cos# max = 1). In this way already we have answered the 
first two questions from the list above. 

In Figure [2^ we present similar results from the OCMC which will help us to answer 
whether we know ISR®FSR beyond 0{a 1 ) and inspect in a more detail the dependence on 
cut-offs. In Figure |23|(a-b) we essentially repeat the exercise of Figure |22| finding out the 
ISRcgESR contribution to the angular distribution and A-fb for the same energy-cut using 
/C/CMC instead of KORALZ. As we see the results change slightly, the ISR®FSR effect is 
about 20%-30% smaller. We attribute it mainly to (a) different (better) treatment of the 
ISR in /C/CMC (b) exponentiation of the ISR®FSR effect is in /C/CMC. As it is well known, 
in 0(a l ), the ISR®FSR contributes like ^Q e Qf^ In \+ToJe to tne an gular distribution - 
this even causes the angular distribution to be negative close to cos# = — 1. In the CEEX 
exponentiation the above singularity is summed up to infinite order and the angular 
distribution near | cos#| = 1 is not singular any more. (This kind of exponentiation will 
be implemented in the next version of ZFITTER, see J|5| for first numerical results.) 
The typical experimental cut | cos#| < 0.9 eliminates most of the above trouble anyway 
- what is probably more important is the correct "convolution" of the IR-finite 0{a 1 ) 
ISR(g>FSR with the 0(a 2 ) ISR. In the /C/CMC this is done in a maximally clean way from 
the theoretical/physical point of view (at the amplitude level) while in the semianalytical 
programs like ZFITTER || this is done in a more "ad hoc" manner. Let us remind the 
reader that we still lack the genuine IR-finite 0(a 2 ) corrections in the ISR®FSR class 
from diagrams like 2-boxes and 5-boxes, see Section |3|. These contributions are most likely 
negligible, of order O^a 2 ) at most. 
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Figure 23: Results from <D(a 2 ) /C7CMC. for fi + fi final state at ^/s=189Ge\/. The energy cut is on 
s'/s, where s' = frv^j- The angular cut is | cos#| < cos0 max . Scattering angle is 6 = 9* of ref. [55|. 



In Figure £3|(c-d) we make the energy cut more loose, t> max = 0.9, thus admitting the 
ZRR into the available phase-space. As a result the relative ISR®FSR decreases by a 
factor 3, simply because it gets "diluted" in the factor 3 bigger integrated cross section, 
while ZRR does not contribute to ISR®FSR because of its narrow-resonance character, 
already discussed in Section |3] at length. The fact that the ZRR does not contribute to the 
ISR®FSR can be seen explicitly in Figure |2l where we plot the ISR®FSR contribution to 
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Figure 24: ISR&FSR contribution to A^b "bin-per-bin" . 



A-p-g "bin-per-bin" , that is calculated in each bin separately. As we see the contribution 
from the ZRR which at this energy (189GeV) is located at v — 0.75 is very small, smaller 
than from all other v's where there is no Z resonance. 
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Figure 25: Results from 0(a 1 ) from KXMC for final state at A/s=189GeV. The energy cut 

is on v p = 1 — s'/ s ' where s' is estimate of the "propagator eff. mass" as defined by ALEPH. The 



angular cut is | cos6>| < cos# max . Scattering angle is 9 = 6' of ref. |55] 
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In the above exercises and also in the following we use always the energy cut on the 
v = 1 — s' / s variable defined in terms of the effective mass of the "bare" final fermions, 
that is without any attempt of combining them with the collinear FSR photons. This is 
experimentally well justified for the /i-pair final states but not for r-pairs or quarks. It is 
possible and in fact rather easy to define a "propagator" or "reduced" s' p which takes into 
account the loss of energy due to ISR but not FSR. In other words the s' p effective mass 
squared sums up FSR photons. One can ask a legitimate question: If we would cut not 
on the "bare" final fermion variable s', but instead on the "propagator" s' p then perhaps 
the estimate of the ISR® FSR contribution would be dramatically different, for instance 



it would be much smaller? In Figure |26] we show a numerical exercise in which we employ 
the energy cut in terms of v p — 1 — s'j s. One can construct such a s' looking into angles 

In Figure ^ we use 



of the outgoing fermions. This type of variable was used in ref |55 



the definition of s' p of ALEPH ||56|| . As we see in this Figure, the result is not dramatically 
different from what we have seen in Figure The magnitude of ISRcgESR contribution 
is close to what we could see if we applied the same value of the energy cut for the "bare" 
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We shall now examine the dependence of the ISR®FSR contribution on the energy- 
cut v max in a more detail. In Figure we show the ISR®FSR contribution to A FB 
as a function of energy-cut v max at two energies (a) 189GeV and (b) at the Z peak 
y/s = Mz- No cut on cos# is applied. In addition to /C/CMC results we show the results 
from 0(a l ) mode of KORALZ and from ZFITTER. At 189GeV and for the typical 
energy cut 0.2 < v max < 0.3 all three programs agree very well. This cut is relatively 
"inclusive" such that exponentiation effects are not so important and ISR is eliminated 
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in a "gentle" way (the total cross section is close to the Born value). For stronger cuts 
u m£Cf < 0.2 we see a large (factor 2) discrepancy among the results from /C/CMC and 
both KORALZ and ZFITTER, because of the lack of exponentiation in KORALZ and 
ZFITTER. (in ZFITTER ISRtgFSR is taken without exponentiation and combined with 
ISR "additively"). We also observe the discrepancy of about 0.2% among /C/CMC on one 
hand and both KORALZ and ZFITTER on the other hand for the ZRR. Our guess is 
that it is due to the difference in the method of combining ISR®FSR with the second 
order ISR (of course, we believe that the CEEX method of doing it at the amplitude level 
is the best one can do). In Figure |26|(b) we see, first of all, the well known phenomenon 
of the strong suppression of the ISR®FSR contribution at the resonance, especially for 
the loose cut-off. Even for a strong cut, v ma , x = 0.1, the ISR®FSR contribution is about 
0.01, about factor 30 smaller than in the off-resonance case. Here, /C/CMC agrees rather 
well with KORALZ and ZFITTER. The differences are generally^ up to 0.0015. 
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Figure 27: s'-cut dependence of 5a, No 9-cut. 



In Figure |27|(a) we examine the ISR®FSR contribution to integrated cross section as a 
function of the energy cut v max . At 189GeV and for the typical energy cut 0.1 < f max < 0.6 
all three programs agree reasonably well, KORALZ and ZFITTER are generally closer to 
each other than to /C/CMC. After admitting the ZRR, t> max > 0.8 all three programs agree 
even better. For a very strong cut, t> max < 0.1 KORALZ and ZFITTER disagree dramat- 
ically with /C/CMC because of the lack of exponentiation in KORALZ and ZFITTER. In 
Figure |27|(b) we see, again the strong suppression of the ISR®FSR contribution at the 
resonance, especially for the loose cut-off. Suppression is cut-off dependent and generally 

25 The difference between KORALZ and ZFITTER should be perhaps smaller because both are 0(a 1 )? 
May be the difference is due to the angle definition? 
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stronger for KORALZ and ZFITTER than for /C/CMC. Most of the comments which we 
made on ISR®FSR contribution to A^b apply also here. 
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Figure 28: Back on Z peak. 



Finally in Figure ^8] we go back to a vicinity of the Z peak (LEP1) and we show the 
magnitude of the ISR®FSR contribution to the integrated cross section as a function of 
the CMS energy, for the \x~ fi + final state and for all five quark final states taken together 
(the so-called hadronic cross section) from the /C/CMonte Carlo. No angular cut or energy 
cut is applied (full phase space). For the final state we also include results from the 

KORALZ 0(a l ) and ZFITTER/TOPAZO Results from quarks are multiplied by a 

factor 10 to be visible, because ISR®FSR contribution in this case is small. It is not only 
suppressed by the smallness of the quark charge, but we also have partial cancellation 
among up- and down- type quarks, see ref p5fl . However, the ISR®FSR contribution 
to hadronic cross section has to be known much more precisely (factor ~ 3) because it 
is measured much more precisely, due to higher statistics. In Fig. ^ we see that the 
suppression of ISR®FSR is much weaker as we go away from the centre of the resonance, 
and it changes the resonance curve in such a way that it affects the fitted mass of the 



Z. The actual size of the shift of Mz was studied in ref. |57j and it was found to be 
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0.15MeV. Results of the /OCMC are smaller about 10-20% than the Ola 1 ) estimates of 
KORALZ 0{a 1 ) and ZFITTER, away from the Z peak. This is compatible with the 10- 
20% size of the 0(L 2 a 2 ) ISR corrections with respect to O^^a 1 ) corrections, which are 
included in /C/CMC and are not included in KORALZ and apparently also not included in 
ZFITTER/TOPAZ0 (which agree very well with KORALZ). Our last comment concerns 
the reliability of our estimate for the ISR®FSR contribution in the absence of the correct 
implementation of the simultaneous emission of the FSR photon and FSR gluon. We 
think that through the usual arguments, see ref. H57|, we can neglect from consideration 



emission of the FSR single gluon, as long as we stick to very a inclusive cross section, like 



the total cross section in Figure [28|. For stronger angular cuts, or events with definite jet 



multiplicity, we would need to improve our calculation. 

We summarize the results of this section on ISR®FSR as follows: 

• For a typical expt. energy cut of 0.3 ISR®FSR int. is about 1.5% in a tot and A FB . 

• For the energy cut 0.1 it is a factor 2 bigger. 

• The cut | cos#| < 0.9 makes it 25% smaller. 

• The 0{a 1 ) ISR®FSR int. is under total control using KORALZ and KK Monte 
Carlo for arbitrary cuts. 

• Effects beyond 0{a 1 ) are neglig ible, (<20% of 0(a 1 )), except when the energy cut 
is stronger than 0.1. 

• ISR® FSR int. at the Z radiative return is very small, as expected. 

• Changing from s' to Q 2 -propagator in the energy cut has no effect. 

6.6 Total theoretical precision 

Let us summarize the total theoretical precision: 

• For the most typical cut-off range 0.1 < w max < 0.3 excluding the Z radiative return 
we quote for CEEX the total precision 0.2% for LEP2 and for the LC at 0.5TeV. 

• For a cut-off including ZRR we quote 0.2% total precision for LEP2 and 0.4% total 
precision for the LC at 0.5TeV. 

• For 77* we quote 0.3% at LEP2 (no firm result for LC). 

In the above estimates the technical component was significantly below the physical one. 
Restrictions apply: No light fermion pairs (pure photonic QED), no EW component. 
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7 Outlook and summary 



The most important new features in the present CEEX are the ISR-FSR interference, the 
second-order subleading corrections, and the exact matrix element for two hard photons. 
This makes CEEX already a unique source of SM predictions for the LEP2 physics pro- 
gram and for the LC physics program. Note that for these the electroweak correction 
library has to be reexamined at LC energies. The most important omission in the present 
version is the lack of neutrino and electron channels. Let us stress that the present pro- 
gram is an excellent starting platform for the construction of the second-order Bhabha 
MC generator based on CEEX exponentiation. We hope to be able to include the Bhabha 
and neutrino channels soon, possibly in the next version. The other important directions 
for the development are the inclusion of the exact matrix element for three hard photons, 
together with virtual corrections up to 0(a 3 L 3 ) and the emission of the light fermion 
pairs. The inclusion of the W + W~ and ti final states is still in a farther perspective. 
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8 Appendix A: Basic KS/GPS spinors and photon 
polarizations 

The arbitrary massless spinor u\(p) of momentum p and chirality A is defined according 
to KS methods fl24], |25[| . fn the following we follow closely the notation of ref. || (in 
particular we also use C = CO- m the above framework every spinor is transformed out 
of the two constant basic spinors u A (£), of opposite chirality A = ±, as follows 

u x (p) = -^L= piu A (C), u + (C)=?MC), r? = -l, (0 = 0. (209) 

The usual relations hold: (u x (() = 0, u; A u A (C) = u A (C), u A (C)u A (C) =j£wx, &u\(p) = 0, 
u\U\(p) = u\(p), U\(p)u\(p) =pto\, where u x = KI+A75). Spinors for the massive particle 
with four-momentum p (with p 2 = m 2 ) and spin projection A/2 are defined similarly 

u(p, A) = - ? =L= (i> + m) u_ A (C), v(p, A) = -=L= (jrf - m) u x ((). (210) 
V 2p • C V 2p ■ C 

or, equivalently, in terms of massless spinors 

u(p,\) = u x (p ( ) + -j=u^ x ((), u(p,A) = u_ A (p f )--=u A (C), (211) 

where p^ = p = p — ( m 2 / (2(p) is the light-cone projection (p^ = 0) of the p obtained 
with the help of the constant auxiliary vector (. 

The above definition is supplemented in ref. with the precise prescription on spin 
quantization axes, translation from spin amplitudes to density matrices (also in vector 
notation) and the methodology of connecting production and decay for unstable fermions. 
We collectively call these rules global positioning of spin (GPS). Thanks to these we are 
able to easily introduce polarizations for beams and implement polarization effects for 
final fermion decays (of r-leptons, t-quarks), for the first time also in the presence of 
emission of many ISR and FSR photons! 

The GPS rules determining the spin quantization frame for u(p, ±) and v(p, ±) of 
eq. (pilj) are summarized as follows: 



(a) In the rest frame of the fermion, take the z-axis along — £. 

(b) Place the x-axis in the plane defined by the z-axis from the previous point and the 
vector fj, in the same half-plane as ff. 

(c) With the y-&xis, complete the right-handed system of coordinates. The rest frame 
defined in this way we call the GPS frame of the particular fermion. 

See ref. for more details. In the following we shall assume that polarization vectors of 
beams and of outgoing fermions are defined in their corresponding GPS frames. 
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The inner product of the two massless spinors is defined as follows 

s+(pi,p 2 ) = u+(Pi) u -(P2), s_(pi,p 2 ) = m_(pi)m+(p 2 ) = -(s+(Pi,P2))*- (212) 
The above inner product can be evaluated using the Kleiss-Stirling expression 

s+(p, g) = 2 (2K)" 1/2 (2gC)" V2 [(K)(w) - (pv)(qC) ~ * wC W<f ] (213) 

in any reference frame. In particular, in the laboratory frame we typically use ( = 
(1, 1, 0, 0) and rj = (0, 0, 1, 0), which leads to the following "massless" inner product 

s+(p,q) = ~{q 2 + iq 3 )V(P° ~ P l )/{<f ~ I 1 ) + (P 2 + *P 3 )V(<1 ~ ^/(P ~ P 1 )- (214) 
Equation ( |211| ) immediately provides us also with the inner product for massive spinors 



ufa, Xi)u(p 2 , A 2 ) = S(p 1 ,m 1 , \i,p 2 ,m 2 , A 2 ), 
«(pi, Ai)v (pa, A 2 ) = S(pi,m 1 ,X 1 ,p 2 ,-m 2 ,-X 2 ), 
v(pi, Xi)u(p 2 , A 2 ) = ^(pi, -mi, -Ai,p 2 ,m 2 , A 2 ), 
v(pi, Xi)v(p 2 , A 2 ) = S(pi, -mi, -Ai,p 2 , -m2, -A 2 ) 



(215) 



where 



S(p 1 ,mi, Xi,p 2 ,m 2 , X 2 ) = 5 Xl ,-x 2 sx 1 {pi P2 <: ) + S Xu x 2 [ mi \/^^ + m2 y I ' ( 216 ) 

In our spinor algebra we shall exploit the completeness relations 

p' + m = ^2u(p,X)u(p,X), $-171 = ^2v(p, X)v(p, A), 

A A 

1/t = ^^u(k, X)u(k, A), k 2 = 0. 



(217) 



For a circularly polarized photon with four- momentum k and helicity a = ±1 we adopt 
the KS choice (see also ref. [f)8|D of polarization vector^] 

V2 u^ a (k)u a (f3) V2 m_ ct (/c)u ct (C) 

where /3 is an arbitrary light-like four- vector (3 2 = 0. The second choice with u a (Q (not 
exploited in P^fl ) often leads to simplifications in the resulting photon emission amplitudes. 
Using the Chisholm identity^] 

u a {k)^u a {(3) r = 2u a {(3) u a {k) + 2u^{k) u_ a (P), (219) 



26 Contrary to other papers on Weyl spinor techniques |39| we keep here the explicitly complex 
conjugation in e. This conjugation is cancelled by another conjugation following from Feynman rules, 
but only for outgoing photons, not for a beam photon, as in the Compton process, see ref. fi0|| . 

27 For (3 = £ the identity is slightly different because of the additional minus sign in the "line-reversal" 
rule, i.e. M (T (fc)7 M u (T (C) = —U- CT (C)^ ,J 'U- a (k), in contrast to the usual u a (k) r y ll u a (l3) = +u_ cr (/3)7 A 'u_ (T (fc). 
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Ua{k)-ff,u a {C) Y = 2u CT (C) u a (k) - 2u^ a {k) u_ CT (C), 
we get two useful expressions, equivalent to eq. (|218|) : 

V2 



(220) 



u- a (k)u a (f3) 



[u a {(5)u a (k) + u-. a (k)u- a ((3)] 



V2 



(221) 



VKk 



[u ff (C)u«r(fc) - u_ CT (A;)u_ CT (C)] 



In the evaluation of photon emission spin amplitudes we shall use the following im- 
portant building block - the elements of the "transition matrices" U and V defined as 
follows 



tifa, AO «(P2, A 2 ) = 17 © Ea 2 J = U^ X2 (k, Pl ,m 1 ,p 2 ,m 2 ) 

v( P i,Xi) C(k,C)v(p 2 ,X 2 ) =V{t)[ZZ] = VZ lM {k, Pu m uP2 ,m 2 ). 

In the case of u a (() the above transition matrices are rather simple^: 



(222) 



U + (k,p 1 ,m 1 ,p 2 ,m 2 ) = V2 



U Xli x 2 (k,p u m 1 ,p 2 ,m 2 ) = [-U^ Xi (k,p 2 ,m 2 ,p u m 1 )]* , 
^,A a (*>Pi» m i»P2,m 2 ) = Ul Xi ^ X2 (k,p u -m!,p 2 , -m 2 ). 
The more general case with u a (j3) looks a little bit more complicated: 

v/2 



(223) 

(224) 
(225) 



U + (k,p 1 ,m 1 ,p 2 ,m 2 ) 



s-(k,/3) 



X 



s+(pi, fc)s_(/3,p 2 ) + m i m 2^/^"|J, m 1 ^f i s + (k,p 2 ) + m 2 ^^s+{p 1 , k) 



(226) 



with the same relations ( |224j ) and ( 225|) . 

In the above the following numbering of elements in matrices U and V was adopted 



{(Ai,A 2 )} 



-+ 



(227) 



When analysing (multi-) bremsstrahlung amplitudes we shall also often employ the fol- 
lowing compact notation 



U HSU = KM^P^P2,m 2 ), V [ p Xi t P X2 \ = ^ iA2 (fc,p limi)P2) m 2 ) 



(228) 



28 Our U and V matrices are not the same as the M-matrices of ref. [£5) , but rather products of several 
of those. 
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When analysing the soft real photon limit we shall exploit the following important 
diagonality property^ 

U [V*l] =^E,"J =h(k,p)6 XlX „ (229) 

(230) 

which also holds in the general case of u a (j3), where 

y/2 ( ,„ m 2 



MM = ^ p) (s_ a ((3,P>a(p,k) + — y/(20C) (2C*0 ) ■ (231) 



29 Let us also keep in mind the relation b-„(k,p) = — (b a (k,p))* , which can save time in the numerical 
calculations. 
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